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Introduction 



In the thirties, W. Pauli put forward the idea of a new neutral particle that appeared together with 
the electron in the final state of nuclear beta decay. In its origin, this was "a desperate solution" 
to preserve the idea of energy conservation in physical processes. Apparently energy disappeared 
in beta decay reactions without leaving any trace. In fact, that was the case since energy was 
carried away by the invisible neutrino. More than twenty years elapsed, till C. Cowan and F. 
Reines observed for the first time electron antincutrinos coming from a nuclear reactor. Despite 
the multiple applications neutrinos had in Particle Physics through these years, we still know very 
little about these elusive particles. Almost a decade ago, we had strong experimental evidence that 
atmospheric neutrinos do oscillate (they change their leptonic flavour) in their path to the detector. 
A similar mechanism is at the heart of the so-known solar neutrino problem. This set of results 
was beautifully confirmed through experiments done with man-made neutrinos (those coming from 
reactors and accelerators). We are now convinced that neutrinos are massive particles, however 
with oscillation experiments, we cannot tell anything about absolute neutrino masses. There are 
other experiments that carefully measure the kinematics of reactions were neutrinos are involved 
(like the tritium beta decay) that are able to restrict the absolute mass of the electron neutrino 
below 2 eV. Although massive, neutrinos have a very tiny mass compared to their corresponding 
charged leptons. 

The evidence for neutrino masses gave rise to a true revolution in Particle Physics, since this 
was the first solid clue about the existence of new Physics beyond the Standard Model. Within 
this model, the neutrino was exclusively left-handed and therefore could not be given a Dirac mass 
term. Experimentally we have never observed right-handed neutrinos, but there is no fundamental 
principle that forbids those neutrino states. In fact, many extensions of the Standard Model do 
contemplate their existence. 

After a brief introduction of the Standard Model in Chapter [TJ we will see how to build mass 
terms for neutrinos. We will deal with fundamental questions like the Dirac or Majorana nature of 
neutrinos, the hierarchy problem, the see-saw mechanism or with the possible existence of heavy 
neutrinos at the TeV scale. With rather generality, we introduce the neutrino mixing matrix 
that allows us to describe the results obtained in neutrino oscillation experiments. Since many 
theoretical models propose the existence of heavy neutrinos, weakly mixed with the light ones, we 
perform a phenomenological study in Sec. ll.4l on the precision that next-generation lepton colliders 
will have in establishing limits to such an admixture and on the heavy neutrino masses. This work 
is thoroughly described in Ref. [l]. 

The rest of the thesis is devoted to assess the potentialities of Liquid Argon (LAr) TPCs for 
the study of neutrino properties. Nowadays there is a tremendous experimental effort going on 
in order to precisely measure the mixing angles, neutrino masses and the existence of a possible 
CP violation in the leptonic sector. In this context, LAr TPCs will allow us to study neutrino 
interactions with matter with unprecedented precisions, since the detector is fully sensitive and 
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acts as a fine-grained track detector, recording very precisely the deposited energies of particles 
that propagate through argon, making thus possible an accurate reconstruction of a large variety 
of events. Moreover, this device records events with an imaging quality similar to those of ancient 
bubble chambers (see Fig. [T]), but with the additional advantage of providing fully- automatized 
three dimensional reconstruction of the events. 




Figure 1: (Left) First observation of a neutrino hitting a proton in a bubble chamber. The collision 
occurred at the point where three tracks emanate on the right of the photograph. (Right) Another 
neutrino-nuclcon interaction, this time measured by one of the read-out planes of a Liquid Argon 
TPC. 

The LAr TPC technology is discussed in Chapter where we discuss all the relevant technical 
details. One of the major challenges we faced during this work was the development of a set of tools 
and procedures that allowed us to perform a three dimensional and calorimetric reconstruction of 
the recorded neutrino interactions. The goal was to exploit the detector capabilities in order to 
obtain an extremely accurate reconstruction of the neutrino event kinematics. To this purpose, we 
have worked with a small LAr TPC prototype having an active volume of 50 liters (see Chapter|3]for 
further details). In Chapter[3]we discuss the reconstruction tools. Thanks to them, in Chapter[5]we 
have carried out a careful study of quasi-elastic muon neutrino interactions recorded in our detector. 
Those multi-GeV neutrinos were produced in the CERN WANF neutrino facility. Despite the small 
size of our prototype, we have carried out a full study of the i/^n quasi-elastic reactions. 

Those reactions were extensively studied during the sixties and the seventies in experiments with 
bubble chambers, however we still lack a precise modeling of the behavior of nuclear matter after 
a neutrino interacts with it. A systematic study of neutrino interactions with LAr TPCs will allow 
to improve sensibly the Monte-Carlo simulations of nuclear effects. The study of quasi-elastic 
interactions, published in Ref. [5], has been completed with the measurement of the charged- 
current quasi-elastic total cross section for multi-GeV neutrinos. This measurement, together with 
a careful evaluation of the associated systematic errors is described in Sec. 15.41 and has given rise 
to the article described in Ref. [3] . 

The thesis also includes, in Chapter [6l a study of the capabilities a LAr TPC would have when 
searching for direct v-r appearance in — Vr oscillations by means of kinematical criteria. We 
carefully analyze the performance of several statistical techniques for pattern recognition, applied 
to look for very few Vr events in a beam mostly composed of i^^. The study has been published in 
Ref. [J. 
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Corrian los anos 30 cuando por vez primera W. Pauli propuso la existencia de una nueva particu- 
la neutra que acompanaba al electron en los procesos de desintegracion beta. En su origen, esto 
no fue mas que una "solucion descsperada" por salvaguardar el sagrado principio de la conserva- 
cion de la energia: parecia ser entonces que en determinados procesos de desintegracion nuclear 
la energia desaparecia, se esfumaba sin dejar rastro. Efectivamente asi ocurn'a, se la llevaba el 
invisible neutrino. Asi las cosas, todavfa tuvieron que pasar mas de 20 aiios hasta que C. Cowan 
y F. Reines consiguieran observar por primera vez neutrinos procedentes de un reactor nuclear. A 
pesar del tiempo transcurrido desde su descubrimiento, hoy en dfa los neutrinos continuan siendo 
un misterio. Hace aproximadamente una decada se produjo la confirmacion experimental de que 
los neutrinos procedentes del sol y la atmosfera terrestre oscilan o cambian de sabor leptonico en el 
transcurso de su camino hacia el detector. En los afios sucesivos, la existencia de estas oscilaciones 
pudo ser corroborada por experimentos de ultima generacion rcalizados con neutrinos producidos 
artificialmente en reactores y aceleradores. Aunque los experimentos de oscilaciones implican nece- 
sariamente que deben existir distintos estados de neutrino cuya diferencia de masas sea finita, nada 
dicen acerca del valor absoluto de las mismas. Existen, sin embargo, otros experimentos basados 
en el estudio cinematico de reacciones de neutrinos (como la desintegracion beta del tritio) que 
restringen la escala de masas por debajo de 2 eV. La conclusion de estos resultados experimentales 
es clara: Los neutrinos tienen, aunque muy pequefia, una masa no nula. 

El descubrimiento de que los neutrinos son masivos supuso una autentica revolucion en la 
Ffsica de Particulas, ya que representaba la primera evidencia solida acerca de la existencia de 
nueva Fisica mas alia del Modelo Estandar. En cL el neutrino es una parti'cula sin las componentes 
quirales diestras necesarias para dotar de masa a cualquier fermion de Dirac. Esto ultimo se infiere 
del hecho de que, dentro de los h'mites de error experimentales, nunca se han observardo neutrinos 
con helicidad positiva, aunque en realidad no existe ningiin principio fi'sico fundamental que lo 
prohiba. De hecho, existen numerosas hipotesis que dan lugar a modelos teoricos que extienden el 
Modelo Estandar e incorporan la posibilidad de tener neutrinos con masas distintas de cero. 

Tras un breve rcpaso del Modelo Estandar, en cl Capitulo[T] veremos como se pucdcn construir 
terminos de masa para los neutrinos. Cuestioncs como el caracter de Dirac o de Majorana de 
los neutrinos, el problema de las jerarquias y el mecanismo del balancm, o la posible existencia 
de neutrinos pesados a la escala del TeV, son tratadas en este Capi'tulo. De una forma efectiva 
y completamente general introducimos la matriz de mezcla de los neutrinos, que en el marco 
del Modelo Estandar minimal con neutrinos masivos y conjuntamente con las masas de estos, 
permite describir los fenomenos de oscilaciones de neutrinos. Puesto que muchos de los modelos 
teoricos proponen en sus solucioncs la existencia de neutrinos pesados, dcbilmcnte mezclados con 
los ligcros, cn la Seccion 11.41 realizamos un estudio fenomcnologico sobre la precision con la que 
los colisonadores de leptones de proxima generacion podran establecer los h'mites de csa mezcla y 
de las masas dc los neutrinos pesados. Este trabajo se dctalla en la publicacion dc la Rcf. [T]. 
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El rcsto de la tesis esta dedicado a mostrar el potencial que poseen las camaras de proyeccion 
temporal (TPC) de Argon Li'quido (LAr) para el estudio de las propiedades de los neutrinos. Hoy 
en dia se esta realizando un enorme esfuerzo experimental para determinar con precision las masas 
de los neutrinos, sus angulos de mezcla y una posible violacion de CP en el sistema leptonico. 
En este contexto, las TPCs de Argon Lfquido van a permitir estudiar las interacciones de los 
neutrinos con la materia con precisiones nunca antes alcanzadas. Todo el volumen sensible de este 
dispositivo actiia como un detector de trazas y de energi'a depositada con una alta precision y 
granularidad, haciendo asi posible una precisa reconstruccion cinematica de una gran variedad de 
sucesos. Ademas, su dispositivo dc lectura es capaz de registrar los sucesos en distintas oricntaciones 
con una calidad de visualizacion similar a la de las antiguas camaras de burbujas (ver Fig. [5]). pcro 
con la ventaja adicional de hacer posible una reconstruccion tridimensional. 




Figura 2: (Izquierda) Primera observacion de un neutrino golpeando a un proton en una camara 
de burbujas. La colision ocurrc a la derecha de la imagen, donde de repente, surgen tres trazas (un 
muon, un proton y un pion). (Derecha) Otra interaccion neutrino- nucleon, esta vcz recogida por 
uno dc los pianos de deteccion de una camara dc deriva dc Argon Li'quido. 

La tecnologi'a del Argon Lfquido se presenta en el Capi'tulo[2l donde se exponen los principios 
de la tccnica dc deteccion, las capacidades de visualizacion y algunas consideraciones tccnicas de 
interes. Uno dc los mayores rctos de nuestra investigacion consisti'a en el desarroUo de un conjunto 
de herramicntas y proccdimicntos que hicieran posible la reconstruccion calorimetrica y tridimen- 
sional dc los sucesos ionizantcs que sc producen en el detector. El objctivo, cxplotar las capacidades 
del detector para lograr una mcdida sin precedentes de la cinematica de las interacciones de neutri- 
nos. Para tal cfccto, se ha trabajado a fondo con un pcquefio prototipo dc Argon Lfquido con unos 
50 litros dc volumen activo, descrito en detalle en el CapftuloO El Capftulo S] esta dedicado a la 
descripicion de las herramientas de reconstruccion, asf como a la discusion de su rendimiento y po- 
sibilidades. Gracias a los progresos obtenidos en los algoritmos de reconstruccion, en el Capftulo [S] 
hemos realizado un estudio de las interacciones quasi-elasticas de neutrinos del muon, procedentes 
de un acelerador de partfculas a energfas de varias decenas de GeV, registradas en el interior de 
la TPC dc 50 litros. Pcse a su limitado tamano, el conjunto de sucesos recogidos ha permitido 
un analisis completo de la cinematica dc las interacciones quasi-elasticas del tipo J^^n —f fJ.~p. 
Aunque este tipo de reaccioncs fueron cxtensamcnte estudiadas en los experimentos con camaras 
de burbujas en los anos 60 y 70, todavfa no cxiste una modelizacion precisa dc los mecanismos 
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de interaccion del neutrino con la materia nuclear. El estudio sistematico de las interacciones de 
neutrinos mediante los detectorcs de Argon Liquido permitirfa una mejora substancial de las simu- 
laciones Monte-Carlo, que tan importantes son siempre en la experimentacion con neutrinos. Este 
trabajo, publicado en la Ref. [2], ha sido completado posteriormente con la medida de la seccion 
eficaz de colision quasi-elastica, via corrientes cargadas, para neutrinos multi-GeV dispersados por 
nucleos de Argon. Dicha medicion, junto con la evaluacion de sus errores sistematicos esta expuesta 
en la Seccion [5^ de esta Tesis, y puede consultarse en la Ref. [3]. 

La memoria de la Tesis concluye en el CapitulolHl dedicado al estudio de las capacidades que 
una camara de Argon Liquido tendria a la hora de realizar una biisqueda directa de la oscilacion 
v^-Vt mediante criterios cinematicos. Allf tambien analizamos el rendimiento de distintas tecnicas 
estadisticas para el reconocimiento de seiiales, aplicadas al problema de biisqueda de Vr en un haz 
compuesto mayoritariamente de v^. Toda la discusion al respecto se encuentra publicada en la 
Ref. H. 



Chapter 1 

Neutrino interactions 



The main goal of this Thesis is to prove the optimal behaviour of Liquid Argon for detection of 
light neutrinos, with the aim of studying neutrino oscillations. Let us then first in this Section 
briefly review what we know about massive neutrinos. Thus, after introducing the Standard Model 
(SM) of Particle Physics, we revise the two types, Dirac and Majorana, of neutrino masses and the 
see-saw mechanism for three generations. Then, we derive the oscillation probabilities for flavor 
transitions between light neutrinos, and discuss the limits on the mixing between light and heavy 
neutrinos. 

1.1 Overview of the Standard Model 

Our present understanding of the elementary particles and their interactions is described by the 
so called Standard Model [S1[H1[3[H]. It contains: 

• Fermions with spin 1/2, which are the particles of which matter is made. 

• Gauge Bosons with spin 1, which act as mediators of the different interactions. 

• Scalar s with spin 0, which are responsible of the spontaneous breaking of the electroweak 
symmetry. 

The fermions are either leptons which are not affected by the strong interactions or quarks which 
arc the fundamental constituents of hadrons, the particles interacting via the strong force. Leptons 
and quarks can be grouped into three generations according to their masses as shown in Tab. 11.11 
Ordinary matter consists of fermions of the first generation. The higher mass fermions of the second 
and third generation are produced by particle accelerators and cosmic rays at higher energies. 
For every fermion there is an associated anti-fermion with the same mass but opposite additive 
quantum numbers. Fermions interact with each other via the exchange of the gauge bosons listed 
in Tab. 11.21 The photon is the gauge boson which mediates the electromagnetic interactions and 
couples to all charged particles. All fermions interact weakly exchanging the massive gauge bosons 
and Z*^. The strong interaction acts only on quarks and is carried by gluons. To complete the 
list of the known forces, gravity has to be mentioned. It affects all particles and is assumed to be 
carried by the graviton, a spin-2 boson. However, since the effects of gravity are negligible when 
considering the interactions of elementary particles, it will not be further discussed. 

The interactions of the Standard Model are based on the principle of local gauge invariance, 
which demands the Lagrangian of the theory to be invariant under transformations of the gauge 
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Quarks 




r 1 fl'\7TmT' 

i ICL V \J Hi. 




Mass (■MrV/r2l 


1st 


u 


+2/3 


1.5 - 4 




d 


-1/3 


4-8 


2nd 


c 


+2/3 


1150- 1350 




s 


-1/3 


80 - 130 


3rd 


t 


+2/3 


171400 ±2100 




b 


-1/3 


4100 - 4400 



Lcptons 


Generation 


Flavour 


Charge (e) 


Mass (McV/c2) 


1st 







< 3 X 10-« 






-1 


0.511 


2nd 







< 0.19 






-1 


105.7 


3rd 







< 18.2 




r~ 


-1 


1777 



Table 1.1: Properties of the SM fermions [5]. Since quarks do not exist as free particles, their 
masses cannot be measured directly, but have to be derived from hadron properties. Therefore, 
the values of the quark masses depend on how they are defined. The constituent quark masses 
are inferred from the hadron mass spectra, while the masses that enter as parameters in the QCD 
Lagrangian are called current masses. The values stated in the Table refer to the latter. 



group 

SUi3)c X SU{2)l X U{1)y ■ (1.1) 

QCD is governed by the symmetry group SU{3)c' and the electroweak theory by SU{2)l x C/(l)y. 
Every fermion in the SM belongs to a representation of the gauge group which uniquely determines 
its dynamics. The different representations where the fermions lay are summarized in Tab. 11.31 
The non-Abelian character of SU{3)c and SU{2)l leads to self-interactions between the gauge 
bosons. In the case of QCD, this self-interaction accounts for the fact that quarks are confined into 
hadrons, i. e., they cannot be observed as individual particles. However, probed at higher energies 
they begin to behave as free particles, a property referred to as asymptotic freedom. 

Whereas the gauge bosons of the electromagnetic and the strong interaction arc massless, the 
theory of the weak interactions must also accommodate the massive and gauge boons. 
In order to achieve this, while retaining the gauge invariant structure of the theory, the Higgs 
mechanism was introduced. This, however, requires at least one additional physical particle, the 
Higgs boson, yet to be found. 

Although the SM is a mathematically consistent, renormalizable field theory, it cannot be 
the final answer [11| . Even if the recently discovered neutrino masses were incorporated, more 
fundamental shortcomings can be noted. For example, the minimal version of the model has more 
than 20 parameters, which, as most physicists believe, are too many for a fundamental theory. 
Although incorporated in the SM, there is no explanation for the fundamental fact of charge 
quantization or for the CP (charge-conjugation and parity) violating interactions. There is no 
profound compelling for the existence of the heavier fermion generations either. Moreover, there is 
no prediction of the charged fermion masses which vary by more than five orders of magnitude. As 
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vzrcLUgt, 


XXClcLbl 








boson(s) 


strength 




fermions 


strong 


gluons (.g) 


1 


10-^5 


quarks 


electromagnetic 


photons (7) 


10-2 


00 


aU charged 


weak 




10-7 


10-18 


aU 


gravitational 


[graviton (G)] 


10-39 


00 


aU 



Table 1.2: SM interactions and corresponding gauge bosons. The relative strengths of the forces 
are roughly given for short distance scales of a few GeV [TUI. Since, it has not been observed, the 
graviton may be considered a hypothetical particle for the time being. 



Field (particle) 


SU{3)c 


SU{2)l 


U{1)y 




[ 


'A 1 


:-). 


1 


2 


-1/2 






Mi? 


TR 


1 


1 


-1 


(;) 


r. 


t). 


(a 


3 


2 


1/6 






CR 


tR 


3 


1 


2/3 


d'R 




4 


b'R 


3 


1 


-1/3 



Table 1.3: Representation assignment of the SM fermions. 



accelerators reach higher energies and results become more precise answers for some of the current 
questions will be certainly found. 

Further details on the structure of the electroweak interactions can be found in Appendix |^ 
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1.2 Massive neutrinos 

Neutrinos introduce a hierarchy problem: why are they so much hghter than their charged partners? 
In the SM this hierarchy is made natural making the neutrino massless. Recent experimental results 
strongly indicate the existence of oscillations for neutrinos produced in the atmosphere [H [131 Ej 
and those coming from the sun [15|, 1161 117) . which can be interpreted introducing finite, non- 
degenerate neutrino masses. The possibility of considering neutrinos as massive particles allows 
for new phenomena such as the already mentioned flavor mixing and oscillations, but also magnetic 
moments and decays, with very important consequences in astrophysics, cosmology and particle 
physics. Moreover, since most extensions of the SM predict the existence of massive neutrinos, this 
should be considered as a potentially rich field to look for physics beyond the SM. 

In this section we briefiy review the formulation of neutrino masses, extending the minimal SM 

[HIiniEolET]. 
Dirac mass 

If as the other SM fermions neutrinos have right-handed counterparts in SU{2)l singlets, the 
standard Higgs mechanism can generate the corresponding Dirac mass term 

= -rnoi^i^ , (1.2) 

where we assume for simplicity only one family to start with. In terms of the chiral fields 
(left-handed neutrino) and i^r (right-handed neutrino) 

£° = -TTlD (i^L '^R + I'R vl) , (1-3) 

mixing fields of opposite chirality. mD = yv/^/2 where y is a dimensionless Yukawa coupling and 
is the vacuum expectation value of the Higgs field. This mass term is invariant under the 
global transformation v — > e'" which implies the conservation of a global quantum number, in 
this case the lepton number L. 

All experiments to date are consistent with the neutrino being left-handed and the anti-neutrino 
being right-handed. Whether the two other states exist in Nature is at present unknown. They 
would not take part in the weak interactions and for this reason are named sterile. The main 
theoretical problem associated with this mechanism is the smallness of the neutrino mass compared 
to the mass of the charged leptons. Since particle masses are directly determined by the particle 
couplings to the Higgs field, this implies couplings differing by several orders of magnitude for 
particles in the same generation. The relation between lepton masses of different generations may 
be thought of as being "accidental" , but a large mass difference between particles belonging to the 
same SU{2)l doublet is more difficult to argue. However, light neutrinos can have also masses in 
the absence of their right-handed counterparts. 

Major ana mass 

In 1937 Majorana |H] discovered that a massive neutral fermion as a neutrino can be described 
by a spinor with only two independent components, imposing the so-called Majorana condition 

1/ = (1.4) 

where v'^ = CV^ = C^'^'^ v* is the operation of charge conjugation, and C the charge-conjugation 
matrix (see |19j). The Majorana condition (|1.4p can be conveniently expressed in terms of its 
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left-handed and right-handed components 

VR^vl , jy^^iy'^, (1.5) 

Thus, the right-handed component vj^ of the Majorana neutrino field v is not independent but 
obtained from the left-handed component i/^ through charge conjugatior0, and the Majorana field 
can be written 

iy = iyL + i^l. (1.6) 
Using the constraint p.Sp in the mass term p.3|) . we obtain the Majorana mass term 

^-\mM(^L^L + i7Ivl) , (1.7) 

where we have inserted the factor 1 /2 to properly normalized the kinetic term in the Lagrangian. 
Hence, a mass term for purely left-handed particles (and right-handed anti-particles) can be con- 
structed (|1.7p . Thus, in contrast to the Dirac case a neutrino mass term is possible even without 
additional (unobserved) right-handed neutrinos. It requires just one helicity state of the particle 
and the opposite helicity state of the anti-particle. 

The Majorana mass term is not invariant under global transformations. Indeed, a lepton 
number violation of Ai = ±2 is produced when considering a Majorana mass term. In fact, a 
renormalizable theory which includes a Majorana mass term for can not be constructed using 
the Higgs sector of the SM. In general, a Majorana mass for the neutrino would be interpreted as 
a hint for new physics at higher energies, whereas a Dirac mass could be accommodated adding to 
the SM the right-handed neutrino states but requiring very small neutrino Yukawa couplings. 



General mass terms and the see-saw mechanism 

The same reasoning as for can be applied to a hypothetical (in principle) and independent 
right-handed Majorana neutrino, which can be equivalently written n ~ + n'j^ d with a mass 
term ^ 

-Cr = --TOR(rr^ri5j + n^ni?) . (1.8) 

The right-handed Majorana neutrino does not need to be the "partner" of the (ordinary) left- 
handed neutrino, but could be, instead, an independent heavier state. In general, if both exist and 
are independent, in addition to the Dirac mass term (jl.Sp also the Majorana mass terms for 
and n_R, in (|1.7|) and (|1.8|) respectively, are allowed. Then, the general mass term 

can be written as 




(1.10) 



It is clear that the chiral fields vl and n^j do not have a definite mass, since they are coupled by 
the Dirac mass term. In order to find the fields with definite masses, it is necessary to diagonalize 
the mass matrix in Eq. (|1.10p . This is done by a unitary transformation 

" " - (7 ■ (1^") 

^We assume the notation = {u^)'^, which is a right-handed field, contrary to which is a left-handed 

field. 

■^Right-handed sterile neutrinos will be denoted by n in what follows, and by A'^ when discussing collider limits. 
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with TOfe, = 1, 2, real and positive and U the unitary mixing matrix relating the and ur fields 
with the mass cigenstates ui l and 1^2 l ■ This diagonalization leads to the following mixing matrix 
and eigenvalues 



U={ ''""'^ , tan 20= ^"^^ 



-ismt/ cos( 



mi? + "IL 



The corresponding mass eigenstates are given by 



vil \ _ / * cos a —I sm U\ I 
V2l) V sin0 COS0 ) In^ 



An important fact to note is that the diagonalized mass term 



(1.12) 



(1.13) 



^Mass^ 1 ^ ^^^e^^^^^jj.C. (1.14) 

k=l,2 

is a sum of Majorana masses for the Major ana fields 

>^k = VkL + vlL, k = l,2. (1.15) 



There arc two interesting physical limits. In the first case, mo 3> mL,mji. Then. i>i and 
1^2 become almost degenerate with a mass mi. 2 — rnjj and mixing angle 9 ~ 45°. Thus, the 
two Majorana cigenstates form a pseudo-Dirac neutrino state, and the associated mass term in 
the Lagrangian conserves lepton number to a good approximation, vl and vr are the active and 
sterile components (in the SM), respectively, and are maximally mixed. However, as we have 
already mentioned (Sec ll.2p . this solution does not explain why neutrino masses are so small. 

The other more interesting case is the see-saw limit |23j . obtained for niR ^ mr)^ rn^, in which 
case the masses for vi and 1/2 are 

mi ~ WLL, m2 ~ niR , (1-16) 

m_R 

differing enormously from each other. If m^ ~ 0, mi increases with decreasing m2 for a fixed 
value of m£). Furthermore, the mixing angle 9 approaches zero, so that 1^1 and 1^2 arc completely 
decoupled. Usually, mu is assumed to be of the order of the quark and charged lepton masses, 
and m/j is a large mass of the order of some unification scale beyond the SM. This forces the 
mass of the ordinary neutrino to be small, and postulates a very heavy right-handed neutrino, 
therefore decoupled from low energy physics. This see-saw mechanism is by far the most attractive 
explanation of the smallness of the ordinary neutrino masses. 



1.2.1 Neutrino mixing 

The hypothesis of neutrino mixing was first proposed by B. Pontccorvo for the —^ transi- 
tions in analogy with /v" — > K^. After the discovery of the muon neutrino [53], neutrino 
mixing between different neutrino flavors was introduced [26j . This phenomenon appears because 
the neutrino flavor cigenstates i^q, a = e,/i, r, which take part in weak interactions, are generally 
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quantum- mechanical superpositions of difFcrent neutrino mass eigcnstates (see Sec. IA.1|) . which 
describe the propagation of neutrinos in space-time. 

From a large variety of experimental data (see [11|). it is well known that there are just three 
neutrinos that participate in the weak interactions. Hence, let us consider three left-handed chiral 
fields i^eL , i^tiL 1 i^tL that describe the three active flavor neutrinos and three corresponding right- 
handed chiral fields ni/j, n2R, n^R that describe three sterile neutrinos [^1, which are blind to weak 
interactions. The corresponding mass term is given by Eq. ()1.10p with 

where Ml, Mr and Md are 3x3 matrices and 

Vf_,L and nR = n2R . (1-18) 

VrL ) ynsR/ 

The 6x6 matrix A4 can be diagonalized by a unitary matrix U, analogous to the one in Eq. ()1.12p . 
The mass term is then written as 



Mass 



I 6 

^ ^ ?«fc t'fcL J^fcL + H.c. , (1-19) 



k=l 



which is a sum of Majorana mass terms for the Majorana neutrino fields. The mixing relations 
between mass and weak interaction eigenstates read 

6 6 
VaL = '^Uakl'kL, a = e,fJ,,T, naR='^UakVkL, a =1,2, 3. (1-20) 

fe=i fc=i 

This mass term allows the implementation of the two former limits. In the case of M^^r = we have 
three light Dirac neutrinos, and when the Mr eigenvalues are much larger than the Ml and Ml) 
entries, the see-saw mechanism is at work. This case is, however, considerably more complicated 
than the one generation case, briefly discussed above. The see-saw mechanism naturally gives 
a neutrino mass scenario with a light and a heavy neutrino sector, which are to a very good 
approximation decoupled. 

In the first case {Ml,r = 0) the mixing matrix U can be obtained in two steps. One diagonalizes 
Mb first, and afterwards splits the Dirac neutrinos in their respective mass eigenstates (as explained 
above for one flavour). The first step is fulfilled hy a, U matrix block-diagonal, with the (unitary) 
3x3 upper-left part U describing the charged current interaction^ In the second case, with the 
eigenvalues Mr much larger than Ml and M^, the matrix is only to a very good approximation 
block-diagonal. Although in neutrino oscillations, which is the main concern of this work, it can 
be in practice also taken to be exactly block-diagonal in the Majorana case. In Section [L4l we will 
briefly review the present as well as the expected limits on the small off-diagonal entries giving the 
mixing between light and heavy neutrinos. 

For low-energy phenomenology, the upper-left submatrix U is sufficient to describe the ef- 
fective mixing of the active flavor neutrinos This scenario, called "three-neutrino mixing" , can 

■^Let us remark, however, that the number of sterile neutrinos is not constrained by terrestrial experimental data, 
because they cannot be detected, and could well be different from three. 

''When we refer to the block-diagonal structure of U is this case, we do to this first step. 

^This matrix is analogous to the Cabibbo-Kobayashi-Maskawa matrix for quarks (see Appendix lA. 1 . 21 . and it 
is usually called Maki-Nakagawa-Sakata (MSN) matrix 1261 . Hence, it diagonalizes the neutrino mass matrix in the 
current eigenstate basis where the charged leptons have a well-defined mass. 
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accommodate the experimental evidence for neutrino oscillations in solar and atmospheric neutrino 
experiments 0. 

Taking into account the experimental constraints from current oscillation experiments, the 
best- fit value for the mixing matrix U is |35| : 

-0.83 0.56 0.00 \ 
JJbf ~ I 0.40 0.59 0.71 . (1.21) 
0.40 0.59 -0.71/ 

In terms of confidence limits at 90% (3cr), the range of variation on the magnitude of the corre- 
sponding matrix elements are |36| : 

/(0.73) 0.79 - 0.86 (0.88) (0.47)0.50- 0.61 (0.67) - 0.16(0.23) \ 

|;7|= (0.17)0.24- 0.52(0.57) (0.37)0.44- 0.69(0.73) (0.56)0.63- 0.79(0.84) . (1.22) 
\(0.20) 0.26 - 0.52 (0.58) (0.40) 0.47 - 0.71 (0.75) (0.54) 0.60 - 0.77 (0.82) / 

This mixing matrix, with all elements except Uez large, is called "bilarge" . It is very different 
from the quark mixing matrix, which has small off-diagonal entries. Such a difference may be 
an important piece of information for our understanding of the physics beyond the SM, which 
presumably involves some sort of quark-lepton unification. 



1.3 Neutrino oscillations 

Due to the different propagation of the neutrino mass eigenstates, the flavour content of the neu- 
trino propagating state changes with time. Since neutrinos are detected through charged current 
processes, which are sensitive to the neutrino flavor, this oscillating behaviour may be observable. 



Oscillations in Vacuum 

Neglecting the small mixing with the heavy sector, the active neutrinos relevant to oscillation 
experiments can be written 

3 

VaL ^'^UaiViL, a^e,fi,T, (1-23) 
i=l 

where vil, i'2L, f^L are the left-handed components of the light neutrino mass eigenstates. Let us 
assume that a flavour eigenstate \va) was produced at a; = and t = 0, with 



3 
i=l 



The space-time dependence of the free mass eigenstates |t'i(i)) of momentum p and energy Ei = 
\jjp^^\^m^^ reads 

\v.m = e*(^^-~-^'*) \v,) . (1.25) 



°In definite models with very heavy right-handed neutrinos the baryon asymmetry observed in the universe can 
originate from the interactions of these neutrinos, and then from an initial Icpton asymmetry |27| (see for example 
Refs. [m [29] IMl ED 132 [33] for specific models and Ref. [34] for a review). 



1.3. Neutrino oscillations 



15 



It follows from Eq. (jl.23p that the space-time evolution of the state \va) is then given by 



i=l 



u* 



iEit\ 



(1.26) 



If the masses rrii are not equal, the three terms of the sum in Eq. (|1.26p get out of phase and 
the state \i'a{t)) acquires components with [3 ^ a. The probability with which the neutrino 
produced as \va) is converted into ji^^) in a different point of space-time is given by 



-i{Eit-px) 



Up, 



(1.27) 



i>j 

where L is the distance between the production and detection points. E = p the neutrino energy 



for rrii = 0, and Am 



m'^. In the last line we have assumed that neutrinos are highly 



relativistic (to^ ^ p and L = x — t), and that the mixing matrix is real. We see that in order 
to have a non-vanishing oscillation probability, it is necessary that at least two neutrino mass 
eigenstates are non-degenerate, and that the matrix U has non-vanishing off-diagonal elements. 

This probability satisfies some important relations. First of all, it becomes evident that neu- 
trino oscillation experiments are not directly sensitive to the scale of neutrino masses, but to the 
difference of their squares. CPT invariance implies Py^_,y^ = Pyc_,yc . For two (or more) flavour 
mixing, we also have P^^^^^^ ~ P^^^^^ if CP is conserved. These relations are valid when 
neutrinos propagate in vacuum. When interactions with matter affect the neutrino propagation 
they may no longer hold because the medium itself is in general not symmetric under CP and 
CPT. Finally, the unitarity of U guarantees that the total probability adds to 1, 



P. 



1 



P 



(1.28) 



The neutrino oscillation phenomena are described by five free parameters. Only two out of the 



three Am? 



appearing in Eq. (jl.27p are independent, since Am|]^ + Ar 



'32 



= Ar 



^3l■ 



On the other 



hand, if we ignore the possibility of CP violation, the matrix U is real and can be parameterized 
by three independent mixing angles, 612, 0i3 and 023- 

For simplicity, it is customary to consider only the case of mixing between two neutrino families. 
In such a case, the mixing matrix U involves only one real mixing angle 9, and the oscillation 
probability simplifies to 



'■Jl 

AE 



Hence, the oscillation amplitude is given by s\t?{29) and the oscillation length (A) by 

A-kE 2.48 E [GeV] 



A 



Am2 Ato2 [cV^] 



[km] . 



(1.29) 



(1.30) 



Thus, in order to observe neutrino oscillations, A must be less or of the order of the source-detector 
distance. This can be rewritten as Am^ > E/L. Therefore, to search for small sc^uared mass 
differences, the distance L should be large and/or the energy of the neutrinos small. 



16 



Chapter 1. Neutrino interactions 



Back to the general three flavour mixing framework, the observed hierarchy Am|]^ ^ ^'^ii — 
Am32 makes the oscillation phenomena to decouple and the two flavour mixing model to be a 
good approximation in limited regions. Each oscillation probability can be written as the sum of 
a "short" and a "long" component 

p _ pshort I plong /-i 01 N 

with 

P'ZX, = -4C/„iC/^i[/a2t/^2sin2(Am2iA) . 

If we assume the mixing matrix to be nearly diagonal, cosOij 3> sin^y, then P^^^y^ and P^^^^^ 
are dominated by the short component and are sensitive to Am|j^ ~ Amgj, whereas P^,^^^^ is 
dominated by the long component and is sensitive to Am^x- A neutrino oscillation experiment 
can be sensitive to cither the short or to both the short and the long components by selecting the 
appropriate L/E range. Tab. 11.41 shows the typical energies, base-lines and minimum detectable 
Am^ for the usual neutrino sources, namely, nuclear reactors, accelerators, cosmic ray collisions at 
the atmosphere and the sun. 





Reactor 


Accelerator 


Atmospheric 


Solar 


Short base-line 


Long base-line 


E (MeV) 

L (m) 
Am^ (eV^) 


< 10 
lOMO^ 

10-2 


30 - W 
102-10^ 
10-1 


30 - 10'^ 
10^-10^ 
10-4 


10^ 
lO^-lO^ 
10-4 


< 14 
10" 

10-11 



Table 1.4: Neutrino energies, baseline distances and minimum testable mass square difference for 
reactor, accelerator, atmospheric and solar neutrinos. 



Oscillations in Matter 

When neutrinos propagate in matter, a subtle but important effect takes place which alters the way 
in which neutrinos oscillate into one another. In constrast with the oscillation in vacuum, where the 
time evolution of the neutrino state is given by the different propagation of the mass eigenstates, 
the presence of matter affects differently the neutrino flavor eigenstates. While all neutrino species 
have the same neutral interactions with matter, the v^. weak eigenstate has a slightly different 
index of refraction than the v^^-r weak eigenstates, due to its charged current interaction with the 
electrons in the medium. This different index of refraction for v,, alters the time evolution of the 
system from that in vacuum. This phenomenon is known as the MSW effect [351 132 • However, it 
is of little relevance to the experimental discussion in the following Chapters. 



1.4 Limits on heavy neutrino masses and mixings 

The neutrino mixing matrix lA in Sec. ll.2.T] can be safely assumed to be block-diagonal in neutrino 
oscillations, at least given the present experimental precision. For a discussion of possible devia- 
tions from unitarity of the upper-left block U , which might be eventually observable in neutrino 
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propagation experiments, see for instance [40j [41]. The small mixing between light and heavy 
neutrinos can be also probed in other experiments, in particular at large colliders [42j . In the 
following, and for completeness, we first discuss present limits and afterwards the expected bounds 
at the International Linear Collider (ILC) [TJH3], because lepton colliders are more sensitive to 
heavy neutrino signals. Let us first introduce a more convenient notation for the 6x6 mixing 
matrix U, 

U^iZ, I \ , (1.32) 



where U ([/') describes the mixing between the light (heavy) neutrinos and V {V) parameterises 
the light-heavy (heavy-light) neutrino mixing. For V — Q the matrix U is the usual 3x3 unitary 
MNS matrix. 



Indirect limits 

The most stringent constraints on light-heavy neutrino mixing result from tree-level contributions 
to processes involving neutrinos as external states like tt —i- £i? and Z vv^ and from new one- 
loop contributions to processes with only external charged Icptons like /i ^ 67 and Z W 
gnilllllSlllllliZlllSlIllllin]. These processes constrain the quantities 

3 3 

1=1 2=1 

because in the former case we must sum over the external light neutrinos (which are not distin- 
guished) and in the latter the sum is over the heavy neutrinos running in the loop. The first type 
of processes in particular tests universality. A global fit to experimental data gives [101 [H] 

ilee < 0.0054 , < 0.0096 , ilrr < 0.016 , (1.34) 

with a 90% confidence level (CL). These limits do not depend on the heavy neutrino masses and 
arc model- independent to a large extent. They imply that heavy neutrino mixing with the known 
charged Icptons is very small, l^iP < 0.0054, 0.0096, 0.016 for £ = e,fi,T, respectively. The 
bound on ^l^e alone does not guarantee that neutrinoless double beta decay is within experimental 
limits for the range of heavy neutrino masses we are interested in (larger or of the order of the 
electroweak scale) ^T\ . 

The second type of processes, involving flavour changing neutral currents (FCNC), get new 
contributions only at the one loop level when the SM is extended only with neutrino singlets, as 
in our case. These contributions, and hence the bounds, depend on the heavy neutrino masses. In 
the limit rriN. ^ Mw, they imply [48j 

|17eM I < 0.0001, |17e^|<0.01, |rj^r|<0.01. (1.35) 

Except in the case of the first two families, for which experimental constraints on lepton flavour 
violation are rather stringent, these limits arc of a similar size as for the diagonal elements. An 
important difference, however, is that (partial) cancellations may operate among heavy neutrino 
contributions. There can be cancellations with other new physics contributions as well. We are 
interested in determining the ILC discovery potential and the limits on neutrino masses and mixings 
which could be eventually established. Then, we must allow for the largest possible neutrino mixing 
and FCNC, although they require cancellations or fine-tuning. 

'^In the following in this Section the subindices £ and i stand for ui and Ui, respectively, and label the corresponding 
current and mass eigenstates. The equality below is a reflection of the unitarity of U. 
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Direct limits 

In order to discuss the ILC limits we concentrate on the hghtest heavy neutrino Ni, and follow 
the analysis in [T] . We omit the subindex 1 and denote the corresponding mixing matrix elements 
ViN- There are two interesting scenarios: (i) the heavy neutrino only mixes with the electron; 
{ii) it mixes with e and either r, or both. A third less interesting possibility is that the heavy 
neutrino does not mix with the electron. We discuss these three cases in turn. In the first two 
cases the heavy neutrino couples directly to the initial state e~^e~ — > Nv, and one must search 
for it in its main decay mode e^e^ £Wv with W —* jj. In Appendix iBl we provide a more 
detailed discussion of this process. Here, it is enough to note that one has to look for a peak in the 
£jj mass distribution. In Fig. 11.11 (a) we plot the ejj invariant mass m^jj for VeN = 0.073. The 
center of mass energy is taken to be 500 GeV. The solid line corresponds to the SM plus a 300 GeV 
Majorana neutrino, being the dotted line the SM prediction. The width of the peak is due to energy 
smearing applied in our Monte Carlo and not to the intrinsic neutrino width r^r = 0.14 GeV. We 
must emphasize that the total cross sections and the derived limits do not depend on the nature 
of the heavy neutrino, but the angular distributions do depend on its Majorana or Dirac mass 
character. In Fig. 11.11 (b) we show the W ^ ev peak which characterizes the main SM background 




500 



Figure 1.1: Kinematical distributions of the ejj invariant mass (a) and the ev invariant mass (b). 
e+e- -> W+W-. 

With convenient cuts on these two distributions in the different leptonic channels fjjV, £ = 
e,fi,T^ one can derive the corresponding combined limits using the signal reconstruction method 
explained in |43| . The integrated luminosity is assumed to be 345 fb^^, as expected after one 
year of running. In Fig. 11.21 we plot the combined limits on VeN and V^^n or Vtn for the same 
heavy neutrino mass of 300 GeV. At the ILC the sensitivities in the muon and electron channels 
are similar, and both are better than in tWv production. This can be clearly observed in both 
plots: a iiNW coupling has little effect on the limits on VeN, but a coupling with the tau decreases 
the sensitivity, because the decays in the tau channel are harder to observe. The direct limit on 
VeN, V^N obtained here improves the indirect one (the solid line in Fig. II. 2[) only for V^i^r < 0.01. 
However, as it has already been remarked, the latter is not general and can be evaded with 
cancellations among heavy neutrino contributions [T]. 

These limits are rather independent of the heavy neutrino mass. The dependence of the total 
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Figure 1.2: Combined limits obtained at ILC on: Vcn and V^at, for Vtn — (a); VeN and Vtn, 
for Vfj^N = (b). The red areas represent the 90% CL limits if no signal is observed. The white 
areas extend up to present bounds, and correspond to the region where a combined statistical 
significance of 5cr or larger is achieved. The indirect limit from /i — e LFV processes is also shown. 



e'^W^v cross section on mjy can be seen in Fig. 11.31 (a), for VeN ~ 0.073, V^jv — VrN — 0. For 
a heavier N the cross sections are smaller and thus the limits on VeN are worse. However, up to 
mjv = 400 GeV this is compensated by the fact that the SM background also decreases for larger 
rriejj. The limits on VeN are shown in Fig. 11.31 (b) as a function of toat, assuming that the heavy 
neutrino only mixes with the electron. 




Figure 1.3: (a) Cross section for e+e^ e^W^v at ILC for VeN = 0.073 and different values of 
niN- (b) Dependence of the discovery and upper limits on VeN on the heavy neutrino mass. Both 
plots assume mixing only with the electron. 

A heavy neutrino signal in the fi or t channels is observable only if the neutrino also mixes with 
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the electron. Wc now quantify this statement. We consider a heavy neutrino coupling only to the 
muon, with V^n = 0.098, or only to the tau, with VrN = 0.13. The beam polarisations P^- = 0.8, 
Pf.+ = —0.6, opposite to the previous ones, are used to enhance the signal and reduce the SM 
background. The cross sections for the SM and SM plus a heavy neutrino are shown in Tab. 11.51 
for these two cases. For N mixing only with the muon, the statistical significance of the signal is 
S/ \fB = 1.85 for one year of running, and 7.3 years would be necessary to observe a ha deviation. 
If the heavy neutrino only mixes with the r, the statistical significance is only S/\/B = 0.76, in 
which case a luminosity 43 times larger is required to achieve a 5(t evidence. 





N - fi 


N -T 


SM 

SM + iV 


1.10 
1.20 


0.389 
0.414 



Table 1.5: e+e — > P^W^u cross sections (in fb) for a heavy neutrino coupling only to the muon 
(first column, £ — fi) or coupling only to the tau (second column, £ ~ t). 



Chapter 2 



The Liquid Argon Time 
Projection Chamber 

During the 1960s and 1970s, bubble chamber detectors (52 EH [H] led to the discovery of many 
new elementary particles. Filled with a superheated liquid, the bubble chamber creates a track 
of small bubbles when a charged particle crosses the chamber, locally bringing the liquid to boil. 
Bubble chambers produce high spatial resolution pictures of ionizing events, and allow a precise 
investigation of the nuclear processes occurring in the medium. However, in recent times, the 
use of bubble chambers has been superseded by electronic detectors, since bubble chambers are 
difficult to trigger, and experiments with high statistics are not really practicable because of the 
time consuming analysis of bubble chamber pictures. The attempt to merge the superb imaging 
capabilities of traditional bubble chambers and the advantages of electronic read-out in a single 
detector led C. Rubbia to propose the Liquid Argon (LAr) time projection chamber (TPC) in 
1977 [SS]. The detector is essentially a cryostat filled with a liquified noble gas and equipped 
with an electronic readout that measures the ionization charge produced by the passage of charged 
particles. This detector provides three-dimensional imaging and, since the ionization charge is 
proportional to the energy deposition, also acts as a calorimeter of very fine granularity and 
high accuracy. Thus, this device is ideal to study particle interactions and docs not present the 
problems of traditional bubble chambers, since the electronic read-out allows the self-triggering 
and automatic processing and analysis of the events. 

ICARUS is a project, proposed in 1985 [5S1 [57], for the installation of a large LAr TPC in 
the Gran Sasso Laboratory (LNGS), Italy, for the study of neutrino physics and matter stability. 
Several years of R&D with prototypes of increasing size have allowed to overcome the major 
technological problems in the establishment of this technique. 



2.1 Principles of the Time Projection Chamber 

The Time Projection Chamber (TPC) is essentially a three dimensional tracking device, capable 
of providing imaging of an ionizing particle track and a measurement of its specific energy loss, 
dE/dx. It consists of a pair of parallel electrodes immersed in an ionization medium (a gas or 
a liquid) and connected to a high voltage power supply (see Fig. 12. 1[) . producing a homogeneous 
electric field perpendicular to the electrodes. 
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2.1.1 Detection principle 

A particle passing through the gas or the Hquid will ionize the medium along the track and create 
ion-electron pairs (see fig. 12. ip . The applied electric field suppresses the immediate recombination 
of the pairs (important mainly in a liquid medium), pulling the electrons and ions apart. The 
ionization electrons of the particle track will drift along the electric field lines to the anode, which 
is composed of an array of position sensitive sensors. 




Sensor 
Planes 



Figure 2.1: Schematic view of a TPC. A charged particle is passing through an ionization medium 
(e.g. LAr) and the electrons drift along the electric field lines to the sensors. 



Imaging of events 

The sensors located in the anode of the detector consist of a set of planes of multiple parallel 
wires (hereafter called wire planes) (see fig. 12. 2p . We can obtain two-dimensional projections of 
the ionization point, constraining one spatial coordinate from the wire collecting the charge (or 
a current being induced) and the other from the drift time. In order to obtain three-dimensional 
coordinates, at least two wire planes with different wire orientations are necessary; therefore, the 
read-out system must be non-destructive. 

The transparency of a wire plane or grid for electrons drifting along the electric field lines is 
a function of the ratio of the fields Ei and , in front and behind the grid respectively, and of 
the ratio p = 2-Kr jp where r is the wire radius and p is the spacing between wires {pitch). Full 
transparency is reached when the following condition is satisfied |58j : 

This condition has to be balanced against the requirement that the grid has to act as an 
electrostatic shielding between the space in front and behind it. If this condition is satisfied, a 
drifting charge will be sensed by the electrodes behind the grid only when the grid itself is crossed. 
Detailed calculations show that the shielding power, cr, of a grid is approximately given by the 
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Charged panicle trajectory 




Figure 2.2: Schematic view of a LAr TPC with three sensing planes. 



foUowing formula [55] : 

, = JLlnfJL) (2.2) 

where d is the distance between the grid and the next electrode. Note, that shielding factors 
deviating a few percent from 1 are acceptable. 

The drifting electrons, following the electric field lines, can cross a succession of several planes 
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of parallel tensed wires oriented in different directions (see Fig. I2.2p . where the condition (|2.ip is 
verified. The electrostatic separation between the sensing wires allows a very sharp localization 
of the drifting charge, hence a good spatial resolution. Each of the wire planes provides a two- 
dimensional projection of the event image with a common coordinate (the drift time), allowing the 
three-dimensional reconstruction of the event. The basic three-dimensional pixel (or "bubble") size 
is determined by the wire pitch and the resolution in the drift direction. The resolution cannot be 
improved by arbitrarily reducing the spacing of the wires. In addition to the practical difficulties of 
precisely stringing wires at a pitch below 1 mm, there is a fundamental limitation: the electrostatic 
force between the wires is balanced by the mechanical tension, which cannot exceed a critical value. 
This gives the following approximate stability condition [SI Eq. 28.14]: 



where L is the wire length, V[fcF] is the voltage of the sense wire and T[g\ is the tension of the 
wire in grams-weight equivalent. 

While approaching a plane, the electrons induce a current only on the wires nearby; when 
moving away, a current of opposite sign is induced. A detailed analysis of the induced signals for 
the geometrical configuration displayed in Fig. 12.21 has been performed by Gatti et al. [59] . Due 
to the shielding of the adjacent grids, the signals have a limited duration. Only the signals on the 
wires of the first grid are "prompt" , i.e. start as soon as the ionizing particle crosses the chamber, 
and can be used to determine the reference time for the electron drift (^o)- The time to can be 
also determined by the coincident detection of the scintillation light produced in the argon, using 
a suitable system of photomultiplicr tubes. 

Selection of the dielectric medium 

The liquid ionization medium has the advantage of a higher density and lower diffusion compared to 
a gaseous medium. The higher density will permit to have an efRcient target for weakly interacting 
particles in a relatively small volume. The high stopping power of the liquid leads to a short 
enough range of the (low energy) charged particles, to be fully contained in the sensitive volume. 
For this type of events the LAr TPC works, in addition to the imaging, as a calorimeter with 
fine granularity. LAr has a radiation length of 14 cm and a nuclear interaction length of 83.6 cm 
giving good electromagnetic and hadronic calorimctric properties. The high ionization density in 
the liquid gives enough electron-ion pairs per unit length to detect directly the ionization charge 
without any multiplication. Thus, the liquid has the advantage to give a high resolution spatial 
tracking information, and, at the same time, acting as a continuous hadronic and electromagnetic 
calorimeter. A summary of Physical and Chemical properties of argon can be consulted in Tab. 12. II 
A big difference for the signal measurement between gaseous and liquid ionization media is 
given by the fact that there is no charge multiplication near the sensing wires in a liquid. As a 
consequence of the short mean free path of drifting electrons in liquid, they don't gain enough 
energy between the collisions to ionize other atoms, and thus, cannot create an avalanche. The 
small amount of charge (about 20 000 electrons for mips and 3 mm wire pitch) reaching the sensors 
must be read by a low noise preamplifier with a high gain (see Sec. 12.1."^ . 

Calorimetry and particle identification 

Charged particles traversing the LAr sensitive volume produce ionization electrons in a number 
proportional to the energy transferred to the LAr. The ionization electrons drift along the elec- 
tric field lines to the wire planes pushed by the electric field perpendicular to the wire planes. 




(2.3) 



2.1. Principles of the Time Projection Chamber 



25 



Atomic number 


18 


Concentration in air 


0.934% 


Naturally occurring isotopes 


'^«Ar = 0.3365(30) % stable 
38 = 0.0632(5) % stable 
""Ar = 99.6006(30) % stable 


Melting point (101325 Pa) 


83.8058 if (-189.3°C) 


Boiling point (101325 Pa) 


87.293 /■s:(-185.8''C) 


Density at boiling point (101325 Pa) 


1.396 A:g/^ liquid 
5.79 q/e gas 


Liquid heat capacity at boiling point (101325 Pa) 


1.078 kJ/kg/K liquid 


Latent energy of fusion at boiling point (101325 Pa) 


161.0 kJ/kg liquid 


dE / dxmin foi' ^ niip 


2.12 MeV/ cm 


Critical energy (electrons) 


31.7 MeV 


Mean excitation potential 


210 eV 


Energy to produce an electron-ion pair 


23.6 eV 


Radiation length Xq 


14.0 cm 


Moliere radius 


9.28 cm 


Nuclear interaction length 


84.0 cm 


Maximal breakdown strength (depending on purity level) 


1.1 - l.iMV/cm 


Diffusion coefficient (89 K) 


A.8cm^ / s 


Recombination factor for mips (/.i) 


0.6 at 0.5 kV/cm 



Table 2.1: Physical and chemical properties of argon. 



inducing a current on the wires near which they are drifting while approaching the different wire 
planes. Therefore, each wire of the read-out plane records the energy deposited in a segment of the 
ionization track. Thus, the combined spatial and calorimetric reconstructions, exploiting the fine 
granularity and imaging capabilities of the detector, allows the precise measurement of the energy 
loss per crossed distance (dE/dx) and hence the identification of the ionizing particle. 

The complete calorimetric reconstruction of the events requires to take into account the effects 
of the electron-ion recombination and charge attenuation by electron attachment to electronegative 
impurities. Both effects reduce the amount of collected charge with respect to the one produced 
during the ionization. The recombination of ion-electron pairs occurs immediately after the ioniza- 
tion, due to their electrostatic attraction. The recombined electron-ion pairs form excited atomic 
states that release the energy in the form of (detectable) light. The magnitude of the recombination 
effect depends on both the electric field and the ionization density (dE/dx). Charge attenuation is 
due to the attachment of the drift electrons to the electronegative impurities, and depends also on 
the electric field and on the concentration of electronegative impurities present on the LAr volume. 
The precise knowledge of these two phenomena allows their off-line unfolding from the calorimetric 
data. 

If the granularity of the read-out electrodes is high enough, a sampling of the energy loss 
per crossed distance (dE/dx) for several points along the ionization tracks is possible. If, on the 
other hand, the momentum of the particle is independently known, e.g. from its curvature in a 
magnetic field, or from its range (for particles stopping in the LAr sensitive volume only), or from 
multiple scattering, we can determine the mass of the ionizing particle by means of the Bethe- 
Bloch formula (see later section which describes the energy loss per crossed distance (dE/dx), 
and hence identify the ionizing particle. The particle identification capability [60] is illustrated 
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Figure 2.3: (Top) Scatter plot showing the regions populated by fully simulated muons, kaons and 
protons on the < dE/dx > vs. kinetic energy plane. (Bottom) Discrimination between kaon and 
pion events. 
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in Fig. 12. 3( where the dE/dx vs. kinetic energy behaviour of stopping muons. kaons and protons 
(top); and a sample of stopping pions and kaons (bottom), arc compared. 

2.1.2 Technical considerations 

The particularities of the detection technique described in the two previous sections, require special 
consideration of certain technical aspects, namely: 

• A system to provide the time of the event occurrence (io) is needed in order to trigger the 
data acquisition (DAQ) system and to determine the absolute position of the event along the 
drift direction. 

• The electronic read-out system must be capable of detecting the signal produced by a few 
thousand electrons. 

• The level of purity of the LAr must be high enough so that the signal is not strongly degraded 
along the drift path. 

The present solutions provided by the ICARUS collaboration for each of these aspects are summa- 
rized below. 

Trigger and to determination 

The trigger system must be suited for the rare events to be detected underground by ICARUS, 
and must provide the absolute time of the event occurrence. In principle, the drift time difference 
between any two points of the projected tracks allows to reconstruct the geometrical topology of 
the events. However, the absolute position along the drift direction must be known when correcting 
(off-line) the measured charge for the effects of the charge attenuation. A precise knowledge of <o 
is thus required to achieve a good energy resolution. 

As we already mentioned in section 12.1.11 the signal on the wire plane facing the drift volume 
starts as soon as the ionizing particle crosses the detector, and the electron-ion pairs are formed. 
This feature can be used to trigger the DAQ system (hence referred as self-triggering capability) 
and to provide the time tg. An alternative system can be provided by a set of photo- multiplier 
tubes (PMT's) immersed in the LAr. which detect the scintillation light produced in the LAr by 
the passage of an ionizing particle. Scintillation photons are produced basically by two different 
processes [6T| : 

• by the direct excitation of an argon atom followed by an excited molecule formation and 
de-excitation: 

Ar* + Ar Ar^ ^ 2 Ar + 7; (2.4) 

• by the formation of a molecular state through recombination processes between electrons and 
molecular ions: 

Ar+ + Ar Ar J + ^ Ky^ ^ 2 Ar 7 (2.5) 

In both cases the average photon wavelength is A = 128 nm, and the emission occurs within a time 
window which depends on the excited state. The de-excitation proceeds from two different states, 
with decay constants of about 5 ns (fast component) and about 1 /is (slow component), respectively. 
The detection of the scintillation light requires the use of PMT's equipped with suitable windows 
(like MgF2, that transmits UV radiation down to 115 nm), or to provide the PMT cathode with a 
proper wavelength shifting system. This trigger system allows the detection of cosmic ray events 
with energies down to a few hundred keV. 
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Read-out system 

One of the main problems which this detection technique must cope with is the relatively small 
amplitude of the signals. This is due to the absence of charge amplification during the electron 
drift process or near the anode. Typically, 1 mm of a minimum ionizing track delivers less than 
10^ electrons in LAr. The signal is even smaller at low electric fields due to the effects of the 
charge attenuation and the electron-ion recombination. The imaging of ionizing events requires, 
therefore, the use of low noise read-out electronics. 

The present solution for the read-out system is the product of several years of tests and experi- 
ence gained by the ICARUS collaboration on small scale prototypes [Hll [Ml [HH ES] . The readout 
system is structured as a multichannel waveform recorder that stores the charge information col- 
lected by each sense wire during the drift of the electrons at a 2.5 MHz sampling rate. Each wire 
module (16 wires) is equipped with a current integrating amplifier feeding a 10 bit flash ADC that 
samples the (multiplexed) signal with a frequency of 40 MHz. 



Argon purification 

In order to make long electron drift paths possible (of the order of one meter or more) LAr must 
be ultra-pure. An electron in the liquid undergoes about 10^^ molecular collisions per second. 
Hence, impurities with large attachment probability must be kept at very low level, on the order 
or less than 0.1 ppb (ppb=:10~^), to achieve electron lifetimes in excess of a few milliseconds. The 
electro-negative impurities are mainly represented by oxygen, water, carbon dioxide and, in minor 
concentrations, some chlorine and fluorine compounds. If the impurity concentration is constant 
over the whole volume, the charge decreases exponentially with the drift time: 

Q{t) = Q{to)e-i (2.6) 

where r is the mean lifetime of the electrons in argon. The lifetime is directly connected to the 
impurity concentration p by an inverse linear relationship [66j : 

r 1 300 

r /is « ^— 2.7 

Thus, the lifetime provides a direct measurement of the LAr impurity concentration (O2 equiva- 
lent). 

Commercial LAr has a contamination level of the order of a few ppm (ppm=:10~^) of oxygen 
equivalent, this corresponds to a mean lifetime of only 0.3 /^s, absolutely insufficient for a TPC. A 
considerable amount of R&D has been performed by the ICARUS collaboration in order to master 
the argon purification process [571 [51] • LAr can also be contaminated inside the cryostat, by the 
degassing process of those parts (walls, electrodes, cables, etc.) covered by the (hotter) argon gas. 
Therefore, it is necessary not only to purify the LAr before the filling of the cryostat, but also 
to ensure a continuous purification inside the cryostat, forcing the LAr recirculation through a 
dedicated purifying unit. 



2.2 The ICARUS detector 

The ICARUS collaboration is planning to build a large LAr TPC as a neutrino observatory and 
a detector to search for nucleon decays |69j- A possible detector would contain 3000 tons of LAr, 
divided into five 600 ton (T600) modules; each T600 module consists of two 300 ton half-modules 
with a common thermal insulation. The first 300 ton module was successfully tested in 2001 
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at the surface of the Earth in Pavia (Italy) [68]; the first T600 module is being installed at the 
Gran Sasso underground laboratory in Italy. The physics program contains the study of neutrino 
oscillations with the long base line v^-hewa from CERN to the Gran Sasso [55] with a search for 
Vr appearance, measurements of solar, atmospheric and supernova neutrinos and a sensitive search 
for nucleon decays. 




In an extensive R&D program several small scale prototypes have been constructed and oper- 
ated, mainly at CERN. Pavia (Italy) and the Gran Sasso Laboratory (Italy), in order to overcome 
the main technological problems in the establishment of the ICARUS technology, and demonstrate 
many of the physics capabilities of the detector. 

ICARUS T600 is a large cryostat divided into two identical, independent, adjacent reservoirs 
{half-modules), containing more than 300 tons LAr each (see Fig. l2.4p . Both half-modules contain 
an internal detector (composed of two TPC's), the field shaping system, monitors, probes and 
PMT's; they are externally surrounded by a set of thermal insulation layers. 

Outside the detector are located the read-out electronics (on top of the cryostat), and the 
cryogenic plant composed of a liquid nitrogen (LN2) cooling circuit and of a system of purifiers 
needed to achieve the required LAr purity. 

The inner detector structure of each half-module consists of two TPCs separated by a common 
cathode. Each TPC is made of three parallel wire planes, 3 mm apart; the wires of the first plane 
are horizontal and the orientation of the two other planes is ± 60° to the horizontal direction; the 
wire pitch is 3 mm. The third wire plane gives redundancy data for better pattern recognition. 
The total number of wires in each half-module is 53 248. The read-out of the signals induced 
on the TPC wires by the electron drift allows the reconstruction of a full three-dimensional (3D) 
image of the event; the 3D reconstruction of the event is obtained by correlating the signals from 
the different planes at the same drift distance. Table [2?2l summarizes some technical data of the 
ICARUS T600 half-module (T300) detector. 

An uniform electric field perpendicular to the wires is established in the LAr volume of each 

^The active volumes are defined as the parallelepipeds with the side surfaces equal to the wire chamber surface 
and with width equal to twice the distance between the first wire plane (0°) and the cathode. 
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Inner dimensions 


19.6 m X 3.6 m x 3.9 m (length, width, height) 


Inner volume 


275.2 m^ 


Activc0dimensions 


17.95 m X 3.0 TO x 3.16 m (length, width, height) 


Active volume 


170.2 m^ 


Sensing planes 


1st Induction, 2nd Induction, Collection, 3 mm spacing 


Channel number 


53248 wires 


Electron drift length 


1.5 m 


Maximal drift time 


1 ms @ 500 V/cm 



Table 2.2: Some technical data of the ICARUS T600 half-module (T300). 



half-module by means of a HV system, as required to allow and guide the drift of the ionization 
electrons. The system is composed of a cathode plane, parallel to the wire planes, placed at the 
centre of the LAr volume of each half- module at a distance of about 1.5 m from the wires of each 
side. The HV system is completed by field shaping electrodes to guarantee the uniformity of the 
field along the drift direction, and by a HV feedthrough to set the required voltage on the cathode. 
At the nominal voltage of 75 kV, corresponding to an electric field of 5001^/ cm, the maximum 
drift time in LAr is about 1 ms. 

The top side of the cryostat hosts the exit flanges equipped with signal feedthroughs for the 
electrical connection of the wires with the read-out electronics, and for all the internal instrumen- 
tation (PMTs, LAr purity monitors, level and temperature probes, etc.). The electronic chain is 
designed to allow for continuous read-out, digitization and waveform recording of the signals from 
each wire of the TPC. It is composed of three basic units serving 32 channels: 

• The decoupling board receives analog signals from the TPC wires via vacuum-tight feedthrough 
flanges and passes them to the analog board. It also provides biasing of the wires and distri- 
bution of calibration signals. 

• The analog board houses the signal amplifiers, performs 16:1 multiplexing and the data 
conversion (10 bit) at a 40 MHz rate. Thus every channel is sampled at 2.5 MHz. 

• The digital board uses custom programmable chips (two per board) specially developed for 
ICARUS, called DAEDALUS, that implement a hit finding algorithm. Each board receives 
the multiplexed digital data via an external serial-link cable. 

Ionization in LAr is accompanied by scintillation light emission. Detection of this light can 
provide an effective method for absolute time measurement of the event as well as an internal 
trigger signal. A system to detect this LAr scintillation light has been implemented based on large 
surface (Sin.) PMTs directly immersed in the LAr. 

The calorimetric measurement of the energy deposited by the ionizing particle in the LAr 
volume is obtained by collecting information from the last of the three wire planes, working in 
charge collection mode. 



Chapter 3 

The 50 liter Liquid Argon TPC 



The 50 liter Liquid Argon (LAr) Time Projection Chamber (TPC) is a detector buih and success- 
fully operated at CERN for R&D purposes within the ICARUS programme. In the year 1997 it 
was exposed to the CERN neutrino beam for the entire SPS neutrino run period as proposed and 
approved at the SPSLC of January 1997 [70]. The detector, complemented with scintillators acting 
as veto, trigger counters and pre-shower counters, was installed in front of the NOMAD detector. 
The collected data brought important information for a better understanding of the performance 
of Liquid Argon TPC's which should be useful for the entire ICARUS programme. 



3.1 The geometrical layout 

The detector structure consists in a stainless steel cylindrical main vessel, 70 cm diameter, 90 cm 
height, whose upper face is an UHV flange housing the feed-through's for vacuum, liquid Ar filling, 
high voltages and read-out electronics (Fig. 13. ip . Inside the main vessel, an ICARUS type Liquid 
Argon TPC is mounted. The TPC has the shape of a parallelepiped whose opposite horizontal 
faces (32.5x32.5 cm^) act as cathode and anode, while the side faces, 47 cm long, support the field- 
shaping electrodes (Fig. l3.2| ). The mass of the Liquid Argon contained in the active volume is 69 kg 
(T = 87 K at 1 atm, p = 1.395 g/cvc?). Ionization electrons produced by the passage of charged 
particles drift vertically toward the anode by means of a constant electric field of 214 V/cm. 

The read-out electrodes, forming the anode, are two parallel wire planes spaced by 4 mm. Each 
plane is made up of 128 stainless steel wires, 100 /im diameter and 2.54 mm pitch. The first plane 
(facing the drift volume) works in induction mode while the second collects the drifting electrons 
(see Sec. I2.1.ip . The wire direction on the induction plane runs orthogonally to that on the 
collection plane. The wire geometry is the simplest version of the ICARUS readout technique [SS] 
since both the screening grid and field wires in between sense wires have been eliminated. The 
wires are soldered on a vetronite frame which supports also the high voltage distribution and the 
de-coupling capacitors. 

The cathode and the lateral field-shaping electrodes are copper strips (5 mm thick and 1.27 cm 
wide) positioned on a vetronite support with printed board techniques. The support was glued 
on a honeycomb structure to ensure rigidity. The distance between two adjacent strips is 10 mm. 
The 10 kV drift voltage is distributed to the strips by means of 100 Mil resistors. In Tab. 13. ll we 
have summarized the main technical data of the 50L detector. 
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Figure 3.1: The 50 liter Liquid Argon Time Projection Chamber. 



3.2 The purification system 

A crucial working parameter for a TPC operating with liquefied rare gases is the lifetime of free 
electrons in the medium. Primary electrons produced by the passage of charged particles drift 
toward the anode crossing macroscopic distances. In the present case, for an ionization electron 
created in the proximity of the cathode, the drift path-length exceeds 46 cm. The effectiveness 
of charge collection at the readout plane is related to the purity of the Argon since electron-ion 
recombination is mainly due to oxygen molecules present in the LAr bulk [67j . The contamination 
of electronegative molecules must be at the level of 0.1 ppb to allow drifts over 0(1) meters. 

This has been achieved using commercial gas purification systems to remove oxygen and polar 
molecules (H2O, CO2, fluorinated and chlorinated compounds) combined with ultra-high-vacuum 
techniques to avoid re-contamination of the liquid through leaks and with the use of low out-gassing 
materials for the detector components. Electron lifetimes higher than 3 ms (0.1 ppb) are easily and 
constantly reachable. Liquid phase purification allows fast filling of large detectors. Continuous 
recirculation of the liquid through the purifier during normal operation is used to keep the lifetime 
stable (against micro- leaks and out-gassing) [H31 E3 EI] . 

As mentioned in previous Sec. 13.11 the active part of the detector is located inside a stainless 
steel cylinder. The connection to the outside area is obtained through a set of UHV flanges housing 
the signal, the high voltage cables and the vacuum feed-through. The cylinder is positioned into 
aim diameter dewar partially filled with low purity Liquid Argon acting as a thermal bath. The 
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Figure 3.2: Sketches of the 50 hter chamber. Side (left) and top (right) views.. 



Inner dimensions 


32.5 cm X 32.5 cm x 47.0 cm (length, width, height) 


Inner volume 


49645 cm^ w 50 1 


Active mass 


67 Kg 


Sensing planes 


Induction & Collection (mutually orthogonal) 


Wire number 


256 (128 per sensing plane) 


Wire spacing (pitch) 


2.54 cm 


Wire diameter 


100 fim 


Max. drift length 


47 cm 


Electric field 


214 V/cm 


Signal sampling rate 


@ 2.5 MHz during 820 ^is (=2048 time samples) 



Table 3.1: Some technical data of the SOL detector. 



ceiling of the external dewar is in direct contact with air. The level of Argon in the dewar is so 
small that the pure Ar in the inner cylinder can evaporate as well. 

The detector is equipped with a standard ICARUS recirculation-purification system [64]. In 
this experiment the Argon purification is carried out in the gaseous Argon phase through the 
recirculation system, while during the initial filling up stage is done in the Liquid Argon phase. 
Once the filling up of the chamber finishes, the recirculation system takes care of increasing the 
purity level first and then to keep it stable during the data taking period. This purification 
procedure in the recirculation phase is outlined in a simple way as follows (see Fig. 13. 3|) : 

• The gaseous Ar in the active volume crosses the feed-through and reaches an Oxisorb filter 
for purification, thanks to pressure gradients. 

• Now, the pure Argon goes through a cooling coil (inside a second buffer cooled by low-purity 
LAr) where it is liquefied again. 

• Finally, the pure LAr is collected back into the active volume. 
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Figure 3.3: Scheme of the 50 hter Liquid Argon purification system. 



The Argon purity is characterized by the hfetimc of the free drifting electrons (r), defined by: 

Q{t) = Qoe-i (3.1) 

Qo being the primary deposited ionization charge and t the electron drift time. The mean free 
path A of these drifting electrons is nothing but the lifetime times the mean drift velocity, which 
have been measured in this experiment (see Sec. 15.2. ip to be 0.91 mm//zs. 



The purity monitor 

The effectiveness of the purification system through direct measurements of r has been monitored 
during data taking by a dedicated setup first developed by the ICARUS Collaboration in 1989 [7^ . 
It consists of a small double-gridded drift chamber located below the readout planes (see Fig. l3.4p . 
Electrons are generated near the cathode via photoelectric effect driven by a 20 ns UV laser pulse. 
Each pulse produces a bunch of lO"^ electrons. They drift toward the anode along the electric field 
lines and cross a 50 cm drift region between the two transparent grids. The ratio of the induced 
signal in the proximity of the grids provides a real time estimate of r during data taking by means 
of the following relation: 

Qa ^ Tc sinhj^) 
Qc Ta sinh{^) 



Td + 



Ja-t--tc 



(3.2) 



3.2. The purification system 



35 



where Td is the time the drifting electrons spend from the first grid to the second, Tc and Ta are 
the time from the cathode to the first grid and to the second grid to the anode respectively. 

The relation (|3.2p can be understood in a very simple way: When the electrons drift between 
cathode and window, the amplifier receives two equal and opposite currents that cancel each other. 
When the electrons drift between the window and the cathode grid a positive current flows in the 
amplifier; as a consequence its charge output increases linearly reaching a value Qc, the charge 
leaving the window. When the electrons drift between the grids, no current flows in the amplifier 
(at least if the grids screen perfectly). Finally when the electrons drift between the anode grid 
and the anode, the current is negative. The output charge signal has a negative step, whose 
height is the charge Qa reaching the anode. If no charge is lost in the drift volume (zero impurity 
concentration), obviously Qa = Qc', if some charge is lost, corresponding to a lifetime r the ratio of 
the charges is given by Eq. (|3.2|) . In the usual case that r >> Tc, T^, we have a simpler expression: 



^' + ^^ (3.3) 



J-, Qa Tc 

An example of the charge output is shown in Fig. 13.41 (right). From a direct measurement of the 
purity monitor output signal (Fig. 13. 4p . the collected charge Qa, Qc and the time intervals T^, 
Tc and Ta can be obtained and thus, the electron lifetime t. The initial lifetime during filling 
was about 100 fj,s. The recirculation/purification system, circulating about 5 liters of LAr/hour 
allowed to increase this value to more than 8 ms in three weeks (see Fig. 13. 5p . The total Argon 
consumption necessary to circulate the pure liquefied gas and to compensate for the heat losses 
was ~ 200 liters per day. 




Figure 3.4: (Left) The purity monitor. (Right) The output charge signal. 



An independent estimate of r can be obtained monitoring the charge attenuation for muons 
crossing the chamber at various heights (sec Sec l5.2.2]) . Both methods provided consistent results 
and the electron lifetime averaged along the whole data taking period (9 months) turned out to be 
always higher than 10 ms. This performance was excellent since the maximum drift time was about 
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Figure 3.5: Lifetime of the drifting electrons in the 50 Uter Liquid Argon TPC. The fihing of the 
chamber with LAr was performed on the fourth of April 1997. 



500 jis and therefore the attenuation of the ionization over the drift distance can be considered 
negligible. 



3.3 The photosensitive dopant 

As we discussed before, an ionizing particle in Liquid Argon produces, along the track, electron-ion 
pairs and excited atoms {Ar*). The excited atoms decay, in few hundred nanoseconds, emitting 
ultraviolet photons with energy distribution peaked at 128 nm (9.7 eV). In the same time scale and 
depending on the ionization density and electric field, some of the pairs recombine also emitting 
UV photons. This produces saturation of the free charge as a function of the deposited energy and 
electric field. A possible solution to avoid this consists in dissolving in Liquid Argon a photosensitive 
dopant able to convert back the scintillation light into free electrons, thus linearizing the dependence 
of the charge on the energy and the electric field. 

During this test, the Argon has been doped with Tetramethyl-Germanium as photosensitive 
dopant because of the following advantages [73] : 

• TMG is not absorbed by the Oxisorb filter in the recirculation system. 

• TMG can be easily purified to an electron lifetime level better than 10 /is. 

• TMG has a large photo-absorption cross section (62 Mbarn) and large quantum efficiency 
(close to 100 %); this implies that small quantities (sa ppm) of TMG are enough to convert 
all the de-excitation photons into electrons in the vicinity of the ionizing track. 
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Figure 3.6: Correlation between collected charge and deposited energy at an electric field of 200 
V/cm for pure Argon and with a TMG concentration of 3.5 ppm. 



• TMG, at the ppm level, dissolves homogeneously in Liquid Argon. 

The relation between the collected charge and the energy deposition was obtained in the R&D 
phase of the 3-ton ICARUS prototype (see [73j) from the analysis of the stopping muon and 
proton events at several electric fields and TMG concentration of 1.3 ppm and 3.5 ppm. Fig. 13.61 
illustrates such relation before and after the doping. 

Argon excited atoms {Ar*) produced by ionizing radiation lead to low excited dimer 
formation through collision with Ar atoms. The rise-time corresponding to excimer formation and 
relaxation are in the sub-nanosecond range. For relatively low concentrations of the photodopant 
TMG the process which drives the recovering of drifting electrons from the scintillation light is a 
photoionization which follows the photocmission: 

Ar2 Ar + Ar + j 

7 + F F+ + e- 

It is evident that the linearity of the detector response, after doping with TMG, is significantly 
improved especially at low fields where the recombination effect in pure Liquid Argon is higher. 
Continuous recirculation of the Argon through the purifier keeps the level of contamination stable 
against micro-leaks or out-gassing. Moreover no degradation of the space resolution has been 
measured, implying that the photo-conversion happens very close to the track. 

3.4 Read-out set-up and data acquisition 

The mean charge per unit pitch (2.54 mm) for a mip crossing the chamber horizontally corresponds 
to about 2.3 x 10^ electrons (^^4 fC) assuming no charge loss along the drift volume. The smallness 
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of the signal and the absence of an ampHfication phase requires very low noise charge pre-amplifiers. 
This has been reached by means of JFET technology. 

The pre-amplifiers 

There is one preamplifier per read-out channel (256 wires), directly placed in the vctronite frame 
which supports also the read-out wires. Since the noise depends on temperature and input ca- 
pacitance, a further improvement is obtained placing the pre-amplifiers into the liquefied gas (at 
-186°C). The polarization voltage (15 V) is distributed separately to each preamplifier to minimize 
the number of dead channels in case of failure of one of the components during data taking. 

Each one of these pre-amplifiers is compound by an integrator with a time constant of ti ~ 
200 fxs and a CR output impedance with time constant of r2 = 1 ms: The signal of the input 
current is integrated through the minimization of the input impedance in conjunction with the 
appropriate frequency response of the preamplifier (this minimization also allows a reduction of 
possible interferences between the signals of the different channels). 

The amplifiers 

The amplifiers gather the signal from the pre-amplifiers with typical gain of 18 mV/fC and an 
output impedance of 120 H.. They are of TL.026 type, with a differential stage at the input and the 
output followed of a buffer stage. There were located outside the chamber at room temperature. 

During the whole data taking period the amplifiers have been operating in current (or "quasi- 
current" ) mode, which allows a rather good visualization of the ionizing events, with a signal to 
noise ratio of 11 for mip signals (see Fig. 13. 7p and no saturation of the electronic response even in 
the occurrence of highly ionizing events (i.e. electromagnetic showers). 

Data acquisition 

The amplified signal is brought up to a set of fast 8-bit ADC (Analog to Digital Converter). The 
signal is sampled with a 2.5 MHz frequency for a duration of 820 /is (2048 time samples), while 
the highest possible drift time (primary ionization at the cathode) amounts up to 500 fis (1250 
time samples). Once a trigger is given, the data stream is recorded into a buffer. It consists on 256 
read-out channels with 2048 time samples; the to signal is also recorded as an squared pulse in a 
virtual additional channel. The arrival of a subsequent trigger causes the switch of the data stream 
to another buffer (up to 8 buffers are available). The multi-buffer writing procedure minimize the 
dead time for signal recording to less than one sampling time (i.e. 500 fis). The events are written 
in raw-data format without any zero suppression; they are stored locally on disk and automatically 
transferred to the main CERN tape facility (CASTOR) using the network. 

The signal shape 

With this read-out configuration the visibility of the events appeared to be optimal in both collec- 
tion and induction views. In Fig. 13.71 the raw data output visualization is shown for a multi-prong 
neutrino event. The optimization of the analog electronic setup leads to rather clear and sharp 
raw images before any additional filtering algorithms were applied. Besides, the induction signal 
has been enhanced varying the ratio of the fields in the drift and gap regions. 

The measured signals are induced by the moving charge of the drift electrons and are propor- 
tional to the drift velocity. Thus, strictly speaking the signal from the induction plane is bipolar: a 
small negative part from the slowly drifting electrons approaching the plane from the drift volume. 
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and a much larger positive part induced when the electrons have passed through the plane and 
drift away much faster in the larger drift field among the induction and collection planes. 

The careful optimization of the instrumental setup leads to a satisfactory level of performance. 
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Figure 3.7: Output visualization of a multi-prong neutrino event in "quasi-current" mode. [Top] 
The raw data from collection and induction wire planes. Below two examples of a single channel 
(wire) electronic read-out in collection [middle] and induction [bottom] modes. 



Chapter 4 



Off-line event reconstruction 



The goal of the reconstruction procedure is to extract the physical information contained in the 
wire output signals, i.e. the energy deposited by the different particles and the point where such a 
deposition has occurred, to build a complete three dimensional spatial and calorimetric picture of 
the event. The basic building block of a track is called hit, defined as the segment of track whose 
energy is read by a given wire of the readout wire planes. Therefore, the spatial and calorimetric 
information of the track segment is contained in the associated hit, and the sensitivity of the 
detector depends entirely on the hit spatial and calorimetric resolutions. 

In the present chapter we describe a method developed to perform the hit spatial and calori- 
metric reconstruction. It is based on the pioneer work done for the stopping muon sample collected 
with the T600 detector [74j , however many new contributions and improvements have been devel- 
oped during this work in order to extend the reconstruction capabilities to low multiplicity neutrino 
events collected with the SOL LAr chamber (see Chapter [SJ. 

The reconstruction algorithm proceeds through the following steps: 

1. Hit identification (Sec. 14. 1( ): The hits arc independently searched for in every wire as output 
signal regions of a certain width above the baseline output value. 

2. Fine hit reconstruction (Sec. 14. 2| ): The parameters defining the hit (position, height, area), 
which contain the physical information, are precisely determined. 

3. Cluster reconstruction fSec. 14.3"]) : Hits are grouped into common charge deposition distribu- 
tions based on their position in the wire/drift coordinate plane. 

4. Two-dimensional (2D) track reconstruction (Sec. 14. 4^ : The 2D projections of the ionizing 
tracks are detected taking care of interaction vertexes and splitting the clusters into chains 
of hits with smooth transitions between them. 

5. Three-dimensional (3D) reconstruction (Sec. 14. Sp : The hit spatial coordinates are recon- 
structed using the association of hits from different views to common 3D track segments. 

6. Calorimetric reconstruction (Sec. l4.6() : The spatial reconstruction of the tracks together with 
the information extracted from the fitting of the hits allows to precisely measure the energy 
deposited by a particle along its path. 

In the next sections, we describe in detail the different steps of the hit reconstruction procedure. 
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4.1 Hit identification 

4.1.1 Hit identification algorithm 

The hit identification aims at distinguishing signals produced by ionization electrons from elec- 
tronic noise. Hits are identified as signal regions of a certain width with output values above the 
local baseline (defined as the average output value in a given signal region). The hit identification 
algorithm typically acts over the whole wire sample of a given readout plane. In the SOL config- 
uration there are 2 wire planes with 128 wires. Each wire is sampled 2048 times per event. The 
hit identification is based on basic geometrical principles of the wire signal and establishes rough 
preliminary bounds to the hit range along the drift coordinate, which will be of utmost importance 
in ulterior steps of the reconstruction procedure. No information from adjacent wires is used at 
this stage. 

The algorithm loops over the selected group of wires (i.e. a whole wire plane), skipping those 
with identified problems, such as large noise conditions or disconnections, where no signal identi- 
fication is possible. Hits are searched for in the wire output signal after a low frequency filter has 
been applied. The different steps are described in the following paragraphs. 

Low &: High frequency noise filter 

Before the use of the hit identification algorithm is better to have the best signal we can for this 
purpose. In order to remove the noise contribution from the wire signals it is worth to perform a 
filter in frequencies. Such a filter is based on the Fourier transform of the raw wire signal. The 
Fourier spectrum reveals that the main contribution of drifting electrons coming from a ionizing 
particle to the wire signal relics on quite low frequencies. The other contributions to the read-out 
signal arc due to electronic noise, which can be partially removed performing a suitable cut in 
frequencies. The Fourier transform of the wire signal in the 2048 output samples is computed 
using the Fast Fourier Transform (FFT) algorithm 

The amplitude of the low frequencies is reduced by convoluting the result with the following 
function of the frequency v ("soft" step function): 

S{f) = 1 - 4« (4.1) 

e " 

where R and a are the parameters determining, respectively, the mean frequency radius and the 
transition region thickness. Setting the parameters to R = 4.5 kHz and a = 0.9 kHz effectively 
cuts low frequencies up to about 6.5 kHz while reducing the diffraction effect. This effect distorts 
the whole output range when a hit of high amplitude and sharp edges is present in the wire. The 
FFT filter efficiently removes the low-frcqucncics modulation of the signal giving as a result more 
uniform and constant noise baselines. 

In the same way, the FFT filter is used to deal with the annoying high-frequency noise removing 
the contribution of waves with a frequency higher than 70 KHz. This make the search of hits much 
easier, above all in the case of small signals which are little bigger than the electronic noise. 

Fig. 14.11 (top) shows the raw data before the application of the filter. Fig. 14.11 (bottom) shows 
the result of the application of the FFT filter with the cuts pointed above. 

Hit search 

The hit search proceeds on the filtered wire output signal. The identification is based on the 
geometrical features of the wire output in a given region, i.e. the presence of a relatively broad 
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Figure 4.1: Effect of the frequency filter on the wire signal. (Top) Raw data with no filter applied. 
(Bottom) The filtered wire. The contribution of plane waves with frequencies below 6.5 KHz and 
above 70 KHz have been removed. 



region of output values above or below the baseline. The algorithm detects the pattern of a hit 
in both Collection and Induction planes and determines the type of the signal shape. This is 
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very important for the Induction plane where different patterns of signals (unipolar or bipolar) 
appears depending on the energy deposition. The behavior of Collection signals is always unipolar. 
Correlations with the signal from neighboring wires are not considered. 
For every wire, the hit search algorithm follows the following steps: 

• Every output sample is compared with the baseline. A value above (below) the baseline by 
more (less) than 3 — 4 ADC counts triggers an unipolar (bipolar) hit candidate. 

• The hit candidate is built with all the subsequent output samples after a threshold is given. 
Depending on the detected type of hit, the algorithm search for an unipolar or bipolar pattern 
in the wire signal. 

• The hit candidate is characterized by its width, i.e. the distance (in drift samples) between 
the hit initial and final points. Rejection of fake (noise) candidates is achieved by imposing 
a minimal width value. 



• Further rejection can be carried out depending on the expected shape of the hit. For Col- 
lection wires, where the hits are expected to have an exponential falling slope, an extra 
requirement on the minimum distance from the peak position to the hit end is imposed. 
For Induction wires, when a bipolar hit is detected we require in addition that the distance 
between the down and the up peak is below a minimum. 

Fig. 14.21 illustrates the two different types of hits which we are dealing with and the parameters 
used to identify them. Fig. 14.31 illustrates the performance of the hit searching algorithm in both 
wire planes. It is the result of an optimized choice of the parameters for the 50L readout. As it 
was desired, the hit search procedure is able to detect the most of the real hits, faking as less hits 
as possible. 
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Figure 4.2: Example of a hit produced by an ionizing track on a Collection wire (left) and on an 
Induction wire (right). Marked are the parameters used in the hit search. 



Resolving close hits 



Close hits may appear when segments of tracks belonging to different ionizing particles produce a 
signal in the same region of the wire/drift plane (see Fig. 14.41 (left)), or when a single ionizing par- 
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Figure 4.3: Example of performance of the hit finding algorithm. (Top): gray scale mapping of 
the output signal in the Collection view (left) and the Induction view (right). (Bottom): in red, 
the hits found with an optimized choice of the hit searching parameters. 



tide travels with small angle with respect to the drift direction (see Fig. 14.41 fright)). Identification 
of close hits results in a better hit calorimetric and spatial reconstruction of those tracks. 

The resolution of close hits is carried out within the hit finding procedure. When a hit candidate 
is found, an embedded algorithm looks for patterns which point out the presence of another close 
hit within the hit being analyze. 

In the case of Collection signals, when a falling and a rising slope are found consecutively within 
a hit candidate, the algorithm splits the candidate into two hits separated by the starting point 
of the rising slope. Although for bipolar signals in Induction wires this criteria is a little bit more 
complicated, the search for close hits is based on the same strategy: the presence of a close hit is 
found when an unexpected change of slope is find out. Once the close hits are found, the selection 
criteria arc applied separately to both candidates. 

Hit parameterization after the hit finding algorithm 

After the search for hits in both wire planes they are stored together with their main geometrical 
features for further processing at the next stage of the reconstruction procedure. 

The hit identification provides the following hit parameters for unipolar signals (sec Fig. [53] (left)): 

• The view and wire indexes of the hit. 
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Figure 4.4: (Left) Example of close hits occurrences caused by close tracks in the same region. 
(Right) Close hits from a single track in the direction of the drift. 




Figure 4.5: Parameters characterizing a hit after the hit finding algorithm for unipolar (left) and 
bipolar signals (right). 

• The drift coordinate of the peak. 

• The drift coordinate range of the hit, defined by the initial and final drift samples computed 
in the hit expansion. 

• The height of the peak taken as the difference (in ADC counts) between the peak and the 
baseline. 

In the case of bipolar signals, in addition to the these parameters, the following ones are also 
kept (see Fig. |13] (right)): 

• The drift coordinate of the down peak. 

• The height of the down peak. 



4-2. Fine hit reconstruction 



47 



• The drift coordinate where the signal cross the baschne. 



4.2 Fine hit reconstruction 



Once the hit has been detected, a fine hit reconstruction is performed aiming at extracting in an 
optimal way the parameters defining the hit (position, height, area) which contain the physical 
information of the original associated track segment. The hit spatial reconstruction is entirely based 
on the determined hit peak position, whereas the hit area in Collection wires is proportional to 
the energy deposited by the ionizing particle, and therefore constitutes the base of the calorimetric 
reconstruction. In the present section we describe the fine hit reconstruction algorithm applied to 
the typical unipolar signals from Collection plane. No fine hit reconstruction is performed over 
bipolar Induction hits, and just the geometrical information from the hit identification phase is 
enough for the next stages of the reconstruction method. 



4.2.1 Hit fit reconstruction method 

The adopted approach consists of performing a fit of the output signal (the raw one) around the hit 
region with an analytical function that reproduces well the hit shape, and extracting the parameter 
values from the fitted function. The aim is to extract from them, in the most accurate possible 
way, the hit peak position, height and area. 

The output signal within the hit window is fitted using the following analytical function of the 
drift time t: 

m^B+A ' (4.2) 

1 + e ^2 

where B is the baseline, A the amplitude, to the point for which the height of the function with re- 
spect to the baseline is equal to A/2, and ti and T2 are related to the falling and rising characteristic 
times, respectively (see Fig. 14. 6p . 

For a given hit, the values of the fit parameters are those obtained by minimizing the of the 
fit of f{t) to the output signal within the hit window. The minimization is performed using the 
MINUIT package [T^. 

The hit area is obtained by numerical integration of the fit function in the hit window, while 
the peak position (tmax) and height {/{tmax)) are obtained analytically through the following 
analytical expressions: 



to + At 

e ^1 



(4.3) 

f{tmax) — B + A TaT 

1 + e ^2 



with 

n - T2 



At = T2 In ■ 

Fig. 14.71 shows some examples of fitted hits corresponding to different events occurring inside the 
LAr active volume. 



Hit parameterization after the hit fine reconstruction algorithm 

The hit fine reconstruction adds the following quantities to the hits' parameterization: 
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Figure 4.6: Function used in the hit fit (see Eq. (|4.2p ) for the following parameter values: 
B = 5 ADC cuts, A = 20 ADC cnts, to = 14 ^J-s, n = 3 ^is, T2 = 1 ^is. 

• The local mean value around the hit. 

• The bounds of the hit window, where the different reconstruction procedures are carried out. 

• The peak position, the peak height and the hit area as computed by the direct ADC counting 
method. 

• The peak position, the peak height and the hit area as computed by the hit fit method. 

• The values of the fit parameters values. 

• The x^/n-d.f. of the fit. 



4.3 Cluster reconstruction 

A cluster is defined as a group of adjacent hits within the wire/drift coordinate plane. The goal 
of the cluster reconstruction is to perform a first grouping of hits belonging to common charge 
deposition distributions, such as tracks or showers. Clusters provide identification criteria for the 
different patterns, and thus determine which reconstruction procedure must be followed. Clusters 
also provide a criterion for discrimination between signal and noise hits as well as a method to 
recover undetected hits through an exhaustive search. 
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Figure 4.7: Fitted signal from a hit given by a proton near its stopping point (top-left). The 
fitted hit of another stopping proton together with two hits from mips (top-right). An example of 
detected and fitted hits of a parallel shower of mips (bottom). 



4.3.1 Cluster reconstruction algorithm 
Linking the hits 

Clusters are groupings of hits which are close to each other, and two hits are close when exists a 
link between them. A link is defined as a pair of hits coming from consecutive wires and having 
overlapping drift coordinate ranges, where the drift coordinate range of a hit is defined as the 
one computed by the hit finding algorithm (see section [4. ip . Therefore, clusters can be defined as 
groups of hits with links between them. The links allow to define the topological features of a hit 
or a hit grouping, which will be of utmost importance during the 2D and 3D reconstruction stages 
(Sec. 1131 and Sec. SI]). 
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To this purpose, a search for all the hits close to each other is performed, obtaining for any 
hit the set of those hits which are surrounding it. In such a way the links arc built as logical 
connections between a hit and its close hits. An example of how the links are established is shown 
in Fig. 14.81 In terms of these links, the bounds of a cluster are defined as those hits having one 
link at most. 




Figure 4.8: The pictures show the signal of a neutrino interaction in a region of the Collection 
plane. The drift coordinate is in vertical and the wires in horizontal. (Left) The raw data of a 
vertex from a neutrino interaction. (Right) The found hits in red; the formed links between them 
are in yellow. Because all the hits are connected by links to each other, they form a cluster. 



Building the clusters 

Clusters are built in a two step approach. In the first step (pass I), preliminary clusters are built as 
groups of hits being connected by links. The second step (pass II), attempts to expand and merge 
together the clusters found in pass I, in order to reduce the number of the resulting clusters. The 
reconstruction algorithm takes as input a list of hits from a given wire plane together with their 
links, and constructs all possible clusters out of them, using only the information obtained by the 
hit finding algorithm (see section H?T|) . 

Pass I It is essentially established above: A preliminary grouping of the hits into clusters is per- 
formed: Hits are grouped into a cluster if any of them is connected by a link with at least 
another one belonging to the cluster. 

Pass II The goal of cluster reconstruction pass II is to expand and join cluster fragments produced 
by pass I. It proceeds as follows: the regions around the cluster borders are examined (where 
a cluster border is defined as a hit linked at most with one hit within the cluster) , searching 
for new hits which could be added to the cluster, hence reducing the cluster fragmentation 
caused by the misidcntification of hits during the hit search phase. Attempts for clusters 
merging/expansion are performed both along the wire and the drift time coordinate direc- 
tions. 

The algorithm for the cluster expansion proceeds in the following sequential steps: 
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• Cluster borders are identified as those hits with at most one link to the cluster. 

• The direction of the search is determined as wire coordinate increasing (decreasing) for bor- 
ders linked with a hit from the previous (next) wire. For single hit clusters, both directions 
are consecutively inspected. 

• The slope of the search is determined as the difference between the drift time coordinates of 
the cluster border hit and its linked hit. For single hit clusters, this slope is equal to zero. 

• New hits are searched for in the consecutive wire along the determined direction using the 
same hit finding algorithm, but with looser parameters. Those new hits which were close to 
a previously existing one in the cluster are automatically joined to that cluster, creating the 
corresponding new link. In this way a cluster is expanded finding new misidentified hits in 
a clever and effective way: this exhaustive search (with loose parameters) is performed just 
in the region where is expected to possibly have some new hits, but not in the overall data 
sample like in Sec. 14.11 

• The latter expansion can be continued until no more hits are found or when it reaches to 
another existing cluster. In this case both clusters are merged. 

• The borders of the new resulting cluster are identified and the whole process is repeated until 
no further expansion is possible. 

In Fig. l4.9l the reconstruction of cluster is exemplified. In Pass I, the algorithm establishes links 
between the hits and construct preliminary clusters; In Pass II, the algorithm starts the exhaustive 
search for hits in the directions pointed by the yellow arrows (Fig. 14.91 (top)). The final result is 
shown in (Fig. 14.91 (bottom)), where the bigger final cluster is marked in yellow (other clusters are 
not shown). It is worth to underline how misidentified hits are recovered in this second pass. 

Hit parameterization after cluster reconstruction 

Cluster reconstruction adds the following quantities to the hit parameters : 

• The cluster identification number. 

• The number of links associated to the hit. 

• The list of hits to which the hit is linked. 

These data are of utmost importance during the next step of the reconstruction procedure (see 
Sec. EH and Sec. EH). 

4.4 2D track reconstruction 

The 2D track reconstruction algorithms work separately in each wire plane in order to identify the 
2D projections of the real ionizing tracks. The goal is to detect smooth chains of hits which can 
be identify as track segments. To this purpose two different approaches have been implemented: 
The Tree Algorithm [77] and the Neural Network based approach [78l [79l [80j . A proper recognition 
of 2D track projections in both (Collection and Induction) views will lead to vertex identification 
and to an efficient 3D track reconstruction procedure (Sec. 14.5^ . 

The description of the 2D track recognition methods and how they perform over the 50L TPC 
data is presented in the following sections. 
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Figure 4.9: The pictures show the signal of a neutrino interaction in a region of the Collection 
plane. The drift coordinate is in the vertical axis and the wires in the horizontal axis. (Top): The 
raw data of some prong particles from a neutrino (DIS) interaction. We represent in red the hits 
found before the cluster reconstruction procedure. Yellow circles and arrows show the border hits 
and the direction of the exhaustive hit search. (Bottom): The found hits (in red) after the cluster 
reconstruction method. One of the built cluster is singled out with its hits in yellow. 

4.4.1 Chains of links 

Both 2D finding methods operate into each one of the found clusters (see Sec. 14. 3|) . retrieving the 
information of the hits and the links between them. Since a track segment can be thought of as an 
unbroken chain of consecutive hits, all possible links between a hit and its adjacent are collected 
into chaint0. In principle there are many possible chains of links inside a cluster and only some of 
them properly define a 2D track projection. The goal of the track finding methods is to efficiently 
selects the real track projections from the whole sample of chains. An example of some of the 
chains which could be formed from the links in a cluster is shown in Fig. 14.101 The two chains at 



chain is defined as a Collection of Units consecutively connected 
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the bottom would represent real 2D track projections. 




Figure 4.10: Examples of some different chains (in green) which could be constructed from the 
links established (in yellow) in a cluster (see Fig. l4.8p . The four chains shown are track candidates, 
but only the two ones at the bottom would be real 2D track segments. 



4.4.2 The Tree Algorithm 

The Tree algorithm was successfully applied to the pattern recognition problem in multi-layer 
track chambers (see [77]). Because of the hits from a ionizing particle in a LAr TPC also have a 
layered structure (the wires readout), the approach to the track finding problem developed in [77] 
is directly extensible to our purposes. 

The Tree algorithm works on the hits of a certain cluster using the information provided by 
the links and constructs from them all the possible chains (candidates to track segments) in the 
following way: 

• The algorithm starts from a link which contains a border hit (root link) and determines its 
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elementary tree. The elementary tree of any given link is composed of this link, which is then 
called the trunk of the elementary tree (see Fig. 14.111 (top)), and all links having the middle 
hit in common and with approximately the same values of track parameters (in this case we 
have only one track parameter: the slope of the link). These links are then called branches 
of the elementary tree. Besides, all the possible elementary trees that can be consecutively 
connected with the starting one are also determined (see Fig. 14.111 (right)). 

• The collected elementary trees are connected forming the full tree associated to the root link. 
The links which define the limits of a full tree are called leaves (see Fig. 14.111 (bottom)). The 
algorithm recursively climbs the full tree from the root link, locating all chains from the root 
to any leaf but never finding the same chain twice. In this way tracks (or more precisely, 
reliable track candidates) can be recognized. Because of the branches of any elementary tree 
must be compatible with the slope of the trunk, any chain built in this way also have smooth 
transitions between their links. 

• The real track segment is selected from the chains belonging to the full tree by means of a 
selection criteria. In this case, a combination of two requirements is imposed: The longest 
and smoothest chain is selected as the track segment. 

• After the selection of the track, their links are stored, tagged and they do not participate 
any more in the algorithm (under the assumption that a link is never shared by two different 
track segments). Then, the algorithm continues processing in the same way the remaining 
links, starting from another root and detecting new track segments. 

Examples of elementary trees are shown in Fig. 14.111 (right) The elementary trees can be 
combined to form a full tree as shown in Fig. 14.1 If bottom). The correct track is immediately 
recognized as the smoothest chain of links in the full tree (chain in red in Fig. 14.111 (bottom)). 



Links: 



Elementary Trees: 
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A Full Tree: 




Figure 4.11: (Top) The starting point: a Collection of links in a cluster. (Right) All the possible 
elementary trees. (Bottom) The full tree is built chaining the elementary trees. 
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4.4.3 The Neural Network approach 

The Neural Network based computing techniques are powerful tools for pattern recognition prob- 
Icmtl^, which have been successfully applied to many problems in experimental particle physic^. 
The track finding problem is a task of this type: given a set of hits (signals) in space, the aim is to 
reconstruct particle trajectories subject to smoothness constraints. The Neural Network approach 
for the track finding in the LAr TPC here presented is based on the work developed in [78] and 
|80j . where an algorithm based on the Hopficld-style recurrent neural network [55] is discussed. 

A neural network consists of processing units, which are called neurons, and connections be- 
tween them. The state (activation) of a neuron is a function of the input it receives from other 
neurons. The output of a neuron is usually equal to its activation. Therefore, the state of each 
neuron depends on the state of the neurons from which it receives input and on the strength 
(weight) of these connections. Neural nets can learn representations by changing the weights of 
the connections between the neurons. In the following, the method specifically developed for the 
LAr TPC is described. 



The neural network algorithm 

The links between the hits in a cluster play the role of the neurons, which can be in two possible 
states Sij = 1 if the link i j is part of the track segment and Sij = if this is not the case. 
The connection between the neurons (links) is established through a term called synaptic strength 
Tijki, which involves the (ij) link with the (kl) one. The dynamics of such a Hopficld type neural 
system is governed by means of the minimization of the following energy function: 



^^'^ijkiSijSki (4.4) 



E -- 

2 

ijkl 

The following step is to find a suitable form of the synaptic strengths Tijki which performs the 
required task. 

Because of a track with N hits can be considered as a set of iV — 1 links with a smooth shape and 
without bifurcation, the synaptic strenghts Tijki have been chosen to take care of this requirement: 

rp rpicost) rpiconstraint) i. 

J-i]kl — J-ijkl ^ijkl \^-^) 

where Tj^'^^f*'' is such that short adjacent links with small relative angles are favored: 

Tl-k^'^^S.k^ (4.6) 

where m is an odd exponent, 9iji is the relative angle between the link (ij) and the (jl), rij is the 
length of the {ij) link and n an integer (see Fig. I4.12p . With this choice the energy of (|4.4p will 
blow up for long links and large angles. 

The constraint term consists of two parts: 

rp{constraint) rp(l) . rp{2) i. —\ 

^ijkl ~ ^ ijkl ijkl {^■') 

The first term in (|4.7p takes care of the requirement that there should be no bifurcating tracks: 
two segments with one identical index should not be "on" at the same time. This is done with the 

^See I81| for a general introduction on neural networks foundations 
^See |79j for an specific discussion on particle physics problems 
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Figure 4.12: Definition of link lengths rij and angles 9iji between links. 



following choice for T^jli : 



(1) 



ijkl 



(4.8) 



(2) 1 
The other term, T^^!;, ensures that the number of neurons that are "on" is roughly equal to the 

number of hits A^. It has the generic form: 



(2) 



1 



ijkl ~ 2^ 



(4.9) 



This configuration of the synaptic strenghts obviously is symmetric {Tijki = Tkuj) and with van- 
ishing diagonal elements (T!yij = 0). 

The final expression for the energy function is: 



E 



StjSki + -i^ce'^SijSii + i^a'^SijSkj + -/3 Ski 



N 



1^3 



k^i 



(4.10) 



The a and /3 factors are Lagrange multipliers entering in the definition of the energy function. At 
a given time the state of the system is given by the vector S = (..., 5*^, ...), while the dynamics of 
the system is determined by the T matrix of synaptic connections. By construction, the solution 
of the problem of track finding is the state of the neurons S which minimizes the energy function 
(|4.10p . The local updating rule which brings the system to a local minimum through the gradient 
descent method of the energy function is [S^] : 



Si = signiy^T^jS. 



(4.11) 



However, in general is necessary to reach the global minimum of the system for a better perfor- 
mance. One way to get this is to expose the system to a noisy environment, in which the the state 
vector S obeys the Boltzmann distribution: 



P(S) 



4cxp-^(«)/^ 



where the partition function Z 



Z = ^ exp 



-£;(S)/T 



(4.12) 
(4.13) 



is the sum over all possible states and T is the temperature of the system. 
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For an efficient and fast perform of tfie system relaxation to the minimum, is usual to make 
use of the mean- field theory (MFT) approximation [83] and get for the thermal average of Sij : 



Vij = {Sij)T = tanh 



1 dE 
'TdV~ 



tanhf i ^TyfciVfc; 



(4.14) 



Note that the local update rule in (|4.1ip can be obtained from (|4.14p taking the T ^ limit. V^j 
can take any value between -1 and 1, and is more convenient to limit the range to the interval [0, 1] 
in other to get no contribution from deactivated neurons. To this purpose, the following linear 
transformation is taken 



1 + tanh 



1 dE 
TdV~. 



The final update rule for the neurons is obtained combining Eqs. ()4.15p and (14.10p : 



(4.15) 
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(4.16) 



Given an initial state of the neurons, the dynamics of the system is fully determined by (|4.16p . 



Algorithm operation 

The neural network described above can be applied to the 2D track recognition problem in a certain 
cluster. It takes as input the links between the hits (neurons) as they were introduced in Sec. 14.3.11 
(each hit is just linked with its close hits) and, after the iteration process described by the update 
rule in (|4.16p . the set of "active" links will define the track projections. The algorithm operates as 
follows: 

• All the neurons (links) are first initialized to an active state, Sij = 1 (|4.13l (top)). This 
configuration performs better than a random initialization in the particular case of the in- 
teractions we have in the LAr TPC. This is because of such state is closer (in general) to the 
optimum state, which minimizes (|4.10p . 

• The update of the neurons is performed asynchronously through Eq. (|4.16p . This means that 
a neuron is updated after the other: the current updating neuron "sees" the other neurons 
in their actual state. 

• The process is iterated until a certain level of convergence is reached. The adopted criteria for 
convergence is based on the changes of the activation state of all neurons after each iteration: 



• After convergence, the state of a neuron Sij is activated if its mean value Vij is bigger or 
equal than 0.5 and deactivated if lower (|4.13l (middle)). 

• Finally, the 2D tracks projection are identified as chains of consecutive activated links 
(Fig.EIIKbottom)). 

There are some free parameters in the definition of the update rule in (j4.16|l : The a and (3 
factors can be modified in order to give more or less relative importance to the different terms in 
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(|4.16p . The exponents in the T^^l^j^ term of (|4.6p rule the smoothness requirement for the track 
segments. These parameters have to be defined in such a way that short and almost straight links 
are favored. The configuration adopted in this work is the following: m — 5, n ~ 2, a ^ 0.4, 
/3 = 0.2, which gives good quality of the results and fast convergence. However, the behavior of the 
network does not depend strongly on the parameters and other configurations are equally good. On 
the other hand, the temperature T determines the speed of the convergence. High temperatures 
means slower convergence but gives the system more chances to scape from local minima than in 
the case of low temperatures. 
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Figure 4.13: (Top) The starting point: All ncm-ons enabled. (Middle) The optimum configuration 
by means of the updating rule in (|4.16p . (Bottom) 2D track segments are found as chains of 
consecutive enabled links. 



Gathered information after 2D track finding 

Both methods introduced above efficiently detect 2D track segments inside a cluster of hits: tracks 
are recognized as chains of hits with smooth transitions (links) between them. This association of 
hits into 2D track structures provides a method to detect interaction vertexes in each wire plane, 
which will be of utmost importance at 3D reconstruction stage. Although track projections are 
not necessarily straight, they are fitted to a straight line which gives logical "support" and allows 
to extrapolate the track segments. In addition, the 2D track detection allows to assign to each 
hit a track identification number as well as topological information related with their surrounding 
hits and the track to which they belong. In such a way, a track bound can be easily defined in 
the same way as was done for cluster bounds (see Sec. 14.3.1^ : The bounds of a track arc those 
hits having one link at most with another hit belonging to the same track (it can have however 
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links with another hits of other tracks). Besides, a track is an unidimensional object which can be 
parameterized hke a curve: every hit occupies a given position in the chain. In Fig. 14.141 (top-left) 
two tracks segments have been found in the cluster: the hits colored in green (blue) belong to the 
same track, in light green (light blue) the bounds of the tracks have been marked as well. 

Hit parameterization after 2D track reconstruction 

The information gathered by the 2D track finding procedure is stored as the following parameters 
in the hit structure: 

• The track identification number. 

• The number of links with other hits within the same track segment. 

• The list of hits to which the hit is linked within the same track segment. 

• The position of the hit in the track chain. 

Track expanding and vertex finding 

This information gathered after the 2D track finding procedure can be used to achieve two im- 
portant tasks: Track expanding and vertex finding. A vertex, in a certain wire projection, can be 
defined as the point where two or more tracks meet. Because it is difficult to resolve the different 
hits in the vicinity of a vertex, in most of the cases, the track segments do not reach to the inter- 
section of the tracks (see Fig. 14.141 (top-left)). This technical issue can be tackled in the following 
way: 

• The fitted track is extrapolated from the bound hits trying to find another track bound in 
its extrapolated path (Fig. I4.lil (top-right)). 

• If such a bound is found, the track is expanded by means of virtual hits until that bound 
hit (Fig. 14.141 (bottom-left)). Virtual hits (represented in pink in Figs. 14.141 bottom-right 
and bottom-left) are only logical and do not provide any calorimetric information, they are 
created to complete the geometrical reconstruction of the tracks. 

The latter procedure completes the 2D track finding approach in a proper way, providing the 
right coordinates for the track segments. Finally, vertexes are easily identified as those border hits 
in the track which have a neighbor which is also a border hit of another track. In Fig. 14.151 the 
track finding method performance is shown: The eight ionizing tracks of a deep-inelastic (DIS) 
neutrino event have been detected, expanded and the interaction vertex have been recognized. 

4.5 3D reconstruction 

The goal of the three dimensional (3D) reconstruction is the determination of the space coordinates 
of the different energy deposition segments produced by the ionizing tracks traversing the LAr 
sensitive volume. Each wire plane constrains two spatial degrees of freedom of the hits, one 
common to all the wire planes (the drift coordinate) and one specific for each plane (the wire 
coordinate). The redundancy on the drift coordinate allows the association of hits from different 
planes to a common energy deposition, and together with the wire coordinates from at least two 
planes, allows the hit spatial reconstruction of the hit. 

The 3D reconstruction algorithm developed for the LAr TPC has to different modes of opera- 
tion: 
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Tracks segments 





Irteractton Vertex 




Figure 4.14: The track finding procedure finds two tracks in the cluster of Fig. 14.81 (Top- left) 
The two found tracks arc highlighted in green and blue respectively. (Top-right) The tracks are 
fitted to a straight line. (Top-left) If the extrapolation of a track connects with the bound hit of 
another one, the track is expanded until that hit. (Bottom-right) Final step, the track segments 
arc properly reconstructed and the interaction vertex is recognized. 



• The first one collects the hits within a cluster and tries to match them with the hits of 
complementary views. 

• The second one does the same, but its search of complementary hits is based on the detected 
2D track segments (Sec. l4.4() of different views: it tries to match first the 2D tracks projections 
and thereafter match the hits belonging to them. 

In the following sections we describe both approaches to the 3D reconstruction. For latter 
usage, we start summarizing the expressions to transform from the wire/drift coordinates system 
into the Cartesian system for the 50L TPC. 
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Figure 4.15: Example of the track finding performance over a DIS neutrino interaction cluster 
with eight prong particles. Two of the tracks are highlighted in green and blue respectively. The 
fitted lines are also shown, as well as the virtual hits (in pink) used to expand the tracks till the 
interaction vertex. 

4.5.1 Transforming wire/drift coordinates into 3D Cartesian coordi- 
nates 

We define the local Cartesian reference frame (see Fig. I4.16P centered at the bottom left corner 
of the electronic readout of the LAr sensitive volume. Because both the wire direction on the 
Induction and Collection planes run orthogonally each other, the transformation to the Cartesian 
reference frame is very simple. Thus, coordinates x, y and z are associated directly to the Induction 
wires (in decreasing number direction), to drift time and to the Collection wires (increasing number 
direction) respectively. We explicitly quote below these transformations: 

X ^ {N^ires - 1 - {Wind + 0.5)) p 

y = s (vd/f) (4.18) 
z = {wcol + 0.5) p 

where the constants entering the formulae are the following: A^Mures (= 128 wires), the number of 
read wires in Induction plane; p {— 2.54 mm), the wire pitch at LAr temperature; / (= 2.5 MHz), 
the readout sampling frequency; and Vd (= 0.91 mm/fis), is the electron drift velocity (see 
Sec. I5.2.ip . Finally, Wmd, Wcoi and s are related to the indexes of the wires recording the hit 
in Induction and Collection wire planes, and the drift time coordinate (in number of samples), 
respectively. 
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Figure 4.16: Sketch of the SOL TPC and the Cartesian local reference frame. 
4.5.2 3D reconstruction algorithm hit by hit 

The hit 3D reconstruction algorithm takes as input the sample of hits belonging to any recon- 
structed cluster from an arbitrarily selected view (called main view), and searches for hits from 
complementary secondary views corresponding to the same ionizing track segment. Such a hit is 
called associated hit. The combined information from a hit and its associated hit from one or more 
complementary views allows the determination of the hit spatial coordinates. 

The algorithm proceeds in three steps (passes) . In pass I, all possible associations between hits 
from the main and the complementary views are carried out. The associations are established based 
on the redundancy of the drift coordinate and the relations between hits within the clusters (links). 
Pass II searches for deficient associated hits, redefining their position based on the neighboring hits. 
Finally, in pass III, the optimal spatial coordinates for unassociated hits are computed. 

Pass I 

The goal of pass I is to identify the maximum possible number of associations for hits from the 
main view, combining single-view (links) and multi-view (associations) relations between hits. 

The first and most important step is the search for associations of the cluster bound hits (i.e. 
hits with at most one link to the cluster). This provides a starting reference point for the association 
of the remaining hits of the cluster. Cluster bound hits are required to be associated to cluster 
bound hits from the complementary view. The remaining hits are looped over following the cluster 
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links. In first approximation, a given hit hi linked to a hit /12 is associated to a hit from the 
complementary view only if is linked to a hit /12 associated to If the previous condition 
cannot be fulfilled, hi is associated to the hit association candidate h^ closest to h^. 

The hit association candidates are those hits from a complementary view with equal drift 
coordinates (within a tolerance of 5 drift samples) and compatible wire coordinates. 

Pass II 

The inefficiencies of the hit finding and cluster reconstruction algorithms (see sections |4 . 1 1 and 
together with the presence of ambiguities in the determination of the hit association result in a 
non negligible level of bad hit associations and left-over hits after pass I. 

In order to increase the hit association rate and quality, pass II searches and repairs failures 
on hit associations. For this, we loop over the hits starting from cluster bounds and following the 
links within the cluster. The following algorithm is applied: 

• For every two consecutive (linked) hits, ha and /if,, having both an associated hit in a com- 
plementary view, the position in space of the next hits are compared to the extrapolation 
defined by the straight line joining ha and /if,. 

• All hits (hi, ...,/iri) found until the first hit hn+i whose distance to the extrapolated line is 
smaller than 1 cm are flagged for an ulterior 3D coordinates (re)assignment. Flagged hits 
are freed when no /i„+i hit fulfilling the previous condition is found before a cluster bound 
is reached. 

• Spatial coordinates of the flagged hits are (re)computed in the following way: coordinate y 
is computed by using the corresponding expression from Eq. ()4.18p ; coordinates x and z are 
interpolated using the straight line parameters deflned by the projection of hb and h„+i in 
the (x, z) plane. 

This algorithm repairs in an efficient way the position of those hits coming from straight line 
or small angle tracks, yielding no effect on hits coming from large angle scattering tracks. A third 
pass is needed in order to reduce to zero the rate of left-over hits. 

Pass III 

Hits remaining unassociatcd after pass II are reconstructed in space during pass III, in order to 
reach the required 100% efficiency on the hit 3D reconstruction. Pass III algorithm proceeds as 
follows: 

• Starting from every unassociated hit /iq, the hits in the cluster are looped over, following the 
link path in both increasing and decreasing wire coordinate directions. For each direction, 
the loop continues until the first associated hit (reference hit) is found. 

• If reference hits at both sides of /iq (ha and /if,) arc found before the cluster bounds are 
reached, the spatial coordinates of all unassociated hits between them are computed as in 
pass II. 

• If one of the reference hits is not found, the spatial coordinates of the unassociated hits are 
computed using a straight line extrapolation. The straight line parameters are defined by 
the found reference {ha or /if,) and the closest associated hit met after following the links in 
the same direction we used to find the reference hit. 
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4.5.3 3D reconstruction by 2D track matching 

This approach to the 3D reconstruction of hits is completely based on the information provided by 
the 2D track projection finding method. The effort spent in finding the 2D tracks leads to a very 
clear, simple and efficient 3D reconstruction algorithm. It works as follows: 

• Starting from every track segment ti in the main view, the algorithm tries to match the track 
bounds with the bounds of any another track segment in the complementary view 0. 

• If a complementary track is found, both projections arc associated to each other and the 
algorithm begins the matching of the hits within the tracks. 

• Taking advantage of the track segments structure, the matching of the hits is performed 
exploiting the position of the hits within the track [f]; The association of hits from com- 
plementary tracks is based on the drift coordinate like in Sec. 14.5.21 and on the relative hit 
position in the track. Given a hit hi belonging to a certain track ti in the main view, its 
associated hit hi is searched from the hits in the associated track if in the following way: 
All the hits in t^ whose drift time is compatible with the drift time of hi are collected; if 
there are more than a single candidate, the hit /if whose relative position in if is closer to 
the relative position of hi in ti is selected. 

In Fig. 14.171 we show an example of how this algorithm works on a low multiplicity neutrino 
event: Tracks projections are correctly determined in both views and then matched by time coin- 
cidences of their borders. 

This simple algorithm has an excellent performance when the track projection segments are 
correctly associated. For the rest of the cases where any track segment can not be properly matched, 
the 3D reconstruction is performed using the hit by hit method (see Sec. 14.5. 2|) . 

Hit parameterization after 3D reconstruction 

The picture of the hit is completed by adding the following information to the hit parameterization: 

• The associated hit(s) from the complementary view(s) and its complementary track segment 
(if any). 

• The space coordinates, x, y and z in the Cartesian reference frame. 

4.5.4 Reconstruction performance 

The 3D reconstruction of hits performs efficiently on events of low multiplicity composed of several 
well-defined tracks (not necessarily straight). Since it is purely based on geometrical criteria, 
the reconstruction method is not meant to be carried out on events with more complex energy 
deposition patterns, such as electromagnetic or hadronic showers. For these kind of events, a 
method combining spatial and calorimetric relations among hits needs to be developed. 

This section is essentially dedicated to present some examples of the performance of the spatial 
reconstruction algorithm operating on some neutrino events collected during the running of the 
50L TPC exposed to the CERN neutrino beam (Chapter [5]). The 3D reconstruction of a i/^ CC 

^Obviously, the different projections of the same ionizing track must have equal drift coordinates of their bounds 
(up to a difference of 5 drift samples). 

^ Every track has a structure defined by the chain of hits and every hit occupies a certain position in that chain 
(see Sec. I4.4.3|l . 
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Figure 4.17: Example of the 2D track matching procedure over a low-multiphcity neutrino inter- 
action. The tracks projection arc correctly found and matched. 

event with four charged particles in final state is shown in Fig. 14.181 The tools developed in this 
work have been successfully applied to more complex topologies, as the case of the CC event 
shown in Fig. l4.19| where a deep-inelastic neutrino event with eight charged particles in final state 
is fully reconstructed. 
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Figure 4.18: Example of the 3D reconstruction of a 4-prong CC event. (Top) The two raw data 
views of a z^^ CC interaction which breaks a nucleon into four visible charged particles. (Bottom) 
Three dimensional view of the reconstructed event. 
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Figure 4.19: Example of the 3D reconstruction of a multiprong deep-inelastic I'f^ CC event. (Top) 
The raw data from Collection and Induction views. (Bottom) The three dimensional trajectories 
of the eight primary particles are successfully recovered. 
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4.6 Calorimetric reconstruction 

A proper spatial reconstruction of the tracks together with the information extracted from the 
fitting of the hits ahows to precisely measure the energy deposited by a particle along its path. 
This capability of LAr detectors is of utmost importance since allows the energetic reconstruction 
of every stopping particle, as well as its identification based on the analysis of its energy loss 
pattern. 

4.6.1 The average energy loss (A) 

The average particle energy loss —dE/dx is usually given in units of M eV / {g / crn^) . The length 
unit dx (in g/cm^) is commonly used, because the energy loss per surface mass density dx ~ p ■ ds 
with p (in g/cw?) being the density of absorbing material and ds the length (in cm) is largely 
independent of the properties of the material. This length unit dx consequently gives the surface 
mass density of the material. 

According to the Bethe-Bloch formula, the average energy loss is given by the following Eq. [SH 

US]: 



- 2Tr NAr^rUeC^ 
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(4.19) 



with 



iiiNArlm^c^ = 0.1535 MeV cmVg; 

z: charge of the incident particle in units of e [z^l for muons); 

Z,A: atomic number and atomic weight of the absorber (Z=18, A=39.948 for argon); 
me,re: electron mass (0.511 McV) and classical electron radius (2.817 • 10"^'^ cm); 
Na: Avogadro number (6.022 • lO^^ mor^); 

/: Ionization constant, characteristic of the absorber material (188 eV for argon); 



- /? = w/c of incident particle (7 = 1/ \/l — P'^) ; 

- p: density of absorbing material (1.4 g/cm'^ for liquid argon); 

- Wmax- maximum energy transfer in a single collision (see below); 

- S: density correction (see below); 

For all practical purposes, Eq. (j4.19p in a given material is a function only of (3 and only a 
minor dependence on toq (the mass of the ionizing particle) at the highest energies is introduced 
through Wmax ■ The expression as computed for muons in liquid argon is shown by the top curve 
in Fig. 14.201 exhibiting a minimum of 2.1 MeV/cm around = 360 MeV. 

The maximum transferable kinetic energy to an electron depends on the mass and the momen- 
tum of the incident particle. It can be probed that for relativistic particles {Ekin ^ E and p ^ E, 
with c = 1): 
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Wmax (7) 
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(4.20) 



Values for the density correction S can be approximated by [84j : 
X <Xq 

4.6052 X + C + a{Xi - X)"' Xa < X < Xi 
4.6052X + C X > Xi 

where X = logio{f3^). The quantities Xo,Xi,to, a and C depend on the absorbing material. For 
liquid argon we assume the following values [51]: Xo=0.201, Xi=3, m=3, a=0.196 and — C=5.212. 

It is also of interest to consider the mean energy loss excluding energy transfers greater than 
some cutoff T^ut {truncated mean energy loss) [S]: 
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so that the equation reproduces the Bethe- 
In practice, the cut-off Tcut is given by the energy threshold for 



upper 

Bloch equation for Tcut > W, 
S-TSty detection. 

Finally, in practical situations, the mean or average value is given by 
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(4.22) 



where dE/dx comes from the Bethe-Bloch formula (Eq. (|4.2ip ). and Ax is the distance within we 
are measuring the energy loss. 



4.6.2 The most probable energy loss (Amp) 

The energy loss distribution for thin absorbers is strongly asymmetric, due to the intrinsic statistical 
nature of the ionization process. The fluctuation of the energy loss around the mean value are 
well described by the Landau distribution [ST] . According to the Landau's theory, the probability 
/l(x. A) that a singly charged particle loses an energy A when traversing an absorber of thickness 
X is expressed as: 



/L(x,A) = i.0(A) 



where 
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X is the deviation from the most probable energy loss Amp ■ 
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where A is the actual energy loss in a layer of thickness x and 

^ = 0.1536 • f ^ • P a; [MeV] 

{p in g/cm^ and x in cm) 



(4.23) 
(4.24) 

(4.25) 
(4.26) 
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For the particular case of LAr TPC's, x (= Az) is the effective portion of the track exposed to 
the wire, also called track pitch length. 

The most probable energy loss is calculated to be [87] : 
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where je = 0.577 is Euler's constant. 6 is the density effect and: 
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(4.28) 



Therefore, for a given material the most probable value energy loss only depends on the muon 
momentum and on the track pitch length. Eq. (|4.27p represents a very good approximation for 
the most probable energy loss of charged particles in thin absorbers on the vicinity {Pj ^ 4) and 
beyond the "minimum-ionizing region". In practice, A„ip is extracted by fitting the energy loss 
distribution with a Landau-Gaussian convoluted function, where the Gaussian accounts for the 
resolution in the charge measurement, mainly due to the electronic noise. The Landau-Gaussian 
convolution is defined by the following expression: 
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(4.29) 



where K is an overall normalization factor depending on the statistics, a is the width of the con- 
voluted Gaussian distribution, e the energy and is the Landau distribution, given by Eq. (|4.23p . 



4.6.3 Average vs. most probable energy loss 

In many practical situations, is necessary to parameterize the measured charge of a certain ionizing 
particle. In principle, we could do that by means of either A or A^p, however, the use of the mean 
value presents two technical problems, one intrinsic to the definition of this quantity, the other 
related to the limited available statistics. On one hand, the fluctuations of the energy loss are 
determined by a few high-energy 5-rays, and thus A depends on the (5-ray detection threshold (Tcut) 
entering expression (|4.2ip . This quantity is not a priori known since it has a strong dependence 
on the performance of the spatial reconstruction algorithm. On the other hand, since A strongly 
depends on the Landau tail, large statistics are needed in order to significantly populate the tail 
and hence obtain a reliable estimation of A. None of the previous considerations affects the 
determination of Amp- 

On the other hand, the use of A^p as energy loss estimator presents an additional problem: 
while A is linear with the crossed distance (Ax), as shown in Eqs. (|4.19p and (|4.22p . A^p behaves 
as AxlnAa;, as seen from Eqs. (|4.27p and (|4.26p . This means that the energy loss per crossed 
distance measured at hits from track segments of equal kinetic energy but different associated 
crossed distances yields Landau distributions with equal A but different A^p values. The effect 
is shown in Fig. 14.201 where the mean and most probable energy losses per crossed distance are 
plotted as a function of the muon momentum for various Az values. 

4.6.4 Particle identification 

The function (|4.19p as computed for protons, kaons, pious and muons in liquid Argon, is shown 
by the curves in Fig. 14.211 The lines follow the expected analytical "average" behaviour only. The 
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Figure 4.20: Average energy loss —dE/dx from Bethe-Bloch formula (top line) compared to the 
most probable energy loss for three values of the track pitch length Az (bottom lines). 



experimental points spread around the corresponding curve for each type of particle because of the 
fluctuations of the energy loss by catastrophic energy loss processes, i.e. by interactions with high 
energy transfers, and because of the detector resolution on the kinetic energy and dx measurement. 

At non-relativistic energies, dE/dx is dominated by the overall 1/0^ factor and decreases with 
increasing velocity until about k, A (or (3 « 0.96), where a broad minimum of ionization is 
reached (a particle at this point is known as a mip). 

From the figure, it is clear that for energies below the minimum ionizing value, each particle 
exhibits a dE/dx curve which, in most cases, is distinct from the other particle types because of 
the mass difference. This is precisely the characteristic that has been exploited in this work as a 
means for identifying particles in this energy range. 

The minimum value of dE/dx is almost the same for all particles of the same charge and, as 
the energy increases beyond this point, the term \/0^ becomes almost constant and all curves 
seem to converge together to very similar dE/dx values. As soon as the detector resolution on the 
measured kinetic energy and on the track pitch are taken into account, the fact that all particles 
behave so close above the m.i.p. region independently on the mass, make them indistinguishable 
by means of the dE/dx technique. 



72 



Chapter 4. Off-line event reconstruction 




4.6.5 Energy calibration 

The energy calibration consists of transforming the charge associated to a hit into energy units. 
The ionization charge is precisely measured in the Collection wire plane (see Sec. I4.2p . The energy 
Ei associated to a given hit is: 

CW 

= — Qm (4.30) 

where C is the calibration factor; W is the average energy needed for the creation of an electron-ion 
pair (i.e. 23.6 eV); is the electron- ion recombination factor and Qi{t) is the corrected charge 
by the electron lifetime from Eq. (j3.ip in Sec. 13.21 The calibration factor C converts the detector 
measuring units (ADC counts) into charge units (fC), while the recombination factor TZ depends 
on the absorber medium, on the applied electric field and on the density of released charge, i.e. on 
dE/dx. 

Thus, the energy calibration consist on determining C and finding a suitable model which 
describes the recombination factor TZ for a given detector configuration. 



The electron recombination factor 

Ionization electrons produced by the passage of a charged particle through the LAr volume have 
a non negligible probability to recombine with one of the positively ionized atoms. In such a 
case, the energy released in the ionization appears in the form of photons and we say that the 
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drift electrons have been quenched. The recombination probabihty depends on the absorber, on 
the appHed electric field, and on the density of released charge, i.e. on dE/dx. This dependence 
may be intuitively understood as follows: the higher the field the faster the electrons and ions are 
separated, and the less time is available for the recombination; on the other hand, the larger the 
dE/dx, the higher the number of neighbouring positive ions which are available for an electron to 
recombine. 

The dependence ofTZ with dE/dx can be modeled by Birk's law [S^ which has been successfully 
applied to LAr detectors in [5^ : 

dx 1 + ks ^ ' 



dx 



dE/dx being the pure Bcthc-Bloch ionization loss, a and fc^ are constants depending on the applied 
electric field and fcs is the Birks coefficient accounting for quenching due to recombination. The 
recombination factor is the ratio between the measured and the produced (theoretical) charge: 



AQ 



nieas 



^ = ^TWr (4-32) 



Therefore, Eq. (j4.3ip can be rewritten in terms of the recombination factor by noting that a 
has units of charge per energy, so that it can be expressed as the inverse of the factor CW (see 
Eq. (|4.30p ) times a dimensionless multiplicative constant (5 = aCW) ranging between and 1. 
Multiplying both sides of the equation by CW and dividing by dE/dx we obtain: 

A precise knowledge of the recombination factor is of the utmost importance since it determines 
the absolute energy calibration of the detector. 

In Sec. 15.2.^ we present the analyses for the energy calibration for the SOL LAr TPC. 



4.6.6 Measurement method 

From the previous discussions, it is clear that the determination of the energy loss per crossed 
distance can be obtained by means of the measurement of the deposited energy for every hit 
together with the calculation of its associated track pitch length. In the following we discuss some 
aspects of the procedure to achieve this task in the framework or LAr TPC detectors. 



The track pitch length 

The length of the track portion associated to a hit hi (hereafter called the track pitch length: Azi) 
depends on the wire pitch and the orientation of the track with respect to the wire plane. It can 
range from the wire pitch length (for track pieces parallel to the wire pitch direction) to distances 
of the order of the chamber's side (for tracks perpendicular to the pitch direction). Thus, the track 
pitch length can be written in terms of the track direction in a given hit: 

Az, = ^ (4.34) 
cos 0, 



being p the wire pitch and 0i the angle between the track at the hit hi position respect to the 
direction of the wires (see Fig. 14.22"]) . 
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Track segmentation 

In order to make an estimation of the track pitch length at any given hit hi belonging to a track, 
a partition of the track into segments of a certain number of hits is performed. Since a track is 
not straight in general, the aim of this procedure is to obtain a parameterization of the trajectory 
of the track in terms of straight line segments, which locally approximate the real track direction. 
Given a certain track segment with n hits belonging to it, its direction is determined by fitting 
a three dimensional straight line to these n points. Besides, an additional constraint is imposed 
to this fit: any track segment uses the result to the fit of the previous calculated one in order to 
have smooth transitions between the different parameterized lines, i.e. the fitted lines belonging 
to consecutive segments must have one point in common. As a result, we obtain a partition of the 
track trajectory into straight segments, which are consecutively connected defining a trajectory 
(see Fig. 11211). 




Figure 4.22: Illustration of the track segmentation procedure and the track pitch determination. 

The direction associated to any track segment is determined by the parameters entering in the 
linear fitting, given in terms of a unit vector d;, or in terms of the angles 9 and tj) 

9i = arccos((i^) 

(f,, = arctan (^^^ (4.35) 

where the index i may refer to a hit or a segment, and = {df , c?|, c?f } is the corresponding 
associated direction vector. 

Hence, the energy deposition of a given hit is calculated from the fitted area in Sec. 14.2.11 while 
its associated pitch length is obtained from the track segmentation procedure through Eq. (|4.34p . 

Finally, the mean energy loss per crossed distance (Eq. ()4.19p ) measured for every hit is ap- 
proximated by: 

/ dE\ E, 
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with Ei and Azi the energy (Eq. (|4.19p ) and track pitch (Eq. (|4.34p ) associated to hit hi. 



Range determination 

The segmentation procedure also provides an optimal method to compute the track ran(;^. From 
a theoretical point of view, the mean range for a stopping particle might be calculated through the 
integration of the Bethe-Bloch formula (|4.19p along the particle's path: 

dE 

Jq dE/dx ^ ^ 

where integrations limits go from the initial point with kinetic energy Eq until the stopping point. 
If one assumes that the energy loss is continuous, the distance a particle can penetrate before it 
loses all its energy is a well defined number, the same for all identical particles with the same initial 
energy. Fig. 14.231 shows the energy-range curves for different particles calculated by a numerical 
integration of the Bethe-Bloch formula, which yields the approximate path length traveled. As in 
the dE/dx vs. kinetic energy case (see Fig. l4.21( ). the curves give the expected "average" behavior. 




Figure 4.23: Kinetic energy vs. range from Bethe-Block. 



The definition of a range is complicated because of the multiple Coulomb scattering in the 
material and catastrophic energy loss processes (both ignored in Eq. (|4.37p ). which cause the 
particle to follow a zigzag path through the absorber and lead to substantial range straggling. 



®The distance a particle can penetrate before it loses all its energy. 
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The range, defined as the straight-hne length, will generally be smaller than the total zigzag path 
length. The track segmentation method leads to a better estimation of the "real" range: 

Ns 

i? = ^ (4.38) 

being the straight-line length of each one of the segments composing the track, and Ng the total 
number of such segments. 



Chapter 5 



The SOL TPC exposed to The 
CERN Neutrino Beam 

Within the broad physics programme of the ICARUS project, an important issue is the study of 
oscillations through the detection of accelerator neutrino interactions. The possibility to observe 
ly^ — > Ve and — *■ Vt by means of kinematic criteria is known to be limited, among others, 
by the knowledge we have of nuclear effects (Fermi motion, nuclear re-scattering and absorption, 
etc.). Therefore it is very important to acquire experimental data in order to tune the existing 
Monte-Carlo models. In addition, the exposure of a LAr TPC to a neutrino beam is mandatory to 
demonstrate the high recognition capability of the technique and gain experience with real neutrino 
events. 

In 1997, the ICARUS collaboration together with a group from INFN and Milano University [90] 
proposed to expose a 50 liter LAr TPC to the multi-GeV wide band neutrino beam of the CERN 
West Area Neutrino Facihty (WANF) [H], during the NOMAD ^ and CHORUS [53] data taking. 
The test was part of an R&D program for a medium baseline Vr appearance experiment |94j . The 
idea was to collect a substantial sample of quasi-elastic (QE) interactions (t'^ -I- n — > p -I- to 
study the following physics items: 

• Measurement of the acoplanarity and missing transverse momentum in events with the ^-p 
topology in the final state, in order to assess Fermi motion and proton re-scattering inside 
the nucleus. 

• Appearance of nuclear fragments (short tracks and blobs around the primary interaction 
vertex) in quasi-elastic events. 

• A preliminary evaluation of e/7 and e/n'^ discrimination capability by means of the specific 
ionization measured on the wires at the beginning of the candidate track. This measurement 
is limited by the size of the chamber. 

The data collected in 1997 offer the unique opportunity to study nuclear effects in the Ar- 
gon nucleus and to assess the identification and reconstruction capability of a LAr TPC for 
low-multiplicity neutrino events. The work presented in this thesis shows, for the first time, a 
comprehensive set of results from the 1997 test (see [SSI [Ml [IZl [HS for preliminary results). The 
structure of the chapter is as follows: the experimental setup and the calibration measurements of 
the 50 liter TPC are discussed in Sec. 15.11 and Sec. 15.21 respectively. At the beginning of Sec. 15.31 
the event reconstruction and particle identification are detailed, thereafter the analysis of a golden 
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sample of quasi-elastic charge current (CC) interactions is presented together with a comparison 
with theoretical expectations. Finally, in Sec. 15.41 we measure the quasi-elastic cross section for 
CC interactions. 



5.1 The experimental setup 

The LAr TPC was placed on a platform 4.5 meters above ground, right in between the CHORUS 
and NOMAD detectors (Fig. 15. ip . Full details of the geometrical and technical features of the 
50L TPC have been treated in Sec. 13.11 As mentioned there, the active part of the detector is 
located inside a stainless steel cylinder. The connection to the outside area is obtained through 
a set of UHV flanges housing the signal, the high voltage cables and the vacuum feed-through. 
The cylinder is positioned into aim diameter dewar partially filled with low purity Liquid Argon 
acting as a thermal bath and it is rotated 30° along the vertical axis with respect to the nominal 
beam direction to reduce the number of particles crossing just one readout wire. 

The modest size of the LAr TPC fiducial volume (~50 liter), coupled with the high energy of 
the WANF V beam (Sec. I5.1.ip . made necessary a muon spectrometer downstream the TPC. The 
NOMAD detector was the perfect choice to achieve this task: A coincidence with the NOMAD 
DAQ was set up to use the detectors located into the NOMAD magnetic dipole as a magnetic 
spectrometer (Sec. 15.1^ . The experimental setup was completed with additional counters for the 
trigger and veto systems (Sec. I5.1.'5| . 



ML TPC & NOMAD aatup 




Figure 5.1: A sketch of the experimental setup (top view). 



5.1. The experimental setup 



79 



5.1.1 The Neutrino Beam 

The chamber has been exposed to the v beam produced at the CERN West Area Neutrino Facihty 
(WANF [UlIM])- The geometry of the WANE is shown in Fig.ES The Super Proton Synchrotron 
(SPS) at CERN suppUed a beam of 450 GeV protons which was directed onto a berylhum target, 
thereby producing a large number of secondary hadrons. In every accelerator cycle (14.4 s) the 
protons used to produce neutrinos were extracted from the SPS in two spills, which were separated 
by 2.7 s. Each spill had a full width at half maximum of 3 ms and contained about 1 .8 x 10^^ protons. 
The number of protons incident on the target was monitored by a beam current transformer (BCT). 
The number of produced neutrinos was often expressed in units of protons on target (pot). 

TOP VIEW of neutrino cave 
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Figure 5.2: Schematic layout of the WANF beam line (all lengths in meters) 



The multiplicity and charge asymmetry of the secondary particles were measured by secondary 
emission monitors (SEMs) placed before and after the target. The beam focusing elements, down- 
stream of the target, were large magnets, designed to collimate the positively charged secondary 
particles, while deflecting the negatively charged particles away from the beam axis. Most of the 
focusing was done by the first element, called the horn. The second element, the reflector, focused 
the smaller angle particles which had been missed by the horn. Two helium bags were installed 
in order to increase the neutrino flux by reducing multiple scattering and secondary interactions 
along the beam. Downstream of the focusing elements was a vacuum decay tunnel of 290 m length 
in which the positively charged particles could decay and, thereby, produce the neutrinos. The 
major contributions to the flux came from the reactions: 

7r+ ^+ + 
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The charge conjugate reactions also occurred, giving rise to a component in the beam. Earth and 
iron shielding directly after the decay tunnel filtered out all but neutrinos and some muons. Silicon 
detectors in several pits within the shielding region measured the muon multiplicity and could be 
used to monitor the beam and determine the absolute flux. An additional toroidal magnet, which 
was placed after the muon pits, deflected most of the remaining muons. The resulting neutrino 
beam reached the NOMAD detector 835 m downstream of the proton target with the composition 
and properties shown in Fig. 15.31 and Tab. 15.11 [99] . The predictions of the neutrino fluxes are 
based on Monte Carlo simulations. For further reference on the WANF layout and the alignment 
procedures, see |100j . 




50 100 150 200 50 100 150 200 
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Figure 5.3: Composition of the neutrino beam and the v^, v^, and Ve energy spectra at NOMAD, 
within the transverse fiducial area of 260 x 260 cw? . The Figures are taken from 



The mean energy of the reaching the detectors is 24.3 GcV when integrating over the active 
NOMAD area (2.6 x 2.6 m?), while contamination from other flavors are below 7% for and 
~ 1% for [55] ■ Being the TPC centered on the beam axis and covering a smaller surface, a 
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Tabic 5.1: Composition of the neutrino beam and the sources of the different components 



harder neutrino flux is expected, with a mean energy of 29.5 GeV (see Fig. 
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Figure 5.4: energy spectra on SOL (NOMAD) transverse fiducial area. 



5.1.2 The NOMAD detector 

The NOMAD detector [32] was composed of several independently working sub-detectors. A 
schematic view of the detector is shown in Figure 15.51 along with the coordinate system. The 
sub-detectors were placed along the beam-line with positional sensitivity transverse to the beam 
direction. The veto system upstream of the detector allowed a fast identification of muons entering 
the detector. Following the muon veto, there was the forward calorimeter (FCAL) which was 
used as a target for studies requiring very high statistics while sacrificing some resolution in the 
vertex finding. The drift chambers (DCH) in the center of the detector served both as an active 
target and a very precise tracking device; hence the most interesting of the detected interactions 
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occurred there. Further downstream, the trigger planes ^ the transition radiation detector (TRD), 
the preshower (PRS) and the electromagnetic calorimeter (ECAL) were situated. These inner 
detectors (DCH to ECAL) were all surrounded by the old UAl dipole magnet which provided a 
constant magnetic field of 0.4 Tesla. The hadronic calorimeter (HCAL) followed downstream of 
the inner detectors and at the very end there were the muon chambers. The x and y axes of the 
co-ordinate system were centered in the middle of the first drift chamber. The neutrino beam 
entered with an angle of 42 mrad with respect to the z axis. 

In what follows, we describe some of the NOMAD subdetectors that were key for the experiment 
and the analyses described in this work. 




Figure 5.5: Side view of the NOMAD detector. 



The Veto Counters 

Being the most upstream detector, the veto consisted of 59 scintillator counters and covered an 
area of 5.4 x 5 m? . It was used to detect muons traveling with the neutrino beam and charged par- 
ticles from upstream interactions, and prevented them from causing valid triggers. An additional 
veto plane was positioned between the FCAL and the DCH. Charged particles from interactions 
downstream of the veto, e. g. in the magnet coil or the support structures of the drift chambers, 
were removed during analysis by constraints on the vertex position. 

The scintillators were read out by photo-multipliers and a veto signal was constructed from a 
logical OR of all the individual counters. The efficiency of the veto was constantly monitored and 
found to be stable at about 96 % — 97 %. The contribution of the veto to the overall dead time of 
the experiment during the two neutrino spills amounted to about 4 %. A detailed description of 
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the veto system can be found in [lOlj . 
The Forward Calorimeter (FCAL) 

The forward calorhiieter, being the instrumented front pillar of the magnet, was used as a massive 
active target and allowed neutrino physics studies including di-muon production and neutral heavy 
lepton searches. The FCAL consisted of 23 iron plates, which were 4.9 cm thick and separated by 
1.8 cm air gaps. The first 20 gaps were instrumented with scintillator counters. The dimensions 
of the scintillators were 175 x 18.5 x 0.6 cw?. Five consecutive counters along the beam direction 
were bunched together via light guides to form a module which was read out at both ends by 
photo-multipliers (see Figure 15. 6p . Ten such modules were arranged vertically to form a stack. 
There were four stacks aligned along the beam axis. The active region of the FCAL had a mass 
of 17 tons, was about 5 nuclear interaction lengths in depth and had an active area transverse to 
the beam of 175 x 190 cm?. Further details can be found in 



V -beam 




175 cm 



Figure 5.6: Top view of the FCAL. 



The Drift Chambers (DCH) 

A crucial part of the detector were the drift chambers which provided at the same time the target 
material and the track and vertex reconstruction of the charged particles. The full drift chamber 
region consisted of 11 modules, each made up of four individual drift chambers (composed of low 
density, low atomic number materials). In addition to these 44 drift chambers, there were another 
5 chambers in the TRD region, providing additional tracking. The chambers were filled with a 
mixture of argon (40 %) and ethane (60 %), and equipped with sense wires. The volume with 
maximal acceptance and reconstruction efficiencies was referred to as fiducial volume and had the 
dimensions of 2.6 x 2.6 x 4 m? with a mass of 2700 kg. The overall density of 100 kg/m? was about 
equal to the density of liquid helium. 

The sense wires in the three gaps of a drift chamber were oriented with -5, and -1-5 degrees with 
respect to the magnetic field direction, corresponding to the NOMAD x axis. The three ionization 
signals produced by a particle traversing a drift chamber allowed a hit positioning accurate to 
about 1.5 mm in the x direction, due to the angles between the wires. Transverse to the wires {y 
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direction) a position accuracy of 150 was reached by the separation of the wires and a constant 
gas drift velocity. The position along the beam line (z direction) was determined by the position 
of the wire planes which was constantly monitored. The efficiency of each wire to record a hit 
was typically 97 %. The three dimensional information from each drift chamber could be used to 
reconstruct the helical path of charged particles through the magnetic field. The curvature allowed 
the measurement of the momentum of the particles through the relation 

pcosOx = qBr (5-1) 

where p is the momentum and q the charge of the particle, 9x the pitch angle of the helix, r the 
radius of curvature and B the magnetic field. The momentum resolution was a function of the 
particle momentum and track length. For muons and charged hadrons with normal incidence to 
the measuring planes the momentum resolution could be parametrized by 



with 

a = 0.05 to1/2 

b = 0.008 m5/2(Gey/c)- 



(5.2) 
(5.3) 

and L denoting the length of the track. The first term in Eq. (|5.2p corresponds to the error due to 
multiple-scattering, the second term arises from the single hit resolution. A detailed description of 
the drift chamber is given in |102j . 



5.1.3 The Trigger and Veto system for the 50 liter set-up 

The trigger is provided by a set of scintillators located between the chamber and the NOMAD 
apparatus (Fig. 15. ip . Each of the 3 trigger scintillator counters has a dimension of 110x27 cm^. 
They were positioned in a half-circle 60 cm beyond the center of the chamber. Incoming charged 
particles are vetoed by 5 scintillators mounted in front of the chamber, 50 cm before its center, and 
by the last scintillator plane of CHORUS. The latter vetoes particles deflected by the CHORUS 
magnetic field entering the chamber at large angles with respect to the nominal beam direction. 

The trigger requires the coincidence of the SPS beam spill and at least one of the trigger 
scintillators, vetoed by the scintillators put in front of the chamber and the CHORUS plane. This 
trigger is put in coincidence with the two trigger scintillator planes of NOMAD (Tl and T2 in 
Fig. 15. ip jl01| . Moreover, a trigger is rejected if the NOMAD acquisition system is in BUSY mode 
or if the delay with respect to the previous trigger is lower than 500 /iS ( "drift protection" 

The request of a charged particle triggering the chamber locally and reaching NOMAD up 
to Tl and T2 inhibits the acquisitions of neutral-current and i>e charged-current events, limiting 
the sample of the present test to i'^ charged-current interactions. The trigger efficiency has been 
monitored during data taking through a dedicated sample of through-going muons. It turned out 
to be 97% averaged over the whole data taking period. 



This condition inhibits the occurrence of event overlaps in the multibufFer readout system. 
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Figure 5.7: Measurement of the electron drift velocity. 



5.2 Calibration of the TPC 

5.2.1 Measurement of the electron drift velocity 

A fundamental parameter for the reconstruction of ionizing events occurring inside the fiducial 
volume of the 50L TPC is the velocity of the drifting electrons. 

An absolute determination of the drift velocity and its uniformity along the fiducial volume of 
the chamber has been obtained exploiting the external scintillators and the additional information 
from NOMAD. A dedicated trigger selecting through-going muons ("mip sample") has been put 
into operation adding two additional scintillators (30x30 cm^) before and after the external dewar. 
The relative position of the chamber with respect to the NOMAD reference frame has been obtained 
by residual minimization of the track parameters recorded both by the TPC and the NOMAD data 
acquisition system. The drift velocity is obtained by fitting the absolute vertical position of the 
muon as reconstructed by NOMAD trackers versus the drift time measured locally by the TPC 
fFig. 15.7)1 . The corresponding drift velocity is 



Vd = 0.905 ± 0.005 mmliis 



(5.4) 



The measurement of the drift velocity is essential because it defines the "pitch" in the drift coor- 
dinate: Since the signal is sampled with a 2.5 MHz frequency for a duration of 820 pjs, (2048 time 
samples), the drift pitch amounts up to 0.36 mm. The maximum drift distance of 47 cm is covered 
by a drifting electron in 520 
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Side View Top View 




Figure 5.8: Example of the 3D reconstruction of parallel through-going mips. 



5.2.2 Measurement of the electron drift lifetime 

The attenuation of the drifting electric charge towards the anode is a crucial working parameter 
for the TPC. It essentially depends on the purity levels in the Liquid Argon and on the applied 
electric field. Thanks to a purity monitor placed within the pure LAr bath (see Sec. 13. 2p . it was 
possible to obtain direct measurements of the drifting electron lifetime t. During the data taking 
period the impurity levels of oxygen kept themselves stable, getting values of r higher than 10 ms. 

In this section, we perform an independent measurement of r based on the spatial and calori- 
metric reconstruction of a sample of muons crossing the chamber at different heights (sec Fig. l5.8p . 
This sample of muons was gathered from the data collected during the exposition to the WANF 
beam, by means of an automatic scanning and reconstruction algorithm. Over the whole data, the 
algorithm was able to collect a big enough number of correctly reconstructed through-going muons. 
They are essentially muons coming from i/^ CC interactions occurring before the fiducial volume of 
the TPC, uniformly distributed in height and crossing the chamber with an approximately constant 
value of their drift coordinates. 

The spatial and calorimetric reconstruction tools allows to measure the energy loss per crossed 
distance (dE/dx) hit by hit for every muon (see Sec. 14.61 and Eq. (|4.36p ). The most probable value 
(Amp) is used as a suitable estimator of the energy loss for each one of those muons (see Sec. 14. 6. 2| . 
As it was already discussed in Sec. 14.6.31 Amp depends on the crossed distance associated to the 
sampling hits of the track, but shows a slight dependence on the kinetic energy. Being the sample 
of mips used in this analysis horizontal and straight muons coming in the beam direction, their 
track pitch length hardly varies from the mean value: Az = 0.29 ± 0.01 cm. The latter ensures 
that the expected Amp is essentially the same for all the muons in the sample (see Fig. 15. 9p . 

In order to measure the charge attenuation effect due to the LAr impurities, we analyze the 
dependence of Amp as a function of the drift coordinate The procedure is as follows: 
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Figure 5.9: Average and most probable energy loss as a function of kinetic energy. 



• A dE/dx histogram of every muon is obtained. 



• Every histogram is fitted to a convoluted Landau-Gaussian distribution, in order to obtain 
the most probable value of the distribution (Amp)- On the other hand, the average drift 
coordinate (Tmean) of the track is computed as well. 

• A two dimensional distribution of the Amp versus the mean drift time is built. 

• Finally, the latter distribution is converted into the profile form of Fig. 15.101 and then, fitted 
to a straight line. 

From the slope of the curve in Fig. 15.101 is possible to obtain an accurate estimation of the 
electron drift lifetime during the data taking stage. When the drift time t is much smaller than r, 
the exponential attenuation of the drifting charge (Eq. (|3.ip ) can be approximated by an straight 
line, being the slope the inverse of the lifetime. Thus, from the fitted data in Fig. 15.101 we have 
a value for the slope of pi = —0.015 ± 0. 007 MeVfis/ cm, which leads to a lifetime of r = 120 fis 
with a relative error of 50%. Such a big error comes from the flatness of pi. This means that we 
are over the sensitivity limit of the measurement and that the attenuation over the drift distance 
is negligible. 
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Figure 5.10: Most probable value of dE/dx vs drift time. 



5.2.3 Beam alignment 

In order to obtain a correct reconstruction of the events is mandatory to determine precisely the 
position of the chamber with respect to the beam. The alignment has been done moving the local 
reference system (RS) of the TPC (see Sec. I4.5.ip to the NOMAD reference system and then, to 
the beam. Since NOMAD alignment with respect to the beam was very precisely measured, this 
error source is less important than the one we will have for the alignment of the TPC with respect 
to NOMAD. 

We call 9 the rotation in the z — x plane and the rotation in the y — z plane. In terms of 
these angles, the alignment of NOMAD with respect to the beam is written as [52] 

9m ~ mrad 

(5-5 

0JV = 42 mrad = 2.41° ^ ' 



Alignment of the TPC respect to NOMAD 

In order to determine the alignment parameters of the SOL TPC respect to the NOMAD RS, we 
made use of the same mip sample of Sec. 15.2.11 used for the electron drift velocity determination. 
There are 268 through-going muons correctly measured and reconstructed both by NOMAD and 
the LAr TPC. The relative angles 9 and are then calculated through the minimization of the 
difference between the direction angles measured by the 50L and NOMAD respectively. Final 
distributions after the minimization are shown in Fig. 15.111 The value of the alignment parameters 
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Figure 5.11: Difference of the measured SOL and NOMAD angles along x and y axes 



arc determined as the ones which make these distributions to be centered at zero: 

6* = 30.08 ± 0.06° 
= -0.34 ±0.03° 



5.2.4 Energy calibration 

As it was discussed in Sec. 14.6.^ the energy calibration consists in expressing the charge associated 
to a hit in energy units. The measured charge is related to the real deposited energy by means 
of a unit conversion constant iCW) and the recombination factor [TZ in Eq. (|4.30p ). In this 
section, we measure the energy calibration factors for the SOL TPC under the assumption that the 
recombination factor is well described by Birk's law in Eq. (|4.3ip . 



The overall calibration factor 

The calibration factor C converts the detector measuring units (ADC counts) into charge units 
(fC), and is usually obtained from a set of dedicated test pulses, simulating signals ranging from 
single to several tens of minimum ionizing particles [1031 [71] . It was determined for the collection 
wires during the set-up phase before the data taking, when read-out was running in "charge" mode. 
As it was commented in Sec. 13.41 the final (and optimal) configuration adopted for the read-out 
worked in "current" mode which, unfortunately, did not have a calibration test pulse to measure 
C. 

Due to this lack of information, what we measure instead is the so called overall calibration 
factor a = CW/R in Eq. ()4.30p . To this purpose, we have collected a sample of around 3000 
through-going muons coming from CC interactions. The spatial and calorimetric reconstruction 
of those muons allows to measure precisely the deposited charge per unit of length {dQ/dx) as 
described in Sec. 14.61 We show in Fig. IS. 121 the measured distribution of dQ/dx for the 3000 
mouns in the sample: The data is fitted to a convoluted Landau-Gaussian function to obtain the 
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Figure 5.12: dQ/dx of crossing muons. The distribution is fitted to a convoluted Landau-Gaussian 
distribution. 

most probable energy loss for mips in terms of detector charge units: A^p = 216.7±0.1 ADC/cmH 
The most probable energy loss slightly depends on the energy for the range of the muons being 
considered (see Fig. 15. 9p . From this measurement, a can be obtained assuming that according to 
the theoretical prediction for muons in the range of energies of the beam, the most probable energy 
loss is A,„p = 1.736 ± 0.002 MeV/cm, therefore: 

= 8.01 ± 0.01 X 10"''^ MeV/ADC (5.7) 



*-rmp 



cw 

R 



A 



mp 



mip 



A, 



The recombination factor 

As we discussed in Sec. 14.6.51 the recombination factor TZ of drifting electrons strongly depends 
on the ionization density and this dependence of TZ with dE/dx can be modeled by Birk's law 
(Eq. Km \. 

dQ «f 



The size of the range of the fully contained protons collected in the golden sample considered 
in Sec. 15.31 offers the opportunity to precisely test the energy loss pattern in the active medium. 



2 The bar in ^rnp means the value is given in detector charge units. 
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Drift velocity 


0.905 ± 0.005 mm/MS 


Electron lifetime 


> 10 ms (over the sensitivity limit) 


Alignment angles 


= 30.08 ±0.06° (/) = -0.34 ±0.03° 




8.01 ± 0.01 X 10-^ MeV/ADC 




0.035 ± 0.001 cm/MeV 



Table 5.2: Calibration parameters of the 50L detector. 



Fig. l5.13l shows the observed charge per unit length (dQ/dx) versus the residual range for protontH. 
Dots are direct measurements from the reconstructed hits of the proton tracks. From this distri- 
bution we estimate the most probable value of dQ/dx for each bin of the residual range. Finally, 
from this data we fit fc^ = 0.035 ± 0.001 cm/MeV ( (|4.3ip ) for the proton hypothesis (continuous 
line), constraining the value of a to make TZ compatible with Eq. (j5.7p at mip region: 

a = (1 + kBAmp\mip)/armp (5.8) 

Based on the knowledge of the Birk's law parameters, a precise determination of the deposited 
energy can be obtained through direct charge measurements by means of Eq. (|4.3ip . Moreover, pure 
Bethe-Bloch ionization loss can be corrected by the detector response in the region of high dQ/dx, 
where quenching effects are more sizable. This correction is important for particle identification in 
LAr, which exploits the different energy loss of particles near the stopping point. 

In summary, we have precisely measured all the calibration parameters needed for a proper 
kinematical reconstruction of the i' events occurring in the TPC. The electron drift velocity allows 
to translate time samples into the vertical coordinate ?/, and the electron lifetime is needed to 
correct the attenuation of the measured charge. The alignment angles are very important to 
properly describe the kinematics on the transverse plane respect to the beam direction. Finally, 
the recombination factor precisely describes the calorimetric picture of the events. In Tab. 15.21 we 
summarize the measured values of the calibration parameters for the 50L detector. 



■^The residual range is the actual range minus the range already covered at the time of the deposition of a given 
dQ/dx. 
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Figure 5.13: dQ/dx as a function of the residual range for protons. Dots are direct measurements 
from the reconstructed hits of the proton tracks, while the points with errors bars are the estimated 
most probable of dQ/dx for each bin of the residual range. Data is fitted with a Bethe-Bloch 
corrected by the detector response function (Birk's law) in the hypothesis of pion (dashed line) 
and proton (full line) 
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5.3 Analysis of quasi-elastic i/^ charge current interactions 

5.3.1 Data taking and event selection 

The data set recorded with the 50L chamber amounts to 1.2 x 10"'^^ protons on target. As mentioned 
in Sec. 15.11 the trigger efficiency was monitored during data taking and its integrated value is 97%. 
Additional losses of statistics are due to the TPC (3%) and NOMAD (15%) dead time and to 
detector faults. These contributions add up to give an effective lifetime of 75%. Over the whole 
data taking period around 70000 triggers were collected in the TPC from which 20000 have at 
least a reconstructed muon possibly transversing the fiducial volume (see Sec. I5.3.3p . From all 
these a priory z/^ CC candidates, a visual scanning reveals that around half of them show a vertex 
in the fiducial volume, the rest are v interactions in the surrounding LAr bath (which also give 
trigger) plus a contamination of crossing muons due to veto inefficiencies. Therefore, we have a set 
of around 10000 z^^ CC events. 

The golden sample 

From the whole sample of ly^ CC events, we have selected a set of 86 events called the golden 
sample . They arc events having a clear two prong topology with a mip leaving the chamber in 
the beam direction and an identified contained proton. The main features of the golden sample 
are the following: 

• An identified proton of kinetic energy larger than 40 MeV, fully contained in the TPC and 
one muon whose direction extrapolated from NOMAD matches the outgoing track in the 
TPC (see Sec l533)) . 

• The distance of the primary interaction vertex to any of the TPC walls has to be greater 
than 1 cm. 

• The muon candidate track projected onto the wire plane must be longer than 12 wire pitches. 

• The presence of other stopping particles is allowed as far as their visible range does not 
exceed the range of a proton of 40 MeV kinetic energy. 

• No other tracks than the identified muon can leave the TPC. 

• Events with at least one converted photon with an energy deposition in the fiducial volume 
greater than 10 MeV are also rejected. 

The tightness of these selections makes the golden sample a very clean topology for visual scanning 
(Fig. Ol (Top)). 

The request of containment is very severe since the chamber volume is small. On the other hand, 
for a kinetic energy (Tp) below 40 MeV the proton range is comparable with the wire pitch and 
neither the proton momentum nor the interaction vertex can be reconstructed with due precision. 
However, for Tp > 40 MeV the 7r*/p misidentification probability is negligible (see Sec. I5.3.2|) . 

5.3.2 Proton reconstruction 

Proton identification and momentum measurement were performed using only the information 
provided by the TPC, following the reconstruction method explained in detail in Chap.Hl 
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CoHeciJon View 
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projections 



Induction View 
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Figure 5.14: Example of the 3D reconstruction of a low-multiplicity (goZden-like) v^^ CC event. 
(Top) The raw image from collection and induction wire planes: Hits and 2D track projections 
have been identified. (Bottom) Three dimensional view of the reconstructed event embedded in a 
3D recreation of the experimental setup described in Sec. 15.11 
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Proton identification 

In Sec. 14.6.41 we discussed that the average energy loss is a function of the mass and the kinetic 
energy of the particle. For the case of particles stopping in the fiducial volume of the chamber, a 
direct relation between their kinetic energy and their covered range (i.e. the distance a particle 
can penetrate before it loses all its energy) is obtained by means of the integration of Eq. (|4.19p 
along the particle's path (see Sec. I4.6.5|) . 

The discrimination between protons and charged pions is performed exploiting their different 
energy loss behavior as a function of the range. For the case of candidates stopping in the fiducial 
volume of the chamber, as the ones of the golden sample , p/tt* separation is unique (see Fig. 15.15)) : 
dQ/dx is measured along the proton candidate track and is compared with the different particle 
hypothesis. Monte-Carlo studies reveal that almost 100% of protons are identified as such, on the 
basis of their dQ/dx shape on the vicinity of the stopping point; in addition, the fraction of pions 
and kaons misidentified as protons is negligible, showing thus the superb particle identification 
capabilities of the Liquid Argon TPC technique. 




Figure 5.15: dQ/dx as a function of the residual range for the reconstructed proton in the event 
shown in Fig. 15.141 continuous (dashed) line represents the expected behavior of protons (pions) 
in the detector. 
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Momentum measurement 

Once we identify a contained proton, its kinetic energy is calculated from range, which only depends 
on the spatial reconstruction of the tracks (see Sec. 14.6^ . In this case, the momentum uncertainty 
is dominated by the finite pitch of the wires (2.54 mm and the performance of the reconstruction 
tools. For an uncertainty of ±1 wire in the determination of the starting and ending points of 
the projection of the track in both wire planes, and neglecting the error in the drift coordinate, a 
rough estimation of the error in the total length of the track is given by: 

L\/2 X 2.54 mm 
Ltot 

being L the projection of the track in the collection/induction plane and Ltot the total track 
length. The error in the range determination is automatically translated into an error in the 
energy measurement. This has been tested using a Monte-Carlo simulation of the detector: The 
resolution on the kinetic energy (Tp) varies from 3.3% for protons with Tp = 50 MeV up to 1% for 
Tp larger than 200 MeV. 

The angular resolution in the collection or induction view depends on the number N of wires 
hit by the particle along its path; it is 

CT ~ 0.36^12/(2.54 (5.10) 
which for a proton covering 10 wires {N = 10) means an error of 15 mrad. 



5.3.3 Muon reconstruction 

As mentioned in Sec. 15.11 the modest size of the 50L TPC made necessary the use of NOMAD as 
an external muon spectrometer in order to be able to measure the prong muons produced in CC 
interactions. The kinematic reconstruction of the outgoing muons exploits the tracking capability 
of the NOMAD drift chambers (see Sec. I5.1.2p . An event triggering the chamber will have at 
least one penetrating track reaching the Tl and T2 trigger scintillators bracketing the TRD of 
NOMAD (Sec. I5.1.3p . The corresponding track, nearly horizontal at the entrance of the NOMAD 
drift chamber volume, is reconstructed with an average momentum precision given by Eq. (|5.2|) in 
Sec. 15.1.21 A 10 GeV horizontal muon crossing all the chambers {L ^ h va) is reconstructed with 
a precision of 2.2%. 

Muon matching 

The reconstructed particle is traced back to the TPC accounting for the magnetic field and the 
presence of the forward NOMAD calorimeter (Sec. I5.1.2p . The latter introduces an additional - 
in fact, dominant - source of imcertainty due to multiple scattering (MS) in the iron plates of the 
forward calorimeter (see Sec. l5.1.2|) . For small scattering angles {9 rad), the MS uncertainty on 
the transverse momentum is independent of p and it turns out to be ~ 140 MeV. The correctness of 
the back-tracing procedure has been cross checked comparing the direction angles of the particles 
belonging to the mip sample, as measured by the TPC, with the corresponding quantity from 
NOMAql. 

*The equivalent pitch in the vertical direction (drift direction) is much smaller due to the high sampling rate of 
the fast ADC (360 /^m). 

^Due to the high sampling rate of the TPC, the angular resolution in the y — z plane, i.e. in the vertical plane 
along the nominal beam direction (^-axis), is dominated by the NOMAD uncertainty. 
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Figure 5.16: Example of the 50L-NOMAD matching for penetrating muons. 



In Fig. 15.161 we show an example of the measurement of a single muon: the traced-back muon 
predicted from the NOMAD measurement is shown in pink, while the trajectory of the same muon 
extrapolated from the TPC data (blue points) is drawn in green. Both measurements are compared 
to each other at the same z position located at the entrance of the NOMAD detector, in order 
to assure that the NOMAD muon is the same that traversed the TPC: Taking into account the 
errors in both measurements, we require that the difference between the traced-back muon and 
the original one is lower than 10 cm at NOMAD entrance point and their difference in angle below 
100 mrad. As we comment in Sec. 15.3.41 the Monte-Carlo simulation reproduces both 50L and 
NOMAD measurements and thus, it allows to test the efficiency of the latter cuts, showing that 
more than 97 % of the muons fuUfil those conditions (Fig. I5.17p . being the rest low-energy muons 
strongly affected by multiple scattering in the front calorimeter. 

The kinematic reconstruction of the outgoing muon from NOMAD allowed us to identify the 
events produced outside the fiducial volume that trigger the chamber: The vast majority are u 
interactions in the dewar or in the thermal bath plus a contamination of crossing muons due to 
veto inefficiencies. In Fig. l5.18l we show the distribution of the entrance point for all the triggering 
muons measured in NOMAD (in green) and just the ones which possibly transversed the active 
volume of the TPC (in red), which amount up to 30 % of the whole sample of triggering muons. 
As we commented at the beginning of Sec. 15.3. H around half of them were effectively originated 
in v interactions inside the chamber. 

5.3.4 The Monte-Carlo data sample 

Geometrical acceptance, background subtraction and QE selection efficiencies have been evaluated 
using a Monte-Carlo simulation. Neutrino- nucleus interactions are generated using NUX |104| . We 
use FLUKA |105j to simulate nuclear effects. NUX-FLUKA has successfully reproduced NOMAD 
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Figure 5.17: Difference between tfie SOL and NOMAD simulated measurements of muon direction 
as a function of the momentum. In real events, we reject events witli a difference higher than 
100 mrad. 
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Figure 5.18: The entrance point to NOMAD of all triggering muons (in green) and just the ones 
coming from the fiducial volume (in red). 
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data and therefore it is a reliable tool for the study of neutrino- nucleus interactions |106j . 100000 
CC QE and non-QE events have been generated simulating the WANE beam for the SOL 
chamber geometrical acceptance (see Fig. 15. 4p and final prompt particles emerging from neutrino- 
nucleus interactions are tracked through the 50 liter TPC and the NOMAD detectors by means of 
a GEANT4 simulation |107j . In Eig. 15.11 we show an output of a quasi-elastic Vfj, CC event: The 
leading proton is fully contained within the TPC, while the prong muon goes out the TPC, passing 
through the front calorimeter, entering in the magnet with the drift chambers and finally giving 
a trigger in the scintillators planes. In order to efficiently reproduce the geometrical acceptance, 
the geometry and materials of the different sub-detectors have been simulated with care, using the 
detailed technical information presented in Sec. 15.11 The energy deposition of all final particles is 
digitized emulating the TPC wire read-out, in order to reproduce detector resolutions and offline 
reconstruction efficiencies. 

Thanks to the simulation, full information about the final state particles and their energy 
deposition in LAr is provided. This information can be accessed both at particle level and at 
hit level, and besides, detailed information of secondary interactions is also available. Outside the 
TPC the simulation just take care of the trajectory of the muon, keeping the incoming position and 
momentum before the front calorimeter and immediately after, when the first of the drift chamber 
is hit. Einally, the impact position (if any) of the muon in the scintillator planes is recorded in 
order to decide if it gives trigger. The momentum measurement by NOMAD is also simulated 
and tracked-back to the LAr chamber taking into account multiple scattering effects in the front 
calorimeter (see Sec. I5.3.3|) . 

To exemplify the results of the simulation and reconstruction tools, we show in Eig. 15.191 the 
visible energy (left) and the transverse missed momentum (right) of the golden sample compared 
to the Monte-Carlo expectation. We see that the resolution on the energy measurement is rather 
accurate given the good precision measurement that both the TPC and NOMAD can achieve for 
the proton and the muon energy respectively. About the kinematical reconstruction in transverse 
plane, the missed momentum is utterly dominated by the muon reconstruction, since it is affected 
by the presence of the front calorimeter (Sec. I5.33|l . On the other hand, the momentum direction 
of fully contained protons are very accurately measured thanks to the excellent 3D reconstruction 
capabilities of the LAr TPC. 

5.3.5 Geometrical acceptance 

Due to the large amount of material between the TPC and the NOMAD trackers (more than 11 
interaction lengths), there is an important fraction of discarded events: The prong muons leaving 
the chamber after a CC interaction can be stopped or deflected by the front calorimeter giving 
as a result a non-triggered or a bad-matched event. Obviously, the latter affects essentially to low 
energy and large angle muons. Neither the energy loss nor the multiple scattering in the front 
calorimeter are enough to disturb more energetic muons. The latter effect is shown in Eig. 15.201 
(right) where we plot the geometrical acceptance as a function of the muon momentum for QE 
and DIS events separately. 

In pure QE interactions, the fraction of the incoming neutrino momentum carried by the prong 
muon is nearly equal to 1 due to the low energy transferred to the nucleon. Thus, the muon spectra 
slightly differs from the spectra of interacting neutrinos, which peaks around 20 GcV (Eig. 15.201 
(left)). On the contrary, non-QE (DIS + RES) interactions are characterized by a high momentum 
transfer with the nucleus, which gives as a result a muon spectra shifted to low momentum values 
(Eig. 15.201 (left)). The convolution of the muon momentum distributions with their respective 
geometrical acceptance leads to a very significant difference between the total acceptance for QE 
events (96%) and non-QE events (72%). 
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Figure 5.19: Visible energy (left) and transverse miss momentum (right) for QE golden events. The 
squared markers correspond to data. The continuous (dashed) line represents the reconstructed 
(generated) events for Monte-Carlo neutrino QE interactions normalized to the golden sample data. 
The dotted line corresponds to the same quantities computed by adding the reconstructed proton 
momentum to the generated muon momentum. 
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Figure 5.20: (Left) Simulated distributions in muon momentum for QE (solid line) and non-QE 
events (dotted line). These distributions are normalized to the same (and arbitrary) number of 
events. (Right) Geometrical acceptance as a function of the muon momentum for QE (solid line) 
and non-QE events (dotted line) 
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5.3.6 Background estimation 

The golden sample contains pure QE interactions (j^^ n — > fjr p) plus an intrinsic background 
dominated by resonant production followed by pion absorption in the nucleus {Vf^ p — > A++/x~ —f 
fi~ p TT"*" ,iyfj^n—> A+/i^ — > /i~ p TT^). There is also an instrumental background due to resonance 
(RES) and deep-inelastic (DIS) events where final state neutral particles (tt^'s, n's and 7's) scape 
undetected. On the other hand, the n n~ — > fi~ n Tr~^ contamination is negligible in this 

tightly selected sample due to the superb tt^ /p identification capabilities of the LAr TPC. 

Thanks to the Monte-Carlo information we can properly simulate the selection cuts imposed to 
the golden events (see Sec. I5.3.ip over the full sample of generated events. Imposing such selection 
criteria, an algorithm efficiently selects the golden events from the whole sample of simulated QE 
interactions. We found that 16.6 % of them are golden, the rest are events with more particles in 
the final state due to nuclear re-interactions or having a not contained proton^ 

On the other hand, we have analyzed the sample of CC deep-inelastic (DIS) and resonant 
(RES) events in order to evaluate inefficiencies of the vetoing selections for n ^ A+/i^ 
li~ p 7r° (i.e. the probability to miss both the decay photons of the 7r° 77 or the e+e~7 system 
in case of tt" Dalitz decay) and for events with several neutral particles (neutrons and energetic 
photons from nuclear interactions) which scape from the detector (see Fig. I5.2ip . The search for 
this irreducible background in the DIS plus RES samples reveals that 0.14% are golden-\ike events. 




Collection wires. \i28 wiras: 32 em.) Induction wires, (ise wir^s; sscm.) 



Figure 5.21: The raw image of a low multiplicity real event in the collection (left) and induction 
plane (right). The event is reconstructed as {v^ n — > /x~A+ /i^ p tt*^) with a mip leaving the 
chamber, an identified stopping proton and a pair of converted photons from the 7r° decay. When 
these photons escape from the chamber, the event is tagged as a golden event. 

Taking into account the relative weight between QE and non-QE eventtQ, we find that the total 
expected contamination of the golden sample is 28% (around 40% of them correspond to events 
with an escaping tt*^). In Table [531 we summarize the parameters used in background estimation. 



"The automatic selection eflicieney was cross-checked by means of visual scanning in a smaller sample of simulated 
data. 

^Respect to the total number of i/^ CC events, in our simulations, 2.3% are QE, 91.5% are DIS and the remaining 
6.2% corresponds to RES events. The simulation of the latest is not done treating each resonance individually, an 
average of them is taken into account instead. 
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The Monte-Carlo prediction for escaping 7r° can be checked using test samples that consist of 
golden events with one (A'^i) or two (-/V2) converted photons pointing to the interaction vertex. 
Assuming the gamma identification probabihty to be uncorrelated for the two photons, we have 
A^i = 2A^(1 — e-y)e-y and N2 = Ne?^, N being the (unknown) overall rate of p n~ 7r° final states and 
the photon identification efficiency. After the scanning of the test samples, = (7Vi/2iV2 + 1)""'^ 
turned out to be (43± 9) %. Hence, the probability of missing both gammas is (32±10) % to be 
compared with the MC calculation of 20.4%. 

The estimated background is of utmost importance in the analyses of QE interactions being 
presented. It is key for a proper estimation the QE cross section calculation (Sec. l5.4p . Furthermore, 
as we will discuss in Sec. I5.3.8| adding this sizable background leads to a much better description 
of the event kinematics. 





QE 


non-QE 


Relative weight 
Geom. acceptance 
Golden fraction 


2.3 % 
96 % 
16.6 % 


97.7 % 
72 % 
0.14 % 



Table 5.3: Summary of parameters for background estimation. 



5.3.7 Event rates 

We use Monte-Carlo to estimate the expected number of events in the detector. The simulation 
of the beam predicts a flux of 2.37 x lO"'^ I'^/cm^/p.o.t. over the TPC exposing area. This flux 
convoluted with the neutrino cross sections and scaled to the flducial mass of the detector (65.3 Kg 
of LAr) gives an event rate of 2.05 x 10~^^ f^CC/p.o.t. Now assuming a total exposure of 1.2 x 10^^ 
p.o.t. and an effective hfetime of 75% (see Sec. l5.3.1"|) . the total number of CC is equal to 18450. 
Out of these 18450 events, only muons above a certain momentum threshold and within angular 
acceptance will trigger the detector. The latter parameters for the experimental configuration are 
summarized in Table [SHI Taking the real data muons that triggered the system, we saw that most 
of them are above 8 GeV in momentum and below 300 mrad in angle. If we apply those cuts to 
the Monte-Carlo samples, we end up with 18450 x 0.023 (QE fraction) x 0.96 (efficiency of the 
acceptance cuts for QE), equal to 400 quasi-elastic events. For DIS-I-RES events: 18450 x 0.977 
(fraction of DIS-I-RES) x 0.72 (efficiency of acceptance cuts) gives 12600. In total we expect 13000 
events to be compared with the 10000 CC we have after visual scanning of good triggers (see 
Sec lOljl . 

The expected number of golden events is obtained taking the 16.6% (golden fraction) of the 
total QE, and adding the corresponding 28% due to background contribution (see Sec. 15. 3. 6| . which 
finally gives 85 ± 9(stot) ± l?>[sys) golden events to be compared with the 86 we observe. The 
systematic uncertainty is dominated by the fraction of QE events (2.3%) and the beam simulation 

(8%) El. 

5.3.8 Analysis of the kinematics of quasi-elastic interactions 

In spite of the limited statistics (86 events), the golden sample provides information on the ba- 
sic mechanisms that modify the kinematics of neutrino-nucleus interactions with respect to the 



detailed discussion about the different sources of systematic errors can be found in Sec. 15.41 where the CC 
QE cross section is calculated. 
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Neutrino Flux 


2.37 X IQ-^ 


I'^/cTO'^/p.O.t. 


Interaction rate 


2.05 X 10"^^ 


v^j, CC/p.o.t. 


Total exposure 


1.2 X 10^9 


p.o.t. 


Effective detector lifetime 


75 


% 



Table 5.4: Summary of the experimental configuration parameters for event counting. 



corresponding neutrino-nucleon process. Nuclear matter perturbs the initial state of the interac- 
tion through Fermi motion; it also affects the formation of the asymptotic states through nuclear 
evaporation, hadronic rc-scattcring or hadronic re-absorption. Several kinematic variables arc only 
marginally affected by nuclear effects; in this case, the corresponding distributions can be repro- 
duced once the ^/-nucleon interaction is corrected for Fermi motion and Pauli blocking. Clearly, 
purely leptonic variables belong to this category. There are also a few hadronic variables whose 
distribution is strongly influenced by the selection cuts but show limited sensitivity to nuclear 
effects. In particular, the proton kinetic energy is bounded by the Tp > 40 MeV cut and the 
requirement of full containment in the fiducial volume. This distribution is shown in Fig. 15.221 to- 
gethcr with the transverse momentum of the proton for the golden sample. Similarly, in Fig. 15.231 
we show the kinetic energy and the transverse momentum distributions of the muon. Figs. 15.221 
and 15.23) can be used as a consistency check. They demonstrate that Monte-Carlo reproduces the 
kinematic selection performed during the scanning and analysis of the golden channel. They also 
show that our selected sample contains a non-negligible contamination from non-QE events. As 
already indicated, we estimate this contamination to be 28 %. 




Figure 5.22: Distribution of the proton kinetic energy and transverse momentum for the golden 
sample (squared marks). Dotted histogram represents the expectation from simulated Monte-Carlo 
QE events, while the hatched one represents the simulated background from non-QE events. Both 
contributions arc summed in the filled histogram and normalized to the golden sample data. 
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Figure 5.23: Distribution of the muon kinetic energy and the muon transverse momentum for the 
golden sample . Dotted histogram represents the expectation from simulated Monte-Carlo QE 
events, while the hatched one represents the simulated background from non-QE events. Both 
contributions are summed in the filled histogram and normalized to the golden sample data. 



Other variables embedding the reconstructed kinematics of the protons are sensitive to genuine 
nuclear effects. In particular, we analyzed two of them: the acollinearity and the missing transverse 
momentum of the event. The former is defined as 



A 



Pxp Pxfj, ~r Pyp Pyf-i 



{^Pxp ^ Pyp^iPxfj 



(5.11) 



Pxp and Pyp being the transverse momentum components of the proton and Pxp, and py^ the 
corresponding quantities for the muon. For a purely QE scattering on a nucleon, the muon and 
the proton ought to be back-to back in the transverse plane so that the acollinearity is zero. 
Fig. 15.241 (right) shows the acollinearity distribution for the golden sample and the expectation 
from simulations. For this particular variable, the inclusion of nuclear effects in the Monte-Carlo 
does not show a striking difference with respect to the case where no nuclear effects are taken into 
account. The selection cuts (in particular, the fact that the proton should be fully contained), 
the resolution NOMAD has for muon reconstruction and the contamination from non-QE events 
are the main reasons for the appearance of events with large acollinearity values. The distortion 
introduced in the tail of the acollinearity distribution by nuclear effects is much smaller than the 
one due to detector effects. Hence this variable is not adequate to show how the event kinematics 
varies in the presence of nuclear matter. 

The influence of nuclear effects on the event kinematics is best seen when we consider the missing 
transverse momentum P^^^^ (see left plot on Fig. I5.24p . We observe that a naive approach that 
takes into account Fermi motion and Pauli blocking disregarding nuclear effects (dashed histograms 
in Fig. 15.241 left) does not reproduce the P^j^^ data. Once nuclear effects are added, simulated 
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events show a good agreement with data. The effect of background events is also relevant. Being 
the non-QE events more unbalanced than QE events in the transverse plane, they have the effect 
of shifting the expected theoretical distributions to higher values of P^igg, just like data do. 




Figure 5.24: Distribution of the missing transverse momentum (left) and acoUinearity (right) for 
the golden sample . Dotted histogram represents the expectation from simulated Monte-Carlo QE 
events, while the hatched one represents the simulated background from non-QE events. Both 
contributions are summed in the filled histogram and normalized to the golden sample data. The 
dashed histograms show the expected distributions in case no nuclear effects arc taken into account. 



In order to quantify the level of agreement between the theoretical expectations and the data, 
we have performed a Kolmogorov test [TSJ 1108] for some of the transverse variables. Such a test 
is based on a comparison between the cumulative distributions of the samples being analyzed (see 
Fig. 15.251 where cumulative distributions for P^^^g and the acoUinearity are shown). Its goal is 
to give a statistical estimation about how compatible is the data sample with being a random 
sampling from a given theoretical distribution. In Tab. 15.51 we quote the results of such a test in 
terms of the Kolmogorov probability. There, we test experimental data against three theoretical 
approaches: golden QE with no nuclear effects, golden QE with nuclear effects included and, finally, 
against golden QE plus the corresponding DIS and RES backgroun Clearly, nuclear effects offer 
a much better description of data than naive approaches commented above, mainly for the case 
of P^i^iss (see Tab. 15. 5p . On the other hand, the effect of adding the background leads to a better 
modeling of the missing transverse momentum and the acoUinearity. 

In summary, we have seen that nuclear effects are an important source of perturbation for the 
kinematics of quasi-elastic neutrino interactions, however the smallness of the accumulated statis- 
tics does not allow a systematic survey of these effects |109( IllOj . Still, these results empirically 
demonstrate the effectiveness of LAr imaging in the characterization of low-multiplicity multi-GeV 
neutrino interactions [2]. In particular, as we discuss in Chap. a good description of trans- 



Theoretical distributions are always normalized to the golden sample data. 
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Variable 


Kolmogorov probability 


QE w/o nucl. off. 


QE with nucl. cff 


QE+DIS+RES with nucl. cff. 




0.35 


0.99 


0.98 


miss 


0.01 


0.28 


0.99 


Acoll 


0.38 


0.61 


0.99 



Table 5.5: Table of Kolmogorov probabilities for transverse variables (see text). 




Figure 5.25: Cumulative distribution of the missing transverse momentum and the acollinearity for 
the golden sample. Continuous line represents the expectation from simulated Monte-Carlo events 
(QE+DIS+RES). 



verse variables is fundamental since they are key to a kinematic-based approach for v^^ Vrive) 
oscillation appearance [1]. 
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5.4 Measurement of the quasi-elastic CC cross section 

The collected sample of golden events gives us the opportunity of measuring for first time the 
quasi-elastic CC cross section using a LAr TPC. To this purpose a precise knowledge of the 
beam composition, the beam profile, the exposition time, trigger efficiencies, acceptance cuts, event 
selection rates and background ratios is needed. In previous sections most of the required quantities 
have been already obtained, specially in Sec. I5.3.7l whcrc an estimation of the expected number of 
QE events was discussed. In fact, the calculation of the cross section is like inverting the process 
we followed there, i.e. given the total number of measured golden events extract the quasi-elastic 
cross section. However, in order to give a competitive measurement instead of just a cross-check, 
a careful treatment of systematic errors has been done. 



Event counting 

As we wrote in Sec. 15.3. 1[ the total number of golden events collected is = 86 events. The data 
taking period corresponds to a total exposure E^p = 1.2 x 10^^ protons on target, which must 
be corrected by the effective lifetime of the detectors Ti^g = 75% (see Sec. I5.3.7p . The events 
contained in the golden sample are pure QE interactions {v^ n fi^ p) plus an intrinsic and 
instrumental background (see Sec. 15.3.^ . To compute the QE cross section, this background must 
be clearly understood and removed. 



Background subtraction 

The golden sample is contaminated with a sizable fraction of deep- inelastic (DIS) and resonance 
(RES) events. To estimate this contamination, we use 

jynQE j:^nQE Acc"'^^ Q^^QE 

^ ^ 1^ = irW AccQ^ 'goE (5-^2) 

]^QE ^j\^nQE-^ ^j^g number of true QE (non-QE) events inside the golden sample {N = N^^ + 
j^nQEy ^ stands for the fraction of a given event class (i.e., in our Monte-Carlo, QE events 
amount to R^^ = 2.3% of the total. For non-QE events ^""3^ = 97.7%). Acc^^ {Acc""^^) is 
the acceptance for QE (non-QE) events and according to our simulations their average values are 
= 0.96 and Acc''^'^ = 0.72. Finally, Q measures the fraction of events of a given class 
that are classified as golden. For QE events, we obtain = 0.166. The fraction of non-QE 
events classified as golden is C/"^^ = 0.00145. Despite the smallness of this number, the final 
contamination in our golden sample is not negligible since, at these high energies, the majority of 
collected events are of non-QE nature. In particular, given the reduced size of our detector, most 
of the background comes from undetected 7r° since, as computed in Sec. 15.3.61 the probability to 
miss the two photons is almost one third. We have estimated that the final contamination inside 
the golden sample due to non-QE events is B = 0.28. 



The measurement 

After background subtraction, the golden sample contains A^"^^ = = 67±8 (stat) genuine QE 
events. The beam simulation predicts a neutrino flux of $ = 2.37 x 10~^ z^^/cm^/p.o.t. Knowing 
the values for the effective detector lifetime and the total exposure, we compute the total cross 
sections using 

'"'^'^ " $ X X Nf,d X Tu., X AccQ^ x QQ^ ^^'^^^ 
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where Nfid = 2.167 x 10^ is the number of target neutrons contained in the fiducial volume. The 
figures used for the rest of the variables have been discussed above. Hence, for an average neutrino 
energy of 29.5 GeV, the total QE cross section amounts to 

GQE = (0.90 ± 0.10 (stat)) X lO^^^^cm^ 

Although a wide-band beam has been used, the reduced statistics accumulated during the test, 
do not allow to obtain a set of total cross sections measurements for several incoming neutrino 
energies. 



Evaluation of systematics 

The systematic error associated to our measurement was evaluated exploiting the full Monte-Carlo 
simulation of the experimental set-up (Sec. 15. 3~^ . To better identify the different contributions to 
the systematic uncertainty, it is much more convenient to express Eq. (|5.13p , after a trivial change 
of variables, as: 

_ N_ 1 

'^'^^ " $ X E^p X Nf,d X TH„e ^ AccQE X QQi^ + (R^Q^/RQE) x Acc^Q^ x G^Q^^ 

We identify three types of systematic errors: The first one involves the knowledge of the exper- 
imental setup, the neutrino flux and the detector performance (see Sec. 15. ip . The second is related 
to the golden event selection criteria (Sec. I5.3.T|) and the geometrical acceptance (Sec. 15.3.5]) . Fi- 
nally, we have the systematics associated with the determination of the fraction of golden events 
and the expected ratio between QE and non-QE events, which strongly depends on the modeling 
of neutrino-nucleus interactions (see Sec. I5.3.4p and the experimental uncertainties in the cross 
sections at these energies. 

In the following, we summarize and describe in detail the different sources of systematics 
entering in the cross section calculation through Eq. 1)5.14^ : 



The uncertainty in the WANE neutrino flux has been taken from [99j and it amounts to 
8%. The systematic error on the detector lifetime is 4%. It is dominated by the systematics 
associated to the NOMAD detector lifetime |92) . 

The uncertainty associated to the fiducial volume has been evaluated varying its dimen- 
sions by ±1 wire in the wire coordinate and ±10 drift samples on the time coordinate. The 
measured resolution on the time coordinate is much better, however we have adopted a con- 
servative approach when evaluating this kind of systematic error. The discussed modification 
of the fiducial volume causes a variation on the number of golden events that amounts to 
3%. We take this value as the systematic error associated to the fiducial cut. 

The criteria used for event selection are another source of systematics. The energy resolution 
for stopping hadrons is estimated to be about 4% and for converted photons we take 11%/\/E. 
This translates into a final systematic error of 6% due to the definition of a golden candidate. 

The main source of systematics for the geometrical acceptance comes from the criteria used 
to match the reconstructed muons in NOMAD with the ones that exit the 50 liter LAr TPC 
(see Sec. I5.3.3| . Due to the small lever-arm of the TPC, the resolution on the measurement 
of the muon direction is, on average, 10 mrad. To evaluate how this uncertainty affects the 
muon acceptance, we have modified the matching criteria described above by ±10 mrad. The 
error on the geometrical acceptance turns out to be 2% for QE and non-QE samples. These 
errors translate into a systematic contribution of 1.5% (0.5%) for QE (non-QE) events. 
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• The fraction of golden events inside the QE and non-QE samples is affected by nuclear effects, 
which modify the topology and kinematics of the final state. Following the discussion in |lllj . 
where the acceptance for identification of QE CC events in a fine grained detector is studied 
as a function of the proton detection threshold and nuclear effects, we conservatively take 
that varies by ±25% due to uncertainties in the modeling of nuclear effects inside the 
Monte-Carlo. With this assumption, the systematic error on the cross section due to the 
fraction of golden events inside the QE sample amounts to 16%. This is the dominant error 
source. A similar error on Q^Q^ (±25%) contributes to 5% on the total systematic. 

• Finally, the ratio of QE versus non-QE events is assumed to be known at the level of 20% 
from experimental uncertanties in the cross section |112| . Its contribution to the cross section 
systematic error is 4%. 

Table summarizes the contributions from the different systematic uncertainties. 

All in all, the estimated total systematic error is 20.6%. This translates onto a measurement 
of the total QE i/^ charged-current cross section: 

aQE = (0.90 ± 0.10 (stat) ± 0.18 (sys)) x lO^^^cm^ . (5.15) 

This value agrees with previous measurements of this reaction reported in the literature (see 
Fig. EM and references [HI [112 El [HTl [IH \m\), and it stands as the first Physics 
measurement ever done with a LAr TPC using a neutrino beam from an accelerator. 



Source Systematic error (%) 

Neutrino flux 8 

Detector lifetime 4 

Vertex detection (fiducial cut) 3 

Golden event selection cuts 6 

Acceptance for QE events, Acc'^^{±2%) 1.5 

Acceptance for non-QE events, Acc^^^{±2%) 0.5 

Fraction of golden events in QE sample, C/*3^ (±25%) 16 

Fraction of golden events in non-QE sample, Q"^^^ (±25%) 5 

Fraction of non-QE to QE events, -^ff (±20%) 4 

Total 20.6 



Table 5.6: Summary of systematic errors. The largest source of error comes from nuclear effects 
affecting the expected fraction of golden events inside the genuine QE sample. 
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Figure 5.26: Muon neutrino- nucleus quasi-elastic cross section data. Continuous line is the theo- 
retical behavior of the QE cross section in free nucleon approximation. The measurement reported 
in this work appears in the figure as a star. The shown vertical error bar is obtained adding in 
quadrature the systematic and statistical uncertainties. The horizontal error bar corresponds to 
the RMS of the simulated neutrino energy spectrum (it amounts to 18 GeV). 



Chapter 6 



Searching for z/^ oscillations 
at CNGS 



The Super-Kamiokandc Collaboration has measured a first evidence of the sinusoidal behaviour 
of neutrino disappearance as dictated by neutrino oscillations [120] . However, although the most 
favoured hypothesis for the observed z/^ disappearance is that of v^J^ — > oscillations, no direct 
evidence for v-^ appearance exists up to date. A long baseline neutrino beam, optimized for the 
parameters favoured by atmospheric oscillations, has been approved and is now running in Europe 
to look for explicit appearance: the CERN-Laboratori Nazionali del Gran Sasso (CNGS) beam 
|121j . The approved experimental program consists of two experiments ICARUS [55] and OPERA 
[122] that will search for Vt oscillations using complementary techniques. 

Given the previous experimental efforts [HH [M] and present interest in direct appearance, 
we assess in this Chapter the performance of several statistical techniques applied to the search 
for v-r using kinematic techniques. Classic statistical methods (like multi-dimensional likelihood 
and Fischer's discriminant schemes) and Neural Networks based ones (like multi-layer perceptron 
and self-organized neural networks) have been applied in order to find the approach that offers the 
best sensitivity. 

The underlying idea consists on finding the kinematical features which characterize the signal 
respect to the dominant background. To this purpose is obvious that optimal reconstruction 
capabilities are of utmost importance. In Chapter [5] we have assessed the excellent performance 
of our reconstruction tools when they are applied to the analysis of the kinematics of i^^ CC 
events. Furthermore, we have also seen that our Monte-Carlo accurately describes neutrino-nucleus 
interactions and nuclear effects. Thanks to these successfully tested capabilities, in this Chapter, 
we can use them to evaluate the performance of a LAr TPC when searching for direct appearance 
of Ur at CNGS. 

6.1 On the problem of appearance at CNGS 

In the CERN-NGS beam [12111123] . the expected i^e and i^r contamination are of the order of 10^^ 
and 10^^ respectively compared with the dominant t'^ composition. These properties allow to 
search for oscillations by looking at the appearance of Ve and charged current events. In these 
cases, the detector must be able to tag efficiently the interaction of i/^s and z/i-'s out of the bulk 
of i>p, events. 



Ill 
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In the case of — > oscillations, the golden channel to look for the Vt appearance is the 
decay of the tau into an electron and a pair neutrino anti-neutrino due to: (a) the excellent electron 
identification capabilities of a LAr TPC; (b) the low background level, since the intrinsic Ve and 
Df. charged current contamination of the beam is at the level of one per cent. 

Kinematical identification of the r decay, which follows the z/,- CC interaction, requires excellent 
detector performance: good calorimetric features together with tracking and topology reconstruc- 
tion capabilities. In order to separate Vt events from the background, two basic criteria can be 
used: 

• An unbalanced total transverse momentum due to neutrinos produced in the r decay. 

• A kinematical isolation of hadronic prongs and missing momentum in the transverse plane. 

6.1.1 Detector configuration and data simulation 

As a reference for this study, we consider a detector configuration consisting of 5 identical T600 
modules (T3000 detector) 0. As the active Argon for one of these modules is 470 tons, the total 
mass amounts to 2.35 ktons. We also assume five years running of the CNGS beam in shared mode 
(4.5 X 10^^ p.o.t. per year), which translates in a total exposure of 5 x 0.47 x 5 = 11.75 Ktonxyear. 

Events have been fully simulated in the whole of each T600 module (« 3.2 ft. x x 18/ m^), 
but only 65% have the primary interaction vertex inside the active part of the LAr volume. 

The study of the capabilities to reconstruct high-energy neutrino events was done using fully 
simulated j/gCC events inside the considered LAr detector |124| . The neutrino cross sections and 
the generation of neutrino interactions is based on the NUX code |104j : final state particles are 
then tracked using the FLUKA package |105j . The angular and energy resolutions used in the 
simulation of final state electrons and individual hadrons are identical to those quoted in |69j : 
a{E)/E equal to 20%/ ® 5% for hadrons and 3%/a/E © 1% for electromagnetic showers {E 
in units of GeV). Concerning the reconstruction of the particle direction, the assumed angular 
resolutions are 0.13//E for electrons and 0.04 for single hadrons. In the analyses, the particle four 
vectors are smeared according to the previous parameters. 

The fiducial volume 

The ability to look for tau appearance events is limited by the containment of high energy 
neutrino events. Energy leakage outside the active image volume creates tails in the kinematical 
variables that fake the presence of neutrinos in the final state. We therefore impose fiducial cuts 
in order to guarantee that on average the events will be sufficiently contained. 

The fiducial volume is defined by looking at the profiles of the total missing transverse momen- 
tum and of the total visible energy of the events. The average value of these variables is a good 
estimator of how much energy is leaking on average. 

The profiles of missing Pt and visible energy as a function of the vertex positions for the T600 
module are shown in Fig. 16.11 We have to leave 20 cm at the top and bottom of the active volume 
(x coordinate), and 25 cm from the sides (y coordinate). Longitudinally, we cut out 1.7 m from the 
total length (z coordinate). This results in a total fiducial volume of 2.7 /i xhw x 16.3/ m^. Thus 
with this conservative approach, we keep 65% of the total number of events occurring in the active 
LAr volume. The total exposure in a 5 T600 setup with 5 years shared running CNGS, taking in 



^See Sec. 12.21 for a brief description of the ICARUS detector or consult Ref. |69| for a much more detailed 
description. 
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account the fiducial cuts, is 7.6 ktonxyear. In what follows, tau selection and background rejection 
eSiciencies are normalized to this figure. 
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Figure 6.1: Performance of the energy flow reconstruction in a T600 module: Transverse momen- 
tum versus x (top plot) and y (middle plot) position. Average visible energy versus z position 
(bottom plot) (a cut E^, < 20 GeV has been imposed). 



6.1.2 Expected rates and generated data 

Tab. [O] summarizes the total event rates expected for 5 T600 modules and 5 years CNGS running. 
It also shows, for different values of Am^, the expected Vr CC rates in case oscillations take place. 

The rate has a strong dependence on the exact value of the Am^ in the parameter region 
favoured by atmospheric data as we see in table [HH] For our central value taken from atmospheric 
neutrino results (Am^ = 3 x 10"^ eV'^), the number of t/^-CC with r e is about 50; the signal 
over background ratio is 50/252 ~ 0.2. 

Tab. 16.21 summarizes the total amount of simulated data. We note that Vt (v,,) CC sample, 
generated in active LAr, is more than a factor 250 (50) larger than the expected number of collected 
events after five years of CNGS running. 

To save CPU power, only i^eCC events with an incoming neutrino energy below 30 GeV have 
been generated, since oscillated events, with a Arri^ = 3 x 10~^eV^, tend to accumulate in the low 
energy region (only 3% of the total VrCC events have an energy above 30 GeV). 
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v^, CC 
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252 
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Ve CC 


17 


11 
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10600 


6890 


V NC 


245 


159 


Vr CC, Am^ (cV^) 




1 X 10-^ 


31 (6) 


20 (4) 


2 X 10~3 


125 (23) 


81 (15) 


3 X 10-3 


280 (50) 


182 (33) 


5 X 10^3 


750 (135) 


488 (88) 



Table 6.1: Expected event rates in active and fiducial LAr for the 5 years of CNGS running and 
5 T600 modules. For standard processes, no oscillations are assumed. For j/^CC, we take two 
neutrino t'p — > Vt oscillations with sin^ 29 = \. In brackets, the number of r — > e events (t decay 
rate is 18%). 



Process 




VrCC 


Total LAr 


22500 


22100 


Active LAr 


14200 


13900 


Fiducial Vol. 


9250 


9000 



Table 6.2: Amount of fully generated data. 

6.2 Statistical pattern recognition applied to oscillation searches 

In the case of a v^^ — > Vt oscillation search with Liquid Argon, the golden channel to look for v^- 
appearance is the decay of the tau into an electron and a pair neutrino anti-neutrino due to: (a) 
the excellent electron identification capabilities; (b) the low background level, since the intrinsic 
Vf, and Vf, charged current contamination of the beam is at the level of one per cent. 

Kinematic identification of the t decay |125] . which follows the i^rCC interaction, requires 
excellent detector performance: good calorimetric features together with tracking and topology 
reconstruction capabilities. In order to separate Vr events from the background, a basic criteria 
can be used: an unbalanced total transverse momentum due to neutrinos produced in the t decay. 

In Fig. 16. 21 we illustrate the difference on kinematics for signal and background events. We plot 
four of the most discriminating variables: 

• Ems'- Visible energy. 

• P™'"^ : Missing momentum in the transverse plane with respect to the direction of the incident 
neutrino beam. 

• p!^^: Transverse momentum of the prompt electron candidate. 
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Figure 6.2: Visible energy (top left), transverse missed momentum (bottom left), transverse elec- 
tron momentum (top right) and pi (bottom right). Histograms have an arbitrary normalization. 



Throughout this work, we take into account only the background due to electron neutrino 
charged current interactions. Due to the low content the beam has on Vf., charged currents in- 
teractions of this type have been observed to give a negligible contribution to the total expected 
background. We are confident that neutral current background can be reduced to a negligible level 
using LAr imaging capabilities and algorithms based on the different energy deposition showed by 
electrons and tt" (see for example |126j ). Therefore it will not be further considered. The contam- 
ination due to charm production and CC events, where the prompt muon is not identified as 
such, was studied by the ICARUS Collaboration [127] and showed to be less important than j/gCC 
background. 
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6.2.1 Direct Background rejection using Sequential Cuts 

The most straightforward way to look for Vr appearance is by means of a series of sequential cuts. 
This procedure does not take into account possible correlations among the selected kinematical 
variables. We take advantage of the soft energy spectrum for oscillated induced events with 
respect to the visible energy spectrum of v^GC Oscillated events tend to accumulate at relatively 
low values of p!^^ , as well, therefore an upper cut on p!^^ efficiently removes the background. Since 
I'eCC events are well balanced and very few populate the high P™*"''* region (missing momentum in 
the plane perpendicular to the incoming neutrino direction), a cut in this variable is also applied. 

Tab. 16.31 summarizes the list of sequential cuts applied to reduce the i^gCC background. The 
expected number of signal events, for Am^ = 3 x lO^^ eV^, after all cuts is 11.7 ± 0.4, for a total 
background of 2.2 ± 0.3 expected eventfU. Fig. 16.31 shows some kinematical variables for signal and 
background events surviving the cuts. The signal to background ratio (S/B) for this approach is 
around five. More sophisticated techniques, taking into account correlations among discriminating 
variables, should be used in order to improve the S/B value. 





l^rCC{T e) 




i^^CC (r -> e) 

Ato^ = 
3 X 10-3 eV^ 


Cuts 


Efficiency 

(%) 


lyeCC 


Initial 


100 


252 


50 


Fiducial volume 


64 


161 


32 


Evisible < 18 GcV 


62 


37.2± 1.1 


30.8 ±0.7 


P^^ < 0.9 GeV 


51 


20.0 ±0.8 


25.6 ±0.6 


Q'^P > 0.3 GeV 


43 


18.5 ±0.8 


21.6 ±0.6 


P^'"^ > 0.6 GeV 


23 


2.2 ±0.3 


11.7±0.4 



Table 6.3: Rejection of the t^gCC background in the r ^ e analysis. 5 years of CNGS nmning and 
5 T600 modules exposure has been considered. We quote statistical errors only. 



6.2.2 Kinematical search using classic statistical methods 
A multi-dimensional likelihood 

The second method adopted for the t ^ e analysis is the construction of a multi-dimensional 
likelihood function (see for example [128j ). which is used as the unique discriminant between 
signal and background. This approach is, a priori, an optimal discrimination tool since it takes 
into account correlations between the chosen variables. 

A complete likelihood function should contain five variables (three providing information of the 
plane normal to the incident neutrino direction and two more providing longitudinal information). 
However, in a first approximation, we limit ourselves to the discrimination information provided 
by the three following variables: Emsibie, P™''"* and pi. 

As we will see later, all the discrimination power is contained in these variables, therefore we 
can largely reduce the complexity of the problem without affecting the sensitivity of the search. 
Two likelihood functions were built, one for r signal (£5) and another for background events [Cb]- 
The discrimination was obtained by taking the ratio of the two likelihoods: 



^Errors in the number of expected events are of statistical nature. The usual Poisson statistical error in probability 
density estimation is rescaled according to the distribution normalization (5 years of CNGS running). 
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T3000 detector. 5 years CNGS running 
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Figure 6.3: Sequential cuts. Distributions of VrCC events (solid histograms) and i^eCC events 
(hatched histograms) after sequential cuts are imposed. 



ln{X) = C{[E„s,Ue,PT''',Pl]) = 
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(6.1) 



In order to avoid a bias in our estimation, half of the generated data was used to build the 
likelihood functions and the other half was used to evaluate overall efficiencies. Full details about 
the multi-dimensional likelihood algorithm can be found elsewhere |126j . However, we want to 
point out here some important features of the method. 

A partition of the hyperspace of input variables is required: The multi-dimensional likelihood 
will be, in principle, defined over a lattice of bins. The number of bins to be filled when constructing 
likelihood tables grows like n'^ where n is the number of bins per variable and d the number of these 
variables. This leads to a "dimensionality" problem when we increment the number of variables, 
since the amount of data required to have a well defined value for InA in each bin of the lattice will 
grow exponentially. 

In order to avoid regions populated with very few events, input variables must be redefined to 
have the signal uniformly distributed in the whole input hyperspace, hence EmsiUe, -P™*^^ and pi 
are replaced by "flat" variables (see Fig. 16. 4p . In general, a variable g{a) defined as 



g{a) = f{x{a)) 



dx 
da 



(6.2) 
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is uniformly distributed ("flat") between and 1 by choosing a{x) = f{x')dx', where a; is a 
variable distributed according to a p.d.f. f{x). With this approach, more bins are automatically 
taken in those regions of faster variation of the signal. The combination of the "flat" variables into 
a single one did not give rise to an uniform three dimensional distribution, due to the presence of 
correlations between the variables. 

Besides, an adequate smoothing algorithm is needed in order to alleviate fluctuations in the 
distributions in the hyperspacc and also, to provide a continuous map from the input variables to 
the multi-dimensional likelihood one (InA). 




Figure 6.4: Comparison for "flat" Euisibie, P-r^ , Pi and P™'''* variables between t signal and 
CC events. Arbitrary normalization has been taken into account when plotting background events. 

Ten bins per variable were used, giving rise to a total of 10"^ bins. Fig. 16.51 shows the likelihood 
distributions for background and tau events assuming five years running of CNGS (total exposure 
of 7.6 kton X year for events occurring inside the fiducial volume). 

Tab. 16.41 shows, for different cuts of InA, the expected number of tau and i/g CC background 
events. As reference for future comparisons, we focus our attention in the cut InA > 1.8. It gives a 
signal selection efficiency around 25% (normalized to the total number of r events in active LAr). 
This r efficiency corresponds to 13 signal events. For this cut, we expect 1.1 ± 0.2 background 
events. After cuts are imposed, this approach predicts a S/B ratio similar to 13. The background 
rejection power is enhanced with respect to the sequential cuts, since correlations among variables 
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Figure 6.5: Multi-dimensional likelihood distributions for i/f, CC and r — > e events. The last bin 
in signal includes the event overflow. Error bars in VeCC + VtCC sample represent statistical 
fluctuations in the expected proflle measurements after 5 years of data taking with shared running 
CNGS and a 5 T600 detector conflguration. 



are taken into account. In Fig. 16.61 we show the distribution of the surviving signal and background 
events for InA > 1.8. The region selected by the multidimensional likelihood approach, in the three 
dimensional space of the chosen variables, docs not present sharp edges like in the case of sequential 
cuts: The topology of the selected region is more complex thanks to the use of variable correlations. 



The Fisher discriminant method 

The Fisher discriminant method [128] is a standard statistical procedure that, starting from a 
large number of input variables, allows us to obtain a single variable that will efficiently distinguish 
among different hypotheses. Like in the likelihood method, the Fisher discriminant will contain 
all the discrimination information. 

The Fisher approach tries to find a linear combination of the following kind 

n 
i=l 

of an initial set of variables {xj} which maximizes 
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Figure 6.6: Distribution of signal and background events in EyisiUe, Pt^, Pi and P™*"** after the 
cut in In A. We see how signal appearance is considerably enhanced respect to background in 
regions defined by the InA > 1.8 cut. 



J{{aA) = 7 _T\ (6-3) 

y'^sig "bkg) 

where t is the mean of the t variable and a its variance. This last expression is nothing but a 
measurement, for the variable t, of how well separated signal and background are. Thus, by maxi- 
mizing Eq. ()6.3p we find the optimal linear combination of initial variables that best discriminates 
signal from background. The parameters aj which maximize Eq. (|6.3p can be obtained analytically 
by (see [HH]) 

^^ = w^\^J^'^^^') (6.4) 

where and /i^'^^ are the mean in the variable Xj for signal and background respectively, and 
W = Vsig + Vbkg- Vsig. back, the covariaucc matrices for the two event classes considered, are defined 
as: 

V^^=j:{x~ti%{x-fi%/N (6.5) 

where k stands for signal or background classes, (xi, Xn) is a vector of data with means 
{fii, and N is the total number of events. 
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In A > 1.8 


25 


1.1 ± 0.2 


12.9 ± 0.5 


In A > 2.0 


23 


0.86 ±0.16 


11.7±0.5 


InA > 2.5 


16 


0.40 ±0.12 


8.1 ±0.4 


InA > 3.0 


10 


0.22 ±0.08 


5.2 ±0.3 


InA > 3.5 


7 


0.12 ±0.06 


3.3 ±0.2 



Table 6.4: Expected number j^eCC background and signal events in the t ^ e analysis. A multi- 
dimensional likelihood function is used as the unique discriminant. Numbers are normalized to 5 
years running of CNGS. Errors in the number of expected events are given by Poisson statistics. 

A Fisher function for r search 

From the distributions of kinematical variables for Vt-GC and I'eCC, we can immediately con- 
struct the Fisher function for a given set of variables. Initially we selected the same set of variables 
we used for the likelihood approach, namely: EyisiUe, P™*^" and p/. We need only the vector 
of means and covariance matrices in order to calculate the optimum Fisher variable (Eq. (|6.4|) ). 
Distributions are shown in Fig. 16.71 where the usual normalization has been assumed. In Tab. 16.51 
values for the expected number of signal and background events arc shown as a function of the cut 
on the Fisher discriminant. Since correlations among variables are taken into account, it is not 
surprising that the result for the Fisher discriminant method is similar to the one obtained using 
a multidimensional likelihood. 
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46 


6.9 ±0.3 


23.1 ±0.4 


Fishcr> 0.0 


33 


2.4 ±0.2 


16.6 ±0.4 


Fisher> -0.27 


25 


1.15 ±0.13 


12.9 ±0.3 


Fisher > -0.5 


20 


0.60 ±0.10 


10.2 ±0.3 


Fisher > -1.0 


10 


0.14 ±0.05 


5.2 ±0.2 



Table 6.5: Expected number i/eCC background and signal events in the t — > e analysis. A Fisher 
variable is used as the unique discriminant. Numbers are normalized to 5 years running of CNGS. 
Errors in the number of expected events are given by Poisson statistics. 
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Fischer discriminant separation capabilities 
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Figure 6.7: The Fischer discriminant variable. Error bars in ly^CC + j^rCC sample represent 
statistical fluctuations in the expected profile measurements after 5 years of data taking with 
shared running CNGS and a 5 T600 detector configuration 



Increasing the number of initial kinematical variables 

Contrary to what happens with a multi-dimensional likelihood (where the increase in the num- 
ber of discriminating variables demands more Monte-Carlo data and therefore it is an extreme 
CPU-consuming process), the application of the Fisher method to a larger number of kinematic 
variables is straightforward, since the main characteristic of the Fisher method is that the final 
discriminant can be obtained algebraically from the initial distributions of kinematic variables. 
For instance, a Fisher discriminant built out of 9 kinematic variables {Eyis, P™'*"*, pi, p!^^ , Eiep, 
Pm, Qt, mr, Q/epjl predicts for 12.9 ± 0.3 taus a background of 1.17 ± 0.14 i/g CC events. We 
conclude that, for the Fisher method, increasing the number of variables does not improve the 
discrimination power we got with the set Ems, P™^^^ , pi and therefore these three variables are 
enough to perform an efficient r appearance search. 

6.2.3 Kinematical search using neural networks 

In the context of signal vs background discrimination, neural networks arise as one of the most 
powerful tools. The crucial point that makes these algorithms so good is their ability to adapt 



^see [125j for a detailed explanation of the variables 
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themselves to the data by means of non-linear functions. 

Artificial Neural Networks techniques have become a promising approach to many computa- 
tional applications. It is a mature and well founded computational technique able to learn the 
natural behaviour of a given data set, in order to give future predictions or take decisions about 
the system that the data represent (see [5T] and [129] for a complete introduction to neural net- 
works). During the last decade, neural networks have been widely used to solve High Energy 
Physics problems (see |130j for a introduction to neural networks techniques and applications 
to HEP). Multilayer Perceptrons efficiently recognize signal features from, an a priori, dominant 
background environment f |131| . [132j ). 

We have evaluated neural network performance when looking for Vt oscillations. As in 

the case of a multidimensional likelihood, a single valued function will be the unique discriminant. 
This is obtained adjusting the free parameters of our neural network model by means of a training 
period. During this process, the neural network is taught to distinguish signal from background 
using a learning data sample. 

Two different neural networks models have been studied: the multi-layer perceptron and the 
learning vector quantization self-organized network. In the following, their foundations and the 
results obtained with both methods are discussed. 



The Multi-layer Perceptron 

The multi-layer perceptron (MLP) function has a topology based on different layers of neurons 
which connect input variables (the variables that define the problem, also called feature variables) 
with the output unit (see Fig. 16. 8p . The value (or "state") a neuron has, is a non- linear function 
of a weighted sum over the values of all neurons in the previous layer plus a constant called bias: 



cjI.s'-'+I^ (6.6) 



where s\ is the value of the neuron i in layer I; uj\^ is the weight associated to the link between 
neuron i in layer / and neuron j in the previous layer (/ — 1); 6' is a bias defined in each neuron and 
g{x) is a non-linear function called transfer function. The transfer function is used to regularize 
the neuron's output to a bounded value between and 1 (or -1,1). An usual choice for the transfer 
function is the sigmoid function (see Fig. I6.9p . 

9t{x) = — ^ = + (6-7) 

This function tries to mimic the response of biological neurons by putting the neuron in a 
high (~ 1) or low (~ 0) state depending on the output of the previous neurons and how the links 
(weights) are set. The temperature parameter T is not determinant in principle and is always set 
to 1. 

These kind of structures for a Neural Network, where the information is propagated from the 
input units to the output ones passing through different ordered layers, are called feed forward 
neural networks. 



Training the MLP 



In a multi-layer perceptron a non-linear function is used to obtain the discriminating variable. 
The MLP is fed with a set of input variables which define the problem. The way the MLP sets its 
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Layer (1-1) Layer 1 

Figure 6.8: A general Multi-Layer Perceptron diagram. The optimal non-linear function of input 
variables (x,) is constructed using a set of basic units called neurons. Each neuron has two free 
parameters that must be adjusted minimizing an error function. 

free parameters is through minimization of an error function. This error function is defined as: 

E= J2 J2 {o,-d,f (6.8) 

eeSample eeSample 

where Oe represents the MLP output for a given event belonging to the training sample, and de is 
the "desired output" . The previous expression is nothing but a distance between the current MLP 
performance and the ideal case. In pattern recognition tasks, we would like to have some kind of 
discrimination function which, in an ideal case, always gives the correct answer, i.e., 1 for signal 
events and for background events. 

Once a topology is chosen for the MLP (number of layers and neurons per layer) , we need a first 
sample of signal and background events for the training stage. The weights and the biases are then 
adapted, minimizing the error function (j6.8|) . following the standard back-propagation algorithm 
(see |81| or [129] ). which is the general procedure for a general feed- forward neural network and 
an arbitrary differentiable activation function. The process of adapting the parameters is repeated 
over the training sample during a number of "epochs" until some condition is fulfilled (usually 
when a stable minimum is reached). 

At the end of the training stage, the network is made, and final weights and biases of the network 
are the ones which provide the best response for the training sample. In order to determine the 
real discrimination capabilities of the MLP, we need to test the neural network performance over 
an statistically independent sample (to avoid biases). 

Freedom of choice 

The construction of a MLP implies that several choices must be made a priori: 
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Figure 6.9: Sigmoid function. This is the usual choice for the transfer function of neurons in the 
MLP modeL It has the abihty of regularizing the initial input variables to a value between and 
1, imitating, from the point of view of biological models, the response of real neurons. 

• Number of input variables: The MLP can be fed, in principle, with as many variables 
as we consider fit in order to fully define our pattern recognition problem. We need to know 
a priori which variables (or combination of them) represent the main features of the classes 
in our pattern selection task. Redundant variables contribute to a worse performance of the 
MLP, hence is better to reduce the number of input variables to avoid a loss of discrimination 
power. 

• Number of hidden layers: It depends on the specific task one wants the network to 
perform, but a general statement is that no more than two hidden layers are needed. In 
pattern classification tasks, regions in input space must be created depending on the class 
they belongs. It could be shown (see [81]) that a MLP with no hidden layers only can made 
decision boundaries consisting in a hyperplane, one hidden layer can perform a single convex 
region consisting in hyperplanes segments, finally, two hidden layers could generate arbitrary 
decision regions. Hence, if your classes are very mixed in the variable space, two hidden 
layers will be able to construct more complex region borders. 

• Number of neurons per layer: The view of pattern classification as a "cut out" procedure 
gives a clue to the minimum number of hidden units needed for a specific classification task. 
In the "hard" case, where the first hidden layer has to define the border of a closed volume in 
the input space, at least (N+1) units (where N is the number of input variables) are needed. 
The second hidden layer usually performs some kind of logical function on the units belonging 
to the first hidden layer, therefore, a small number of neurons (1 or 2) is enoug iB. 

^Logical function as AND or OR can be exactly performed using a perceptron with no more than one hidden 
layer with two neurons. See |130j for an example. 
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• Number of events in each sample: The size of the simulated data set is crucial in order 
to optimize the training algorithm performance. If the training sample is small, it is likely for 
the MLP to adjust itself extremely well to this particular data set, thus losing generalization 
power (when this occurs the MLP is over-learning the data). Over-learning can be tested 
on-line by defining a statistical independent sample to compute the error function in each 
step of the learning procedure. 

• Number of epochs: The number of iterations in the training procedure depends on how low 
and stable is the minimum we get. In addition over-learning must be avoided: by definition, 
the error curve in the learning sample is always macroscopically decreasing, however the error 
curve in test sample could start to rise. At this moment the MLP loses generalization power. 
Hence, we must stop our learning phase before over-learning starts. 

Multi-layer Perceptron setup 

As already mentioned in 16.2.21 five variables (three in the transverse plane and two in longi- 
tudinal direction) utterly describe the event kinematics, provided that we ignore the jet structure; 
therefore, using five independent variables we fully define our tagging problem. Initially we build 
a MLP that contains only three input variables, and in a latter step we incorporate more variables 
to sec how the discrimination power varies. The three chosen variables are Evisihie^ P^^^^ and p/. 
Our election is similar to the one used for the multidimensional likelihood approach, this allows to 
make a direct comparison of the sensitivities provided by the two methods. 

There are plenty of packages available on the market which implement a multi-layer perceptron. 
The simplest choice for us was to use the MLPfit package |133| . interfaced in PAW. 

We have opted for a neural network topology consisting on two hidden layers with four neurons 
in the first hidden layer and one in the second (see Fig. 16.10")) . 

Simulated data was divided in three, statistically independent, subsets of 5000 events each 
(2500 signal events and a similar amount of background). 

The MLP was trained with a first "training" data sample. Likewise, the second "test" data set 
was used as a training sample to check that over-learning does not occur. Once the MLP is set, the 
evaluation of final efficiencies is done using the third independent data sample (namely, a factor 40 
(75) larger than what is expected for background (signal) after five years of CNGS running with a 
T3000 detector). 

Error curves (Eq. (j6.8p ) during learning arc shown in Fig. 16.111 for training and test samples. 
We see that even after 450 epochs, over-learning does not take place. Final distributions in the 
multi-layer perceptron discriminating variable can be seen in Fig. 16.121 

Results of MLP analysis 

Fig. 16.131 shows the number of signal and background expected after 5 years of data taking 
as a function of the cut in the MLP variable. In Fig. 16.141 we represent the probability of an event, 
falling in a region of the input space characterized by MLP output > cut to be a signal event 
(top plot), and the statistical significance as a function of the MLP cut (bottom plot). Like in the 
case of the likelihood, a cut based on the MLP output variable can select regions of complicated 
topology in the kinematical hyperspace (see Fig. 16.15)) . 

Selecting MLP > 0.91 (overall r selection efficiency — 25%), the probability that an event 
falling in this region is signal amounts to ^ 0.95. For 5 years of running CNGS and a T3000 
detector, we expect a total amount of 12.9 ±0.5 VrCC (r e) events and 0.66 ±0.14 i^eCC events. 
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Figure 6.10: The MLP setup. Chosen topology for the MLP. We feed a two layered MLP (4 
neurons in first layer and 1 in second) with input variables: EmsiUe, P™'** and pi. 



Tab. 16.61 summarizes as a function of the applied MLP cut the expected number of signal and 
background events. 

If we compare the outcome of this approach with the one obtained in section [6. 2. 21 wc see that 
for the same r selection efficiency, the multidimensional likelihood expects 1.1 ± 0.2 background 
events. Therefore, for this particular cut, the MLP achieves a 60% reduction in the number of 
expected CC background events. 



Increasing the number of feature variables 



Contrary to the likelihood approach, where the lack of very large Montc-Carlo samples, in- 
creases the complexity of the problem as the number of input variables grows (even for three input 
variables, smoothing algorithms have to be applied in order to get smooth distributions throughout 
the hyper-space), the MLP approach shows no technical difficulties in what concerns an increase 
of the input (feature) variables. 

Given that five variables define the problem under study, we expect, a priori, an enhancement 
of the sensitivity when we move from three to five feature variables. Notwithstanding the increase 
on the number of inputs, we observed that the change in the final sensitivity is negligible. The 
discrimination power is provided by EmsiUe, P™'^* and pi. The surviving background can not 
be further reduced by increasing the dimensionality of the problem. The application of statistical 
techniques able to find complex correlations among the input variables is the only way to enhance 
background rejection capabilities. In this respect, neural networks are an optimal approach. 
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Figure 6.11: Learning curves for the MLP. The neural network is trained for 450 epochs in order 
to reach a stable minimum. The solid line represents the error on training sample, the dashed line 
is the error on the test sample. Both lines run almost parallel: no over-learning occurs. 



Extended fisher discriminant method 

Given that an increase on the number of input variables does not improve the discrimination 
power of the Multi-layer perceptron, we tried to enhance signal efRciency following a different ap- 
proach: optimizing the set of input variables by finding new linear combinations of the original 
ones (or functions of them like squares, cubes, etc). 

To this purpose, using the fast computation capabilities of the Fisher method, we can operate 
in a systematic way in order to find the most relevant feature variables. Starting from an initial set 
of input variables, the algorithm described in |134| tries to gather all the discriminant information 
in an smaller set of optimized variables. These last variables are nothing but successive Fisher 
functions of different combinations of the original ones (see |134j ). In order to allow not only linear 
transformations, we can add non-linear functions of the kinematical variables like independent 
elements of the initial set. 

We performed the analysis described in |134j with 5 initial kinematical variables {E„is , P™*'** , 
pi, P}^^ and Eiep) plus their cubes and their exponentials (in total 15 initial variables). At the 
end, we chose a smaller subset of six optimized Fisher functions that we use like input features 
variables for a new Multi-layer perceptron. 

The MLP analysis with six Fisher variables does not enhance the oscillation search sensitivity 
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MLP variable. T3000 detector, 5 years CNGS 
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Figure 6.12: Multi-layer perceptron output for VrCC (r e) and z^gCC events. We see how signal 
events accumulate around 1 while background peaks at 0. 



Cuts 


i^^CC (r e) 




i^r CC (r ^ e) 




Efficiency 








(%) 




3 X 10-3 


Initial 


100 


252 


50 


Fiducial volume 


65 


164 


33 


MLP > 0.70 


42 


4.0 ±0.4 


21.4±0.6 


MLP > 0.75 


40 


3.0 ±0.3 


19.9 ±0.6 


MLP > 0.80 


37 


2.1 ±0.3 


18.6 ±0.5 


MLP > 0.85 


33 


1.5 ±0.2 


16.4±0.5 


MLP > 0.90 


27 


0.76 ±0.15 


13.5 ±0.5 


MLP > 0.91 


25 


0.66 ±0.14 


12.9 ±0.5 


MLP > 0.95 


19 


0.28 ±0.09 


9.6 ±0.4 


MLP > 0.98 


12 


0.09 ±0.05 


5.8 ±0.3 



Table 6.6: Expected number of background and signal events when a multi-layer perceptron func- 
tion is used as the unique discriminant. Numbers are normalized to 5 years running of CNGS. 
Errors in the number of events expected are of statistical nature. 
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Figure 6.13: Number of signal and background events after 5 years of running CNGS as a function 
of the MLP cut. Shadowed zones correspond to statistical errors. 



that we got with the three usual variables Eyis, P™'**and p;. We therefore conclude that neither 
the increase on the number of features variables nor the use of optimized linear combinations of 
kinematical variables as input, enhances the sensitivity provided by the MLP. 

Self organizing neural networks: LVQ network 

A self- organizing (SO) network operates in a different way than a multilayer perceptron does. 
These networks have the ability to organize themselves according to the "natural structure" of the 
data. They can learn to detect regularities and correlations in their input and adapt their future 
response to that input accordingly. A SO network usually has, besides the input, only one layer 
of neurons that is called competitive layer. 

Competitive self-organization 

Neurons in the competitive layer are able to learn the structure of the data following a simple 
scheme called competitive self- organization (see [130j ). which "moves" the basic units (neurons) in 
the competitive layer in such a way that they imitate the natural structure of the data. It works 
as follows: 

Suppose we have a set of N input variables x^^^ = {x^^\x'^\ ...,a;^'') taking values over the 
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Figure 6.14: (Top) probability of an event belonging to a region in input variable space character- 
ized by MLP > cut of being a signal event. (Bottom) Statistical significance of signal events as a 
function of the cut in MLP. 



total training sample (superscript (p) denotes one element of the sample). These input nodes 
are connected with each of the M neurons in the competitive layer by a vector of weights Wi = 
{wii, Wi2, WiN) (see Fig. I6.16p . From an initial random distribution of the weights, competitive 
self-organization operates computing the following quantity in each neuron: 

where i is an index over the number of neurons in the competitive layer. Once the hi are computed, 
a "winner neuron" is selected: 

hm = min{hj) 

A winner neuron is the one which has its vector of weights closest to the pattern. The output 
of the network is a vector o = (oi, 02, oai) which returns 1 in the m component corresponding to 
the winner neuron and for all the others. The weight vector Wm belonging to the winner neuron 
is then moved closer to x^^-* according to 

Awm = r,(f(f) - w^) (6.9) 

where rj is the learning rate. Thus, the neuron whose weight vector is closest to the input vector 
is promoted to be even closer. The result is that the winner neuron is more likely to win the 




Figure 6.15: Kinematical variables before (left histograms) and after (right histograms) cut are 
applied based on the MLP output. We see how the MLP has learnt that signal events favour low 
Evisibie, high and low values. 



competition the next time a similar vector is presented, and less likely to win when a different 
input vector is presented. As more and more inputs are presented, each neuron in the layer closest 
to a group of vectors soon adjusts its weights vector toward those input vectors. Eventually, if 
there are enough neurons, every cluster of similar input vectors will have a neuron that outputs 
1 when a vector in the cluster is presented, while outputting a at all other times. Thus, the 
competitive network learns to categorize the input vector it sees. 
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Competitive Layer Output 




Figure 6.16: Schematic diagram of the general topology for a self-organized neural network. Neu- 
rons in the competitive layer are connected with each one of the input nodes. As output, it returns 
a vector indicating the winner neuron. 

Learning vector quantization networks 

Although competitive self-organization is an unsupervised learning algorithm, for classification 
purposes one can improve the algorithm with supervised learning in order to fine tune final positions 
of the neurons in the competitive layer. This is called learning vector quantization (LVQ). 

The modifications respect to the previous scheme are quite simple: First we have to associate 
neurons in the competitive layer with the different classes we try to classify, i.e. if we use 8 neurons, 
4 can be associated to signal class and 4 to the background class. In addition to the input vector 
x^P^ we need to provide information about the class it belongs, i.e. signal or background event. 
Now, for a given input vector belonging to a certain class, a winner neuron is selected that could 
belong to the same class than the input vector or not. If the classes match, the vector weight of 
the winner neuron is then adapted according to Eq. (|6.9p : if they mismatch the modified adapting 
rule is simple 

AWm = - Wm) 

As a result, neurons belonging to a given class are favoured to adapt themselves to the regions 
where the patterns of their parent class lie, and all the others tend to go away and trying to find 
their corresponding regions, as well. 

At the end of training period, all the vector weights are placed in an optimized way around 
the input variable space. Afterward, we can test the neural network performance with a statistical 
independent sample. When the LVQ is fed with new patterns, it selects the closest neuron to each 
pattern and, depending on the class label of that neuron, the pattern is classified. An important 
difference with respect to the multi-layer perceptron approach is that in LVQ we always get a 
discrete classification, namely, an event is always classified in one of the classes. The only thing 
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you can estimate is the degree of belief in the LVQ choice. 
A LVQ network for r search 

Once again we use EmsMe, P^"'"* and pi as discriminating variables. A LVQ network with 
10 neurons has been trained with samples of 2500 events for both signal and background. Given 
that, before any cuts, a larger background sample is expected, we have chosen an asymmetric con- 
figuration for the competitive layer. Out of 10 neurons, 6 were assigned to recognize background 
events, and the rest were associated to the signal class. After the neurons are placed by the training 
procedure, the LVQ network is fed with a larger data sample (statistically independent) of 6000 
signal and background events. The output provided by the network is plotted in Fig. 16.171 We 
see how events are classified in two independent classes: signal like events (labeled with 2) and 
background like ones (labeled with 1). 68% of I'rCC (r e) events and 10% of i^eCC events, 
occurring in fiducial volume, are classified to be a signal like event. This means a r efficiency 
around 45% with respect to the tau events generated in active LAr. For the same r efiiciency, the 
multi-layer perceptron only misclassifies around 8% of VeGC events. 

Several additional tests have been performed with LVQ networks, by increasing the number of 
input variables and/or the number of neurons in the competitive layer. No improvement on the 
separation capabilities have been seen. For instance, a topology with 16 feature neurons in the 
competitive layer and 4 input variables (we add the transverse lepton momentum) leads to exactly 
the same result. 

The simple geometrical interpretation of this kind of neural networks supports our statement 
that the addition of new variables to the original set {Eyisibie, P™'*^, pi} does not enhance the 
discrimination power: the bulk of signal and background events are not better separated when we 
increment the dimensionality of the input space. 

Combining MLP with LVQ 

We have seen that LVQ networks returns a discrete output. The whole event sample is classified 
by the LVQ in two classes: signal-like and background-like. We can use the classification of a 
LVQ as a pre-classification for the MLP. A priori, it seems reasonable to expect an increase on the 
oscillation search sensitivity if we combine the LVQ and MLP approaches. The aim is to evaluate 
how much additional background rejection, from the contamination inside the signal-like sample, 
can be obtained by means of a MLP. 

We present in Fig. 16.181 the MLP output for events classified as signal-like by the LVQ network 
(see Fig. I6.17|) . Applying a cut on the MLP output such that we get 12.9 ± 0.5 signal events (our 
usual reference point of 25% r selection efficiency), we get 0.82±0.19 background events, similar to 
what was obtained with the MLP approach alone. This outcome conclusively shows that, contrary 
to our a priori expectations, an event pre-classification, by means of a learning vector quantization 
neural network, does not help improving the discrimination capabilities of a multi-layer perceptron. 
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Figure 6.17: LVQ neural network separation capabilities. In competitive self-organized networks a 
discrete decision is always issued: signal like events are labeled with 2 and background like with 1 . 

6.3 Ur appearance discovery potential 

We have studied several pattern recognition techniques applied to the particular problem of search- 
ing for i/^ Vt oscillations. Based on discovery criteria, similar to the ones proposed in |135| for 
statistical studies of prospective nature, we try to quantify how much the discovery potential varies 
depending on the statistical method used. 

We define fis and /is as the average number of expected signal and background events, respec- 
tively. With this notation, the discovery criteria are defined as: 

1. We require that the probability for a background fluctuation, giving a number of events equal 
or larger than fis + fJ-B, be smaller than e (where e is 5.733 x 10^^, the usual ba criteria 
applied for Gaussian distributions). 

2. Wc also set at which confidence level (1 — S), the distribution of the total number of events 
with mean value fJ.s + fJ-B fulfills the background fluctuation criteria stated above. 

For instance, if S is 0.10 and e is 5.733 x 10~^, we are imposing that 90% of the times we repeat 
this experiment, we will observe a number of events which is 5(t or more above the background 
expectation. 

For all the statistical techniques used, we fix (5^:0.10 and e=5.733 x 10~^. In this way we can 
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Figure 6.18: LVQ and MLP networks combined. Distributions are given in the continuous MLP 
variable. Only events labeled by LVQ with 2 {signal like) have been used for the analysis. 



compute the minimum number of events needed to establish that, in our particular example, a 
direct oscillation has been discovered at the CNGS beam. 

In Tab. [677] we compare the number of signal and background events obtained for the multi-layer 
perceptron and the multi-dimensional likelihood approaches after 5 years of data taking with a 3 
kton detector. We also compare the minimum exposure needed in order to claim a discovery. The 
minimum exposure is expressed in terms of a scale factor a, where a = 1 means a total exposure 
of 11.75 ktonxyear. For the multi-layer perceptron approach (a = 0.86), a statistically significant 
signal can be obtained after a bit more than four years of data taking. On the other hand, the 
multi-dimensional likelihood approach requires 5 full years of data taking to claim a discovery. 
Therefore, when applied to the physics quest for neutrino oscillations, neural network techniques 
are more performing than classic statistical methods. 

Conclusions 

We have considered the general problem of — > Vr oscillation search based on kinematic criteria 
to assess the performance of several statistical pattern recognition methods. 
Two are the main conclusions of this study: 

• An optimal discrimination power is obtained using only the following variables: EmsiUe, 
pmiss g^j^^ their correlations. Increasing the number of variables (or combinations 
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Multi-layer 


Multi-dimensional 




Perceptron 


Likelihood 


# Signal 


12.9 


12.9 


# Background 


0.66 


1.1 


a factor 


0.86 


1.01 



Table 6.7: Number of signal and background events for the multi-layer perceptron and the multi- 
dimensional likelihood approaches. Numbers are normalized to 5 years of data taking in shared 
CNGS running mode and a 3 kton detector configuration. The last row displays the scale factor a 
needed to compute the minimum exposure fulfilling the discovery criteria described in the text. 

of variables) only increases the complexity of the problem, but does not result in a sensible 
change of the expected sensitivity. 

• Among the set of statistical methods considered, the multi-layer perceptron offers the best 
performance. 

As an example, we have considered the case of the CNGS beam and i^r appearance search (for 
the T ^ e decay channel) using a very massive (3 kton) Liquid Argon TPC detector. Fig. 16.191 
compares the discrimination capabilities of multi-dimensional likelihood and multi-layer perceptron 
approaches. We see that, for the low background region, the multi-layer perceptron gives the best 
sensitivity. For instance, choosing a r selection efficiency of 25% as a reference value, we expect 
a total of 12.9 ± 0.5 VtCC (t e) signal and 0.66 ± 0.14 v^CC background. Compared to 
multi-dimensional likelihood predictions, this means a 60% reduction on the number of expected 
background events. Hence, using a multi-layer perceptron, four years of data taking will suffice to 
get a statistically significant signal, while five years are needed when the search approach is based 
on a multi-dimensional likelihood. 
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Figure 6.19: Multi-layer perceptron vs multi-dimensional likelihood. We assume 5 T600 module 
configuration and 5 years of CNGS running (7.6 Kton/year exposure). The shadowed area shows 
the statistical error. 



Conclusions 



The results of this dissertation involve different aspects of Neutrino Physics. After a brief review 
of the main theoretical topics on neutrino interactions and the problem of their masses, a phe- 
nomenological study of the possibilities offered by future lepton colliders for the exploration of the 
limits on the mixing of light and new heavy neutrinos is presented. Afterwards we have focused 
our interest on experimental Neutrino Physics with LAr detectors, which is the main issue of this 
thesis. A sound and innovative work has been carried out in the development of the LAr technique 
and the necessary tools to perform Physics analyses, thus establishing the high potentiality of these 
detectors for the study of Neutrino Physics. 

Assuming the possibility of the existence of heavy neutrinos with masses of few hundreds of 
GeV, we have analysed the contribution of these neutrinos to electroweak process in future lepton 
colliders, exploring the limits on heavy neutrino mixing and masses. We have studied the ILC 
potential for their detection in the process e+e~ tWv, taking into account the SM background 
and the effects of ISR and breamstrahlung, paying special attention to the relevance of the final 
state lepton flavour and initial beam polarization. Using a parton simulation it has been shown 
that it is possible to observe a heavy neutrino signal in this final state if it has a 
mixing with the electron VeN ^ 10^^, and a mass below 400 GeV. Although a mixing 
with the electron of this size is necessary to observe a heavy neutrino at ILC, the signal may be 
more visible in the muon or tau channel if it also has a relatively large coupling to 
them. The production cross sections and then the discovery limits do not depend on the Dirac or 
Majorana nature of the heavy neutrino. These non-decoupled heavy neutrinos do not explain the 
observed light neutrino masses nor the baryon asymmetry of the universe. In this sense this search 
for heavy neutrinos at large lepton colliders is complementary to the joint experimental effort for 
determining the light neutrino properties, and in particular the neutrino mixing matrix (Ref [T]). 

Thereafter, in the bulk of the dissertation, we have been mainly concentrated on the study of the 
Physics potential of LAr TPCs. We have developed a set of tools and procedures which 
makes possible the calorimetric and 3D reconstruction of ionizing events occurring 
within the chamber. These tools have shown an excellent performance working on events of 
low-multiplicity composed of several well defined tracks. 

The technical operation of a 50 liter LAr TPC prototype, exposed to a multi-GcV i/^ neutrino 
beam, was rather good: The LAr purity was kept stable at a high-level during the whole data 
taking period and the electronic read-out configuration supplied an unbeatable signal-background 
ratio with no saturation of the electronic response, allowing an excellent imaging of the ionizing 
events. 

The 50L TPC offered the opportunity to measure and reconstruct for first time 
beam neutrino interactions using the LAr technique. Due to the smallness of the prototype, 
a joint operation with the NOMAD detector acting as a muon spectrometer downstream of the 
TPC has been necessary in order to measure the high energy muons resulting from the neutrino 
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interactions. Before the data analyses, a careful calibration of the TPC was done. All the needed 
parameters (the electron drift velocity, the alignment angles with respect to the beam, the electron 
lifetime and the recombination factors) were successfully obtained. 

During its operation, the 50L TPC gathered a substantial sample of neutrino interactions in the 
multi-GeV region, with a sizable amount of low-multiplicity events. From the whole sample of 
CC events, we have selected a set of events called the golden sample. They are events having a clear 
two prong topology with a muon leaving the chamber and an identified contained proton. This 
sample contains 86 successfully reconstructed events and, in spite of its limited size, 
they have made possible to study the accuracy of our current Monte-Carlo models for 
the interactions of neutrinos with nuclear matter. We expect that most of these events come 
from quasi-clastic neutrino interactions, but there is also a sizable background coming from deep- 
inelastic and resonant events. Thanks to the superb reconstruction capabilities of the LAr TPC, 
the kinematical distributions of the events have been obtained with high resolution. This allows to 
discriminate between models with a full description of nuclear effects and naive approaches that 
only include Fermi motion and Pauli blocking (Ref. [2]). 

On the other hand, the collected sample of golden events gives us the opportunity of measur- 
ing for first time the quasi-elastic CC cross section using a LAr TPC. To this purpose 
a precise knowledge of the beam composition, the beam profile, the exposition time, trigger effi- 
ciencies, acceptance cuts, event selection rates and background ratios was needed. In addition, a 
careful evaluation of the systematic errors entering in the computation has been done exploiting a 
full Monte-Carlo simulation of the neutrino beam, the nuclear effects and the experimental set-up. 
Our measurement of the total QE charged-current cross section (Ref. [3]) is 



C7QE = (0.90 ± 0.10 (stat) ± 0.18 (sys)) x 10"^^ cm^ 



This value agrees with previous measurements of this reaction reported in the literature, and it 
stands as the first Physics measurement ever done with a LAr TPC using a neutrino beam from 
an accelerator. 

Finally, given the excellent reconstruction capabilities of the LAr TPC technique and the ac- 
curacy in the reconstruction of the kinematic variables in the plane transverse to the incoming 
neutrino, we have evaluated our capability to perform a. i/^ ^ Vt- oscillation search 
based on kinematical criteria. A systematic assessment of the performance of several statisti- 
cal pattern recognition methods shows that an optimal discrimination power is obtained using only 
three variables and their correlations: The visible energy {Ems), the transverse missed momentum 
f^pmiss^ and the fraction of transverse momentum carried by the prong Icpton (pi). Moreover, 
among the set of statistical methods considered, the multi-layer pcrceptron offers the best perfor- 
mance. Considering the case of the i^r appearance search at the CNGS neutrino 
beam using a very massive (3 kton) Liquid Argon TPC detector, we expect a total of 
12.9 ± 0.5 Vt-CC (r e) signal and 0.66 ±0.14 z^gCC background after 5 years running. 
This means a 60% reduction of the number of expected background events with respect to other 
broadly use methods in Particle Physics like for example, multi-dimensional likelihoods (Ref. [i]). 



Appendix A 

Electroweak interactions 



The standard electroweak model is based on the gauge symmetry group SU{2)l x U{1)y, and the 
corresponding gauge bosons W^^, i = 1,2,3, and B^. After spontaneus symmetry breaking and 
the usual redefinition of the gauge fields [S] , the Lagrangian which describes the electroweak gauge 
interactions for fcrmions reads: 



- E - 9\l')Z,i^, (A.l) 

i 

i 

where ipi — (^^l) is the usual SU{2) doublet of the i*'' fermion family, being T+ and the 

raising and lowering operators of the same SU{2) representation. The first line of Eq. (|A.1|) 
describes charge current (CC) interactions, while the second one rules the neutral currents (NC); 
the electromagnetic interactions are governed by the third term. 



A.l Charged Current Interactions 

In the SM left-handed fermions are put together in three different families of SU{2)l doublets, 
while the right-handed components are singlets. Thus, the emission of a charged vector boson 

causes a transition between the left-handed fermions within the weak isospin 1/2 doublet ^ 
These processes xpi'^ — > -04^ + and ipi^ ipi^ + are governed by the Lagrangian 

= - ^ Ul^,W+ + j,^, V- ) (A.2) 

with the vector field operator^ and the fermion current j^^ summed over the 3 lepton and 
quark doublets: 

3 3 

j^c = E^^''2^^ " + E^T^'a^^ " ^'^^^ (^-^^ 

t=l 1=1 



^W^t is the operator that annihilates a or creates a W 
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Here, the field operators t/^i have been replaced by the particle symbols {i^iji ^ leptons, Ui,di ^ 
quarks) for simplicity. The projection operator ^(1 — 7^) assures that only leptons with left-handed 
chirality take part in the process 

vir^-^h = vn:rhL- (a.4) 

For massive fermions the weak isospin eigenstates are generally a mixture of the mass eigen- 
states. The basic mechanism by means of which fermions adquire mass in the SM preserving the 
gauge symmetry is through Yukawa couplings with the Higgs boson [8]. It is always possible to 
choose a fermion basis where just the fermions of one of the isospin eigenstates are rotated re- 
spect to their mass eigenstates after symmetry breaking. Thus, without loss of generality, only 
one unitary (3 x 3)-matrix operating in the is needed in order to describe the relation between 
interaction and mass eigenstates: 

-u[tL2] (A.5) 

It could have been equivalently done with the ipj fermions. Therefore, quarks and lepton basis 
are usually chooscn so that charge leptons (e /i r) and up quarks {u c t) coincide with their mass 
eigenstates, while their opposite weak isospin partners are rotated respect to their mass eigenstates. 
The charged fermion current can now be rewritten as a sum over the observed fermions 

Jc^c = (^^^)f//V^(l-7') +[uct)^^\{l-l')U, (A.6) 

The quarks and charged leptons correspond to the mass eigenstates shown in Tab. 11.11 Eq. (|A.6|) 
essentially implies that charged current weak interactions cause transitions between quarks or 
leptons of the three generations. 



A. 1.1 Charged Leptonic Current 

The recently found neutrino oscillations imply that the neutrinos i^e, v^j, and Vt are not the mass 
eigenstates Vmi- This means that the mixing matrix Ui is not diagonal [9j. In fact, current limits 
indicate mixing angles which correspond to almost maximal mixing for z/e <-> and v^^ ^ Vr- Also, 
there might be more than three neutrinos taking part in the mixing (see Sec ll.2p . Because of the 
small neutrino masses and since they interact only weakly, the mass eigenstates are not yet directly 
observable, but have to be determined by means of oscillation rates and cosmological measurements. 
The experimental evidence in favor of oscillations of solar (KamLAND) and atmospheric (K2K) 
neutrinos can be accommodated in the SM with three massive neutrinos. The charged leptonic 
current is often also written in terms of the weak eigenstates: 



Jcc,l 



e • 



(A.7) 
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A. 1.2 Charged Quark Current 

The mixing in the quark sector is described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix: 

(Uud Uus Uub \ 
Ucd Ucs Ucb ■ (A.8) 
Utd Uts Utb I 

The squared matrix elements |C/^yP give the probabihty of the quark Ui to change into the quark 
dj. Conservation of probabihty requires the CKM matrix to be unitary. In principle, the individual 
matrix elements can all be measured in weak decays of the relevant quarks or in deep inelastic 
neutrino-nucleon scattering. The 90 % confidence limits on the magnitude of the matrix elements, 
as given in [9], are 

0.9739- 0.9751 0.221 - 0.227 0.0029 - 0.0045 \ 
0.221 -0.227 0.9730- 0.9744 0.039- 0.044 . (A.9) 
0.0048 - 0.014 0.037- 0.043 0.9990- 0.9992 / 

In contrast with the leptonic sector, the mixing in the quark sector is easily observable, since 
the mass differences are bigger and the quarks interact also via the strong force. 



A. 2 Neutral Current Interactions 

The Z boson, just as the photon, is a mixture of the neutral gauge field of the SU{2)l symmetry 
group (the other two weak bosons being and W~) and the gauge field of the U{'i) symmetry 
group. The form of the weak neutral current, therefore, also depends on whether or not the 
fermions carry charge. The corresponding Lagrangian is given by 

z cos 

where 9w is the Weinberg angle defined by cos 9\\r = mw /mz ■ The neutral current is summed 
over all fermions 

' _ (A.ll) 

= E '/'.7''(r^'(l - 7') - sin^ ew)^i 

i 

and 

„; (A. 12) 

9\ = T'\ 



where is the weak isospin of the left-handed fermion, i. e. ±1/2, and Q' is the charge of the 
fermion in units of e. 



Appendix B 



Calculation of the e^e i W^u 
process with heavy neutrinos 



Following the notation that we introduced in Sec. 11.41 for the 6x6 neutrino mixing matrix 
(Eq. (fT3^ ). we write down the extended (MNS) matrix [HI [HS] as V = (C/ V). This matrix 
parameterises the charged and neutral current gauge interactions. On the other hand, the neu- 
trino mass eigenstates have been renamed as in that Section in order to make explicit the distintion 
between the light (i/l) and the heavy (Nl) sectors. The electroweak Lagrangian is then written as 

£gc = -^Irf'V^;;^) »7+Hx., (B.l) 

In the evaluation of e+e~ i^W^iy, with —i- qq' , we will only consider the contributions 
from the diagrams in Fig. IB.li neglecting diagrams with four fermions e^qq'v in the final state 
but with qq' not resulting from a W decay. At any rate, we have checked that the corresponding 
contributions are negligible in the phase space region of interest. The discovery of a new heavy 
neutrino in e+e^ — + llWv requires its observation as a peak in the invariant mass distribution, 
otherwise the irreducible SM background is overwhelming. This requires to reconstruct the W , 
what justifies to consider tWv production (instead of general four fermion production), with W 
decaying hadronically. 

Let us discuss the contributions and sizes of the different diagrams for e^e'^ — > £~W~^i' in Fig. 
IB. II and what we can learn from this type of processes at ILC. The first four diagrams are SM 
contributions. Diagrams 5,7—9 are present within the SM, mediated by a light neutrino, but they 
can also involve a heavy one. Diagrams 6 and 10 arc exclusive to Majorana neutrino exchange. 
The SM contribution has a substantial part from resonant W^W~ production, diagrams 4 and 8, 
especially for final states with i = /i, t. The heavy neutrino signal is dominated by diagrams 5 
and 6 with N produced on its mass shell, because the Tjv enhancement of the amplitude partially 
cancels the mixing angle factor in the decay vertex, yielding the corresponding branching ratio. It 
must be remarked that the s-channel N production diagram 7 is negligible (few per mille) when 
compared to the t-, u-channel diagrams 5 and 6. This behaviour is general, because the s-channel 
propagator is fixed by the large collider energy and suppresses the contribution of this diagram, 
whereas the t- and w-channel propagators do not have such suppression. Since both diagrams 5 
and 6 involve an eNW vertex to produce a heavy neutrino, only in the presence of this interaction 
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8 9 10 



Figure B.l: Diagrams contributing to e+e — > e W'^v. For I = ^i,t only diagrams 3-8,10 con- 
tribute. 



the signal is observable. Once the heavy neutrino is produced, it can decay to H.W with £ ^ e, fj., r, 
being the corresponding branching ratios in the ratio I VeATp : |V^Arp : I^^tJvP- 

In case that the signal is dominated by t and u channel on-shell N production (the only situation 
in which it is observable) , negative electron polarisation and positive positron polarisation increase 
its statistical significance. For the signal contributions alone we have cr + - : cr + - = 1200 : 1 
(for mpf ~ 300 GeV, VeN ~ 0.073), with a^^- = a^+- ~ 0. For the SM process, : 
(T + - : cr + - = 150 : 7 : 1, cr + - ~ 0. In the limit of perfect beam polarisations -Pg- ~ 
Pg+ = 1, the signal is enhanced with respect to the background by a factor of 1.05 and, what 
is more important, the ratio S/^/~B increases by a factor of two. Using right-handed electrons 
and left-handed positrons decreases the S/B ratio by a factor of 8, and S/^fB by a factor of 50. 
On the other hand, if the neutrino does not mix with the electron but mixes with the muon or 
tau, the behaviour is the opposite. Since the only contribution comes from diagram 7 the use of 
left-handed positrons and right-handed electrons actually increases the signal, while reducing the 
SM cross section for this process |137j . This case is of limited practical interest, since for VeN = 
the signal is barely observable. 

We finally point out that the signal cross section exhibits little dependence on the heavy neutrino 
mass, except close to the kinematical limit |138[ 1139] . and the final results are almost independent 
of m^v within the range 200-400 GeV [13 . For our calculations we take mjv = 300 GcV. In 
contrast with what has been claimed in the literature |140| , we find equal production cross sections 
for Dirac and Majorana neutrinos to a very good approximation. The reason is easy to understand: 
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while in the present case the signal is strongly dominated by diagrams 5 and 6 (which give equal 
contributions to the cross section and do not interfere because light neutrino masses can be safely 
neglected), for a Dirac neutrino only diagram 5 is present. On the other hand, the width of a Dirac 
neutrino is one half of the width of a Majorana neutrino with the same mixing angles |43| . 



Generation of signals 

The matrix elements for e+e^ ?~W'^v l~qq'v are calculated using HELAS |141| . including all 
spin correlations and finite width effects. We sum SM and heavy neutrino-mediated diagrams at 
the amplitude level. The charge conjugate process is included in all our results unless otherwise 
noted. We assume a CM energy of 500 GeV, with electron polarisation P^,- = —0.8 and positron 
polarisation Pg+ = 0.6. The luminosity is taken as 345 fb~^ per year [142] . In our calculations we 
take into account the effects of ISR fl43j and beamstrahlung |1441 1145j . For the design luminosity 
at 500 GeV we use the parameters T = 0.05, N — 1.56 |142j . The actual expressions for ISR and 
beamstrahlung used in our calculation arc collected in Ref. |146j . We also include a beam energy 
spread of 1%. 

In final states with r Icptons, wc select r decays to tt, p and ai mesons (with a combined 
branching fraction of 55% [S]), in which a single v.,- is produced, discarding other hadronic and 
Icptonic decays. Wc simulate the r decay assuming that the meson and r momenta are coUinear 
(what is a good approximation for high r energies) and assigning a random fraction x of the r 
momentum to the meson, according to the probability distributions |147j 

P{x) = 2(1 - x) (B.3) 

for pions, and 

= 2C3-!c^ + i " (^-^^ 

for p and oi mesons, where C = rn^p a^l'^^T- We assume a r jet tagging efficiency of 50%. 

We simulate the calorimeter and tracking resolution of the detector by performing a Gaussian 
smearing of the energies of electrons, muons and jets, using the specifications in Ref. |148j . 

^ = i^®l%, ^=0.02%i.^ ^.^04%, (B.5) 
E"" ^/W Ei^ El y/El 

respectively, where the two terms are added in quadrature and the energies are in GeV. We 
apply kinematical cuts on transverse momenta, px > 10 GeV, and pseudorapidities \r]\ < 2.5, 
the latter corresponding to polar angles 10° < 6* < 170°. To ensure high t momenta (so that 
the meson resulting from its decay is effectively collinear) we require px > 30 GeV for t jets. 
We reject events in which the leptons or jets are not isolated, requiring a "lego-plot" separation 
Ai? = /S.rf -I- > 0.4. For the Monte Carlo integration in 6-body phase space we use RAMBO 

In final states with electrons and muons the light neutrino momentum pi, is determined from 
the missing transverse and longitudinal momentum of the event and the requirement that = 
(despite ISR and beamstrahlung. the missing longitudinal momentum approximates with a reason- 
able accuracy the original neutrino momentum). In final states with r Icptons, the reconstruction 
is more involved, due to the additional neutrino from the t decay. We determine the "first" neu- 
trino momentum and the fraction x of the r momentum retained by the r jet using the kinematical 
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constraints 



E^/ + Ey H — Ej 




(B.6) 



in obvious notation. These constraints only hold if ISR and beamstrahlung are ignored, and in 
the limit of perfect detector resolution. When solving them for the generated Monte Carlo events 
we sometimes obtain x > 1 or x < 0. In the first case we arbitrarily set a; = 1, and in the second 
case we set x = 0.55, which is the average momentum fraction of the r jets. With the procedure 
outlined here, the reconstructed r momentum reproduces with a fair accuracy the original one, 
while the obtained Pi, is often quite different from its original value. 
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